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total  pressure,  static  pressure,  and  the  three  components  of  velocity  to  b* 
determined.  These  data  were  then  Fourier  analyzed  to  provide  a  method  of 
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the  distortion  to  be  evaluated.  A  portion  of  the  data  obtained  in  this 
study  was  analyzed  to  determine  the  total  pressure  losses  associated  with 
a  distorted  inflow. 

The  results  of  this  investigation  revealed  that  the  attenuation  of  the 
distortion  is  a  function  of  the  blade  stagger  angle,  the  distortion  reduce* 
frequency,  and  the  intrablade  frequency.  Little  effect  of  rotor  steady 
blade  loading  or  of  rotor-stator  spacing  was  observed.  It  was  shown  that 
the  attenuation  of  the  distortion  was  maximized  at  values  of  rotor  blade 
spacing  corresponding  to  one-half  the  wavelength  of  the  distortion.  The 
amount  of  attenuation  measured  experimentally  was  found  to  be  much  greater 
than  that  predicted  by  an  unsteady,  two-dimensional,  inviscid  theory. 
However,  the  trends  of  the  predicted  and  measured  attenuation  with  reduced 
frequency  or  intrablade  frequency  were  similar.  This  discrepancy  is 
attributed  to  blade  thickness  and  viscous  effects  which  were  neglected  in 
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blade  was  significant. 
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ABSTRACT 


An  investigation  was  conducted  to  study  the  unsteady  response  of  an 
axial  flow  turbomachinery  rotor  to  spatial  variations  in  the  rotor  inflow 
velocity. ""'Distorted  inlet  flow  is  a  very  realistic  and  prevalent  problem 
in  jet  aircraft  engines,  and  the  consequences  of  subjecting  the  rotating 
blades  of  an  axial  flow  fan  or  compressor  to  an  inlet  distortion  are 
excessive  radiated  noise,  blade  vibration,  and  performance  degradation. 
Therefore,  it  is  important  for  the  turbomachinery  designer  to  be  con¬ 
cerned  with  designs  which  minimize  the  unsteady  response  of  turbomachines 
to  such  spatial  inlet  distortions.  The  study  discussed  in  this  report 
is  a  combination  of  experimental  and  analytical  efforts  which  are  directed 
toward  obtaining  information  relating  the  influence  of  distortion  charac¬ 
teristics  and  rotor  design  and  operating  variables  such  as  solidity, 
stagger  angle,  rotor-stator  spacing,  and  blade  loading  to  the  unsteady 
rotor  response.  The  approach  taken  to  measure  the  response  of  the  rotor 
was  to  make  detailed  circumferential  surveys  of  the  time-mean  flow  field 
upstream  and  downstream  of  the  rotor  using  five-hole  prism  type  probes. 
This  enabled  the  total  pressure,  static  pressure,  and  the  three  components 
of  velocity  to  be  determined.  These  data  were  then  Fourier  analyzed  to 
provide  a  method  of  quantifying  the  size  of  the  distortion,  thus  per¬ 
mitting  the  attenuation  of  the  distortion  to  be  evaluated.  A  portion  of 
the  data  obtained  in  this  study  was  analyzed  to  determine  the  total 
pressure  losses  associated  with  a  distorted  inflow.  • 

The  results  of  this  investigation  revealed  that  the  attenuation  of 


the  distortion  is  a  function  of  the  blade  stagger  angle,  the  distortion 
reduced  frequency,  and  the  intrablade  frequency.  Little  effect  of  rotor 
steady  blade  loading  or  of  rotor-stator  spacing  was  observed.  It  was 
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shown  that  the  attenuation  of  the  distortion  was  maximized  at  values  of 
rotor  blade  spacing  corresponding  to  one-half  the  wavelength  of  the 
distortion.  The  amount  of  attenuation  measured  experimentally  was  found 
to  be  much  greater  than  that  predicted  by  an  unsteady,  two-dimensional, 
inviscid  theory.  However,  the  trends  of  the  predicted  and  measured 
attenuation  with  reduced  frequency  or  intrablade  frequency  were  similar. 
This  discrepancy  is  attributed  to  blade  thickness  and  viscous  effects 
which  were  neglected  in  the  theory.  In  most  cases,  the  unsteady  total 
pressure  losses  were  found  to  be  greater  than  the  losses  in  a  uniform 
flow,  and  the  effect  of  measuring  the  losses  at  radial  locations  other 
than  the  mean  radius  of  the  blade  was  significant. 
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CHAPTER  I 


INTRODUCTION 

1. 1  Significance  of  Unsteady  Flow  in  Turbomachinery 

In  designing  the  blading  of  a  compressor  or  turbine,  the  air  flow  is 
assumed  to  be  steady.  The  existence  of  a  uniform,  steady  flow  is  unusual, 
however,  due  to  the  presence  of  viscous  or  real  fluid  effects.  In  a  real 
flow,  for  example,  viscous  effects  will  result  in  wall  boundary  layers, 
wakes  behind  solid  surfaces,  and  turbulence.  As  a  result,  the  forces  on 
a  blade  will  vary  with  time  as  the  flow  field  over  the  blade  changes. 
Development  of  sophisticated  Instrumentation  techniques  has  confirmed  the 
presence  of  strong  unsteady  effects  which  can  affect  the  aerodynamic  per¬ 
formance,  the  aeroelastic  performance,  and  noise  generation  in  turbo¬ 
machines.  Very  successful  turbomachines  have  been  developed  in  the  past 
by  compensating  for  the  Inadequate  understanding  of  unsteady  flows  with 
extensive  empirical  correlations. 

The  existence  of  time-varying  or  unsteady  forces  and  pressure 
results  in  noise,  vibration,  and  a  time-varying  propulsive  performance 
of  a  turbomachine.  In  general,  these  results  are  undesirable,  and  while 
it  is  not  possible  to  eliminate  their  existence,  it  is  desirable  to 
minimize  their  occurrence. 

In  a  turbomachine,  the  existence  of  upstream  blades  means  that  the 
blades  will  experience  a  time-varying  flow  and  hence  unsteady  forces  due 
to  the  relative  motion  between  blade  rows.  Adjacent  blade  rows  of  turbo¬ 
machines  are  usually  placed  as  closely  together  as  mechanical  constraints 
will  allow.  In  such  an  arrangement,  dictated  by  the  need  to  save  engine 
cost  and  weight,  some  blade  rows  may  be  subjected  to  strong  potential 
pressure  fields  from  the  downstream  rows.  A  rotor  passing  through  a 


potential  field  of  a  stator  will  be  subjected  to  time-dependent  pressure 
fluctuations  which  may  cause  significant  vibratory  stresses.  In  addition 
to  spatial  velocity  variations  from  wakes,  similar  variations  can  arise 
due  to  wall  boundary  layers,  wall  separation,  and  the  existence  of  up¬ 
stream  bends  and  turns. 

In  recent  years,  the  radiated  noise  and  vibration  have  become 
important  operating  characteristics  of  the  turbomachine  requiring  the 
consideration  of  the  designer.  A  portion  of  the  noise  and  vibration 
created  in  a  turbomachine  is  the  result  of  the  time-dependent  interaction 
of  its  blades  with  disturbances  in  the  flow,  both  spatial  and  temporal. 
This  interaction  results  in  an  unsteady  or  time-dependent  lift  on  the 
blades.  Thus,  the  design  of  turbomachines  with  noise  and  vibration  levels 
which  meet  the  established  standards  requires  that  time-dependent  design 
methods  be  available.  Unfortunately,  the  complexity  of  the  flow  through 
a  turbomachine  has  made  this  a  tremendous  problem  for  the  designer.  The 
flow  through  an  axial  flow  compressor,  fan,  or  turbine  stage  is  generally 
far  more  complicated  than  the  flow  over  an  isolated  wing  or  wing-body 
combination.  This  is  because  of  blade-to-blade  interference,  hub-to-tip 
variations  along  and  between  blades,  interactions  between  rotating  and 
non-rotating  stages  due  to  circulation,  thickness  and  wake  effects,  and 
inlet  distortions.  The  design  methods  resulting  from  the  assumption  of 
steady,  time- independent  flow  do  not  allow  for  the  consideration  of  the 
unsteady  lift.  Higher  relative  speeds,  higher  pressure  ratios,  and 
thinner  blades  have  compounded  the  research  and  design  problems  in  recent 
years. 

Distorted  inlet  flow  has  emerged  as  one  of  the  most  serious  problems 
facing  the  turbomachinery  designer.  Distorted  inlet  flow  is  a  term  which 


3 


is  generally  used  to  describe  the  flow  entering  a  turbomachine  where  the 
flow  had  radial  or  circumferential  variations  in  pressure,  velocity,  or 
temperature.  Although  both  radial  and  circumferential  distortions  are 
present  in  compressors  and  fans,  the  circumferential  distortion  poses  an 
additional  complication  since  the  flow  relative  to  a  moving  blade  row  will 
be  unsteady.  Inlet  flow  distortions  are  a  very  realistic  and  prevalent 
problem  in  modern  aircraft  jet  engines  since  the  engines  must  operate 
under  many  different  conditions.  They  can  be  caused  by  many  sources, 
such  as  shock-boundary  layer  interaction  in  a  supersonic  inlet,  crosswind- 
induced  inlet  separation,  and  armament  firing.  These  and  similar  condi¬ 
tions  cause  pressure,  temperature,  and  foreign  gas  type  of  distortions, 
and  the  result  is  generally  reduced  performance  and  stability  of  the 
compressor  system. 

In  the  flight  environment,  a  turbojet  engine  must  contend  with  dis¬ 
turbances  due  to  atmospheric  turbulence  and  the  aircraft  flow  field,  as 
well  as  with  internal  perturbations  caused  by  the  inlet  and  its  controls. 
Therefore,  the  axial  flow  compressor  may  be  required  to  operate  at  flow 
rates  and  pressure  ratios  other  than  the  flow  rate  and  pressure  ratio  for 
which  it  is  designed.  This  especially  occurs  in  the  operation  of  the 
high  pressure  ratio  compressor  in  the  aircraft  gas  turbine  cycle.  The 
operating  conditions  of  the  compressor  at  speeds  other  than  the  design 
speed  are  defined  by  the  performance  of  the  turbine  driving  the  compressor, 
the  area  of  the  outlet  nozzle,  and  the  aircraft  speed  and  altitude.  There 
is  an  unstable  limit  of  operation  known  as  surge,  marked  by  a  complete 
breakdown  of  the  continuous  steady  flow  through  the  whole  compressor, 
resulting  in  large  fluctuations  of  flow  with  time.  When  the  incidence 
angle  and  subsequently  the  loading  on  an  airfoil  is  increased,  a  point  is 


reached  when  the  flow  will  separate  from  the  suction  surface  of  the  air¬ 
foil.  When  this  occurs  in  a  compressor,  surge  can  occur.  Surge  will 
result  in  very  large  fluctuating  forces  on  the  blades  which  may  easily 
lead  to  oscillations  violent  enough  to  destroy  the  compressor  blades. 

Another  important  unstable  condition  connected  with  unsteady  inlet 
flow  which  can  lead  to  compressor  damage  is  rotating  stall.  This  is  a 
phneomenon  related  to  compressors  working  with  throttled  through  flow  and 
usually  at  rotational  speeds  lower  than  nominal  speed.  If  a  blade  row  is 
on  the  verge  of  stalling,  it  is  likely  that  one  particular  blade  will 
stall  before  another  neighboring  blade.  A  single  stalled  blade  results 
in  a  flow  diversion  which  tends  to  overload  one  adjacent  blade  and  unload 
the  other  adjacent  blade.  The  overloaded  blade  will  then  stall  and  the 
resulting  flow  diversion  will  unload  the  originally  stalled  blade.  This 
process  is  repeated,  and  the  stall  cell  (i.e.,  the  region  of  very  low  flow 
rate)  propagates  in  the  opposite  direction  as  the  rotation  of  the  rotor. 

In  actual  cases  there  may  be  more  than  one  stall  cell  rotating  on  a  single 
blade  row.  In  general,  the  zones  of  stall,  which  may  cover  several  blade 
spacings,  rotate  at  about  half  the  rotor  blade  speed.  The  alternate 
loading  and  unloading  of  the  blades  sets  up  an  alternating  stress  on  each 
blade.  This  stress  in  itself  is  not  large  unless  the  forcing  frequency 
happens  to  match  a  blade  vibrational  frequency.  In  this  case,  large 
stresses  are  present  and  fatigue  failure  can  occur,  resulting  in  complete 
destruction  of  an  entire  blade  row.  Little  can  be  done  to  prevent  this 
failure,  other  than  avoiding  operation  at  speeds  corresponding  to 
resonance.  Rotating  stall  is  more  dangerous  than  other  unsteady  phenomena 
to  the  survival  of  the  compressor. 

Both  surge  and  rotating  stall  are  unstable  conditions  which  may 
severely  damage  a  turbomachine.  In  order  to  prevent  this  from  occurring. 


the  stall  limit  as  a  function  of  pressure  rise  and  flow  rate  is  estab¬ 
lished  and  controls  are  designed  to  prevent  the  machine  from  entering 
the  unstable  region  of  its  operating  range.  An  engine  has  to  satisfy  a 
number  of  often  conflicting  requirements  such  as  adequate  stability,  low 
fuel  consumption,  long  life,  low  weight,  and  low  cost.  To  achieve  an 
optimum  performance  compromise,  the  engine  is  designed  to  have  just  enough 
stability  margin  to  function  adequately.  The  required  stability  margin, 
which  is  defined  as  the  margin  between  the  engine  operating  line  and  the 
limit  of  engine  stable  operation,  has  to  allow  for  a  number  of  effects 
including  the  magnitude  and  extent  of  the  nonuniformity  of  the  flow 
entering  the  compressor.  Therefore,  the  problem  caused  by  distorted 
inlet  flow  is  that  the  stability  limit  is  lowered  by  the  distortion,  and 
unstable  conditions  will  occur  at  a  higher  flow  rate  and  lower  pressure 
rise.  It  is  therefore  important  for  the  designer  to  accurately  predict 
the  effect  of  distorted  inlet  flow  on  compressor  performance  and  to 
design  distortion  tolerant  compressors. 

The  phenomena  that  lead  to  performance  degradation  in  compressors 
subjected  to  distorted  inlet  flows  are  essentially  unsteady.  In  the 
absence  of  an  understanding  of  the  unsteady  flow  dynamics,  compressor 
performance  prediction  techniques  are  steady  state  in  concept.  Further, 
predictions  involving  stage  analysis  are  restricted  to  using  steady  state 
characteristics.  Unsteady  flow  modelling  techniques  have  been  applied 
successfully  to  cascades,  but  usually  with  the  limitations  of  lnviscid 
flow  and  small  perturbation  assumptions,  thus  preventing  investigations 
in  highly  nonlinear  areas  of  response  such  as  stall.  There  is  therefore 
considerable  need  for  the  improvement  of  flow  distortion  models,  but 
without  experimental  data  to  quantify  unsteady  response  behavior  of  blade 
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rows,  stages,  and  complete  compressors,  improvements  that  may  take  account 
of  unsteady  phenomena  are  difficult  to  make. 

1. 2  Status  of  Existing  Knowledge 

As  mentioned  in  the  previous  section,  the  existence  of  unsteady  lift 
on  the  moving  or  stationary  blades  of  a  turbomachine  is  caused  by  vibra¬ 
tion  of  the  blades,  potential  flow  interactions  due  to  the  relative  motion 
of  adjacent  blade  rows,  and  inlet  flow  distortions.  There  are  different 
ways  to  experimentally  study  the  unsteady  response  of  a  turbomachine 
rotor.  One  approach  is  to  obtain  measurements  of  the  unsteady  normal 
force  and  moment  acting  on  a  rotor  blade  by  instrumenting  the  blade  with 
strain  gages.  A  second  method  used  to  determine  the  unsteady  response 
of  a  rotor  in  a  spatially  varying  inlet  flow  is  to  make  measurements  of 
the  circumferential  distribution  of  time-mean  total  pressure  at  the  inlet 
and  exit  of  the  rotor.  From  these  measurements,  the  unsteady  circulation 
on  the  blades  caused  by  their  interaction  with  the  inlet  distortion  can 
be  determined. 

Perhaps  the  best  known  method  for  the  design  analysis  of  the 
unsteady  lift  and  moment  in  a  turbomachine  is  that  by  Kemp  and  Sears 
[1,  2].  They  described  the  aerodynamic  interference  between  a  stator 
and  rotor  with  the  simplifications  of  two-dimensional,  incompressible, 
isolated  airfoil  theory.  This  method  is  an  extension  of  the  earlier 
original  work  by  Von  Karman  and  Sears  [3],  and  Sears  [4],  which 
describes  the  unsteady  lift  generated  on  an  isolated,  flat  plate  air¬ 
foil  subjected  to  velocity  fluctuations  normal  to  its  chord.  While 
including  the  steady  interaction  of  the  other  airfoils  in  the  cascade, 
the  unsteady  interaction  and  hence  the  effects  of  cascade  spacing, 
are  neglected  by  Kemp  and  Sears.  They  consider  the  problem  of  mutual 


interference  of  blades  in  cascades  due  to  relative  motion  of  successive 


blade  rows.  They  also  extend  this  analysis  to  include  the  effects  of 
viscous  wakes  on  the  unsteady  forces  on  the  blades  in  cascades.  The 
method  of  Kemp  and  Sears  provides  a  solution  which  expresses  the  unsteady 
lift  as  a  function  of  reduced  frequency,  i.e.,  the  ratio  of  airfoil  chord 
length  to  disturbance  wavelength.  Although  the  unsteady  contribution 
of  the  adjacent  blades  of  the  cascade  is  neglected,  it  is  this  type  of 
solution  which  provides  the  designer  with  the  required  design  information. 

A  further  improvement  to  the  theory  of  unsteady,  incompressible 
aerodynamics  of  an  isolated  airfoil  was  presented  by  Horlock  [5]. 

Ho  lock  analyzes  velocity  disturbances  both  normal  and  parallel  to  a 
flat  plate  airfoil.  The  response  of  the  isolated  airfoil  to  the  chord- 
wise  disturbance  can  be  expressed  in  terms  of  the  Horlock  function  which 
provides  a  design  solution  similar  to  that  of  Sears  for  a  normal  disturb¬ 
ance.  Naumann  and  Yeh  [6]  extend  this  analysis  to  include  the  effect  of 
blade  camber.  A  series  of  design  charts  is  presented  which  show  the 
effects  of  blade  camber,  stagger  angle,  and  reduced  frequency  on  the 
generation  of  unsteady  lift,  assuming  the  turbomachine  blade  can  be 
represented  by  a  single  airfoil. 

The  inclusion  of  the  unsteady  contribution  of  the  neighboring  blades 
of  the  blade  row  has  been  considered  by  Whitehead  in  Reference  [7]  and  by 
Henderson  and  Daneshyar  in  Reference  [8]  for  an  incompressible  flow.  The 
analysis  by  Whitehead  permits  the  unsteady  lift  and  moment  on  a  cascade 
of  rigid  flat  plate  airfoils  interacting  with  upstream  wakes  to  be  deter¬ 
mined,  and  his  results  are  presented  in  the  form  of  extensive  tables  of 
lift  and  moment  coefficients  for  vibrating  airfoils  in  cascade.  This  is 
a  solution  of  the  direct  problem,  that  is,  it  analyzes  the  response  of 
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a  known  cascade  of  airfoils.  The  method  of  Henderson  and  Daneshyar  gives 
a  functional  relationship  for  the  fluctuating  lift  on  a  cascade  of  air¬ 
foils  similar  to  the  analyses  for  an  isolated  airfoil  presented  in 
References  [3],  [A],  [5],  and  [6].  The  analysis  in  Reference  [8] 
includes  the  unsteady  contribution  of  the  entire  blade  row  and  thus, 
represents  an  extension  of  existing  indirect  or  design  solutions  of  the 
unsteady  response  of  an  isolated  airfoil  to  a  cascade  of  airfoils.  How¬ 
ever,  for  the  designer,  the  use  of  this  type  of  solution  is  very  laborious 
and  time  consuming,  since  a  series  of  calculations  must  be  conducted  to 
determine  a  configuration  with  minimised  unsteady  lift. 

All  the  theoretical  analyses  just  discussed  utilize  similar  mathe¬ 
matical  representations  to  develop  the  resulting  solutions.  These  repre¬ 
sentations  replace  the  airfoil  by  a  distribution  of  vorticity  on  the 
blade  and  in  its  wake.  The  strength  of  this  vorticity  and  the  resulting 
induced  velocities  satisfying  the  boundary  conditions  on  the  airfoil 
surface  are  then  determined.  This  specifies  the  unsteady  pressure  dis¬ 
tribution  and  hence  the  unsteady  lift.  For  all  cases,  the  flow  is  assumed 
to  be  inviscid,  incompressible,  and  two-dimensional,  with  the  airfoils 
being  infinitely  thin  and  the  disturbance  velocities  small  in  comparison 
with  the  mean  or  free-stream  velocity. 

Linear  theory  continues  to  play  an  important  role  in  predicting  per¬ 
formance,  stability,  and  flutter  of  turbomachines,  especially  in  the 
difficult  cases  of  compressors  with  supersonic  relative  velocities  and 
subsonic  axial  velocities.  Most  of  the  current  effort  in  unsteady  research 
utilizes  the  two-dimensional  cascade  approximation,  for  which  a  large 
background  of  theory  has  been  developed  over  the  years.  Significant 


differences  from  isolated  airfoils  have  been  found,  and  it  has  been 


demonstrated  that  blade-to-blade  interference  can  have  a  large  unfavor¬ 
able  effect  on  flutter  characteristics. 

McCroskey,  in  his  review  of  research  in  unsteady  fluid  dynamics 

[9] ,  cites  a  number  of  studies  which  are  representative  of  the  most 
recent  activity  in  the  field  of  unsteady  research.  Verdon  and  McCune 

[10] ,  Jones  and  Moore  [11],  and  Rao  and  Jones  [12]  have  developed 
somewhat  different  theories  in  terms  of  the  subsonic  velocity  potential 
for  flat  plate  blades.  They  describe  the  effects  of  varying  reduced 
frequency,  stagger  angle,  interblade  phase  angle,  interblade  spacing, 
and  Mach  number,  and  they  indicate  the  possibility  of  torsional  flutter 
over  a  wide  range  of  conditions.  Ni  and  Sisto  [13]  developed  a  fully 
numerical  time-marching  scheme  for  small  unsteady  perturbations  about 

a  steady  mean  flow.  Both  subsonic  and  supersonic  results  were  presented, 
and  they  found  indications  of  torsional  instability  in  the  supersonic 
case. 

Adamczyk  [14],  and  Goldstein  and  Atassi  [15]  also  considered  small- 
amplitude  fluctuations  about  a  steady  mean  flow,  along  with  finite  blade 
thickness,  camber,  and  mean  incidence.  Their  studies  are  detailed 
mathematical  analyses  of  the  response  to  incompressible  gusts,  and  their 
results  show  a  large  influence  of  blade  thickness  and  mean  incidence 
angle.  The  analysis  of  Horlock,  Greitzer,  and  Henderson  [16]  is  an 
appropriate  method  which  gives  analytical  expressions  for  the  unsteady 
lift  on  a  moving  cascade  of  closely  spaced  airfoils  of  arbitrary  finite 
chord.  It  is  a  much  simpler  analysis  since  it  takes  advantage  of  the 
channel-like  flow  behavior  between  closely  spaced  blades.  It  shows 
explicitly  the  variation  of  the  fluctuating  lift  with  the  various 


parameters,  and  it  clearly  shows  the  limiting  behavior  of  cascades 
exposed  to  low-frequency  gusts. 

Most  of  the  unsteady  cascade  theories  available  at  the  present  time 
are  either  directly  applicable  to  inlet  distortion  calculations  or  are 
applicable  after  proper  modification  of  the  boundary  conditions.  Two 
significant  recent  developments  appear  to  be  the  inclusion  of  three- 
dimensional  effects  using  small  disturbance  theory  (i.e.,  the  blades  are 
assumed  to  be  flat  plates  or  slightly  cambered  blades  set  at  zero  or  small 
incidence  and  the  unsteady  disturbances  are  assumed  to  be  small)  as  demon¬ 
strated  by  Namba  in  Reference  [17]  and  the  incorporation  of  viscous  loss 
hysteresis  and  effects  of  non-uniform  exit  pressure  into  a  multi-segment 
parallel  compressor  model  as  demonstrated  by  Mazzawy  in  Reference  [18]. 
Both  approaches  are  likely  to  contribute  toward  a  better  understanding  of 
the  major  physical  effects  which  are  presently  combined  in  the  empirically 
derived  distortion  indices  used  in  current  design  practice. 

In  contrast  to  the  large  number  of  theoretical  methods  available  for 
the  analysis  of  the  interaction  of  a  turbomachine  blade  row  with  upstream 
disturbances,  there  is  little  systematic  experimental  data  available.  As 
a  result,  the  validity  of  these  various  theoretical  methods  and  the 
assumptions  used  has  not  been  fully  determined.  There  is  a  significant 
need  to  demonstrate  whether  the  representation  of  the  unsteady  response 
of  a  blade  row  as  an  isolated  airfoil  is  justifiable  or  not. 

This  lack  of  experimental  measurements  of  the  unsteady  lift, 
particularly  in  a  rotating  blade  row,  is  due  to  the  complexity  and  size 
of  the  required  instrumentation.  The  development  of  highly  sensitive, 
miniaturized  pressure  transducers  has  resulted  in  the  initiation  of 
studies  in  several  laboratories  to  measure  the  unsteady  pressure 
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distributions  on  both  moving  and  stationary  blade  rows.  The  development 
of  a  simplified  experimental  method  for  the  determination  of  the  response 
of  a  blade  row,  particularly  a  moving  one,  is  desirable.  Such  a  tech¬ 
nique  is  suggested  by  the  analysis  by  Horlock  and  Daneshyar  [19],  and 
Hawthorne  [20]  concerning  the  variations  in  total  pressure  in  an  unsteady 
flow. 

These  analyses  reach  the  identical  conclusion  that  if  a  moving  blade 
row  has  unsteady  lift  on  its  blades,  there  will  exist  a  spatial  variation 
in  the  time-mean  total  pressure  downstream  of  the  blade  row.  This  spatial 
variation  can  be  related  to  the  unsteady  circulation  on  the  blades.  Thus, 
by  conducting  the  relatively  simple  measurement  of  the  time-mean  total 
pressure,  the  unsteady  response  can  be  indirectly  determined  from  a  know¬ 
ledge  of  the  unsteady  circulation.  This  technique  is  particularly  useful 
in  determining  the  unsteady  response  of  the  blade  row  to  variations  in 
blade  geometry  and  disturbance  characteristics.  It  also  suggests  a  method 
for  studying  experimentally  those  turbomachines  whose  size  and  physical 
construction  do  not  permit  the  installation  of  miniaturized  sensors  to 
directly  measure  the  unsteady  lift.  In  addition  to  the  time-mean  total 
pressure,  Hawthorne  also  analyzes  the  instantaneous  total  pressure  and 
shows  this  to  be  directly  related  to  the  unsteady  lift  on  the  blades. 

While  this  approach  would  give  a  direct  indication  of  the  unsteady  lift, 
it  requires  more  complex  instrumentation  to  determine  the  instantaneous 
total  pressure  than  is  required  to  determine  the  time-mean  total  pressure. 
Therefore,  the  most  straightforward  method  for  experimentally  determining 
the  unsteady  response  of  a  moving  turbomachine  blade  row  is  to  make 
measurements  of  the  time-mean  total  pressure.  This  method  was  used  by 
Henderson  and  is  thoroughly  presented  in  Reference  [21].  A  more  recent 
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investigation  which  helps  to  provide  the  kind  of  instantaneous  and  time- 
mean  response  data  that  are  needed  to  determine  the  range  of  validity  of 
various  linear  theories  is  that  of  Bruce  and  Henderson  in  Reference  [22], 
They  have  verified  some  of  the  cascade  theories  for  lift  and  moment  for 
actual  rotor  applications,  but  additional  losses  due  to  unsteady  effects 
were  found  that  are  beyond  the  scope  of  inviscid  linear  theory.  These 
losses  were  not  predicted  because  the  theories  assume  an  inviscid  flow. 

There  is  clearly  a  strong  and  urgent  need  in  both  rotors  and  cascades 
for  more  of  the  detailed  pressure  measurements  of  the  kind  just  mentioned 
as  well  as  of  the  kind  that  Fleeter  et  al.  in  Reference  [23]  and  Carta 
and  St.  Hilaire  in  Reference  [24]  have  begun  to  accumulate.  However,  just 
as  unsteady  cascades  are  more  complicated  than  isolated  airfoils  and 
depend  upon  many  more  parameters,  so  are  complete  turbomachines  still 
more  complex  than  cascades.  Therefore,  McCroskey,  in  Reference  [9], 
states  that  a  unified  plan  should  be  established  in  the  near  future  to 
avoid  a  surplus  of  unwieldy  experimental  and  numerical  data.  Additional 
simplified  analyses  in  the  manner  of  Horlock,  Greitzer,  and  Henderson 
[16]  will  continue  to  be  useful  in  identifying  dominant  flow  features, 
in  checking  elaborate  calculations,  and  in  understanding  the  extensions 
from  cascades  to  rotating  machines. 

Mathematical  models  have  been  developed  to  predict  the  response  of 
axial  compressors  to  distorted  inlet  flow.  The  various  types  of  compressor 
responses  to  inlet  distortions  are  loss  in  performance,  loss  in  stability 
(surge  or  stall),  attenuation  of  the  distortion,  perturbed  flow  fields, 
and  unsteady  blade  forces.  No  mathematical  model  exists  which  would 
predict  all  of  these  responses  for  every  compressor  and  every  distortion. 
Instead,  specialised  models  with  different  prediction  capabilities  are 
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available  for  specific  types  of  inlet  distortions  and  configurations  of 
compressors. 

The  parallel  compressor  model  is  used  to  predict  the  effects  of 
steady  circumferential  total  pressure  and/or  total  temperature  distortions 
on  the  performance  and  stability  of  compressors.  The  model  divides  the 
distorted  compressor  into  sections  which  are  treated  as  hypothetical 
compressors  working  in  parallel  in  a  quasi-steady  flow.  The  parallel 
compressor  model  is  only  suitable  for  circumferential  inlet  distortions 
because  radial  distortions  induce  a  great  amount  of  radial  flow  redistri¬ 
bution  which  contradicts  a  main  assumption  that  there  is  no  connection 
between  the  parallel  compressors. 

Roberts,  Plourde,  and  Smakula  [25]  measured  the  attenuation  of  a 
180  degree  square  wave  total  pressure  inlet  distortion  in  a  single  stage, 
two-stage,  and  three-stage  compressor  configuration.  By  utilizing  the 
parallel  compressor  model  and  comparing  the  experimental  and  analytical 
results,  they  obtained  reasonable  agreement  between  experiment  and 
theory.  In  a  further  experiment,  Roberts,  Plourde,  and  Smakula  measured 
the  performance  characteristics  of  a  three-stage  compressor  subjected 
first  to  a  mild  and  then  to  a  severe  total  pressure  inlet  distortion. 

Mild  and  severe  total  pressure  inlet  distortions  are  characterized  by 
small  and  large  ratios  of  the  distorted  total  pressure  amplitude  to  the 
undlstorted  or  uniform  total  pressure  amplitude.  They  compared  the 
results  with  predictions  obtained  with  the  parallel  compressor  model  and 
obtained  good  agreement  between  experiment  and  theory  for  the  mild  dis¬ 
tortion.  However,  the  agreement  was  found  to  be  poor  for  the  severe 
distortion  near  the  peak  pressure  coefficient  since  the  assumption  of 
parallel  compressor  operation  on  the  undistorted  flow  is  invalid  under 
the  operating  conditions  that  were  used. 
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It  is  not  immediately  apparent  from  the  parallel  compressor  analysis 
what  steps  should  be  taken  in  compressor  design  to  achieve  a  rapid  atten¬ 
uation  of  steady  circumferential  inlet  distortions  through  the  compressor. 
It  is  also  of  interest  to  determine  the  impact  of  the  assumption  that  all 
parallel  compressors  operate  on  the  undistorted  flow  (i.e.,  zero  circum¬ 
ferential  crossflow)  and  also  to  get  an  idea  of  the  flow  redistribution 
upstream  of  the  distorted  compressor.  The  answers  to  these  problems  can 
be  obtained  from  the  two-dimensional  linear  body-force  model  presented  by 
Plourde  and  Stenning  in  Reference  [26].  The  theory  is  restricted  to 
compressors  of  high  hub/ tip  ratio,  so  that  the  flow  field  can  be  regarded 
as  two-dimensional.  Further,  the  assumption  is  made  that  the  flow  is 
inviscid  and  incompressible.  The  compressor  is  assumed  to  consist  of  a 
large  number  of  small  stages,  and  the  effect  of  these  stages  on  the 
increase  in  pressure  is  simulated  by  an  appropriate  distribution  of  axial 
body  forces. 

Compressor  stall  at  high  rotor  speeds  is  generally  initiated  in  the 
later  stages  of  a  multi-stage  machine.  To  deal  with  a  steady  circum¬ 
ferential  total  pressure  inlet  distortion,  it  is  therefore  important  to 
design  the  front  stages  for  rapid  distortion  attenuation.  The  body-force 
analysis  is  not  very  suitable  for  providing  guidelines  for  such  a  design, 
since  it  assumes  all  stages  to  have  an  equal  performance.  Better  guide¬ 
lines  may  be  obtained  from  a  model  which  replaces  the  blade  rows  of  the 
compressor  individually  by  actuator  discs,  as  presented  in  References 
[27],  [28],  and  [29].  Actuator  disc  models  are  used  when  the  parallel 
compressor  model  or  the  body-force  model  is  invalid  and  also  when  the 
distortion  shape  and  individual  blade  row  attenuation  is  of  primary 
concern.  An  actuator  disc  is  a  means  of  representing  a  blade  row  by 
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shrinking  the  chord  length  of  the  blade  to  zero  and  concentrating  the 
effects  of  the  blade  rows  at  a  single  location  while  assuming  a  quasi¬ 
steady  flow.  The  model  can  be  described  for  compressors  with  high  hub/tip 
ratio,  thus  allowing  the  flow  field  to  be  regarded  as  two-dimensional. 

It  is  also  assumed  that  the  flow  outside  the  actuator  discs  is  inviscid 
and  incompressible.  The  assumption  that  the  velocity  perturbations  are 
small  compared  with  the  mean  axial  velocity  permits  the  equations  of 
motion  to  be  linearized.  The  axial  spacing  of  the  discs  corresponds  to 
the  gaps  between  the  blade  rows  since  circumferential  flow  redistribution 
can  only  occur  in  these  gaps.  The  flow  fields  upstream  and  downstream  of 
the  blade  row  can  then  be  solved  either  by  an  analytical  solution  of 
linearized  equations  or  by  numerical  methods.  The  solution  of  the 
linearized  inviscid  flow  equations  is  suitable  for  studying  distortion 
attenuation  if  the  distortion  is  small.  The  nonlinear  equations  and  the 
unsteady  response  must  be  used  for  large  distortions  and  when  stall 
phenomena  is  investigated.  The  parallel  compressor  and  actuator  disc 
theories,  which  were  just  discussed,  do  not  accurately  predict  the  unsteady 
response  of  the  blade  row  since  they  are  quasi-steady  theories.  A  detailed 
description  and  summary  of  these  types  of  overall  performance  and  engine 
stability  prediction  techniques  is  presented  by  Mokelke  in  Reference  [30]. 

The  theoretical  modelling  of  unsteady  flow  effects  in  turbomachines 
is  largely  limited  to  the  analysis  of  lightly  loaded  cascades  (i.e.,  to 
two-dimensional  flow  assumptions).  The  methods  used  involve  direct 
extensions  of  the  external  flow  solutions  for  unsteady  flow  past  single 
airfoils  in  unbounded  flow  to  the  more  difficult  boundary  value  problem 
of  unsteady  flow  past  airfoil  arrays.  Thus,  there  is  a  direct  relation¬ 
ship  with  the  current  status  of  unsteady  external  aerodynamics,  and  the 
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methods  used  in  this  field,  such  as  airfoil  modelling  by  vortex  distri¬ 
butions,  conformal  mapping,  transform  techniques,  etc.,  are  directly 
applicable  to  turbomachinery  unsteady  aerodynamics.  Furthermore,  it 
should  be  noted  that  the  methods  developed  for  the  flutter  and  forced 
vibration  analysis  of  turbomachinery  blading  are  equally  applicable  to 
turbomachinery  noise  studies. 

Within  the  status  of  knowledge,  the  turbomachinery  designer  is 
restricted  to  the  use  of  isolated  airfoil  theories  and  two-dimensional 
analyses  of  lightly  loaded  cascades  to  predict  the  unsteady  response  of 
a  blade  row.  Thus,  the  range  of  applicability  and  the  usefulness  of  the 
current  methods  are  severely  limited.  Also,  the  practical  utility  of 
these  current  methods  is  further  hampered  by  the  lack  of  experimental 
measurements  and  the  inaccessibility  of  the  currently  available  results. 

As  a  result,  the  appreciation  for  the  unsteady  aerodynamic  effects  of 
turbomachines  and  the  access  to  the  available  results  is  limited  to  a 
small  number  of  specialists.  However,  a  simplified  experimental  technique 
which  employs  measurements  of  the  time-mean  total  pressure  suggests  a  more 
expedient  manner  for  obtaining  the  necessary  data. 

Platzer,  in  his  review  on  current  developments  in  unsteady  turbo- 
machinery  aerodynamics  [31],  notes  that  a  significant  change  in  attitude 
and  increase  in  research  effort  has  taken  place  in  the  past  decade. 

He  concludes  by  distinguishing  the  following  major  developments  and 
results: 

1.  Incompressible,  subsonic  and  supersonic  unsteady  flow  past 
lightly  loaded  single  cascades  has  been  studied  in  consid¬ 
erable  detail,  but  relatively  little  comparative  evaluation 
and  documentation  of  the  various  solutions  is  available  in 
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2.  Recently,  solutions  for  incompressible  unsteady  flow 
past  highly  loaded  cascades  have  also  been  published, 
but  little  information  is  yet  available  for  subsonic 
and  supersonic  cascades  with  finite  blade  thickness 
and  loading. 

3.  The  unsteady  transonic  cascade  problem  is  still 
unsolved,  but  some  new  approaches  are  presently  in 
progress  which  are  likely  to  provide  new  insight. 

4.  Work  on  the  analysis  of  three-dimensional  unsteady 
flow  effects  is  in  progress  and  some  valuable  new 
results  for  lightly  loaded  rotors  have  already  been 
published. 

5.  New  measuring  techniques  coupled  with  modern  data 
acquisition  systems  are  being  applied  to  the 
investigation  of  unsteady  cascade  and  rotor  flows 
throughout  the  entire  speed  regime.  A  number  of 
papers  have  recently  been  published  which  provide 
detailed  unsteady  flow  and  pressure  information  for 
comparison  with  the  available  theoretical  models. 

6.  Progress  has  also  been  made  towards  the  prediction 
of  rotating  stall  and  surge  using  nonlinear  unsteady 
analyses,  but  much  more  experimental  and  theoretical 
work  is  required  to  understand  and  describe  highly 
viscous  and  separated  unsteady  flow  phenomena  such 
as  wake  formations  and  interactions,  stall  flutter, 
etc. 
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In  order  to  promote  acceptance  and  use  of  unsteady  flow  information 
in  turbomachinery  research  and  design,  increasing  emphasis  will  have  to 


be  placed  on  the  problem  of  data  access,  evaluation,  and  documentation 
[31].  Therefore,  systematic  comparisons  and  evaluations  of  various 
theoretical  methods  and  experimental  results  should  be  conducted  for 
well-defined  test  configurations.  Further  improvements  in  aircraft 
engine  performance  and  the  avoidance  of  aerodynamic  and  aeroelastic 
stability  problems  which  cause  expensive  engine  modifications  will  depend 
on  the  continued  research  effort  in  unsteady  turbomachinery  aerodynamics 
and  on  the  transfer  of  highly  specialized  unsteady  aerodynamic  and  aero¬ 
elastic  information  to  the  design  community. 

1. 3  Statement  of  the  Problem 

In  dealing  with  the  problem  of  distorted  inlet  flows,  designers  of 
axial  flow  fans  and  compressors  require  a  knowledge  of  the  fluctuating 
forces  on  the  blades,  the  attenuation  of  the  distortion  through  each 
blade  row,  the  loss  in  performance,  and  any  reduction  in  stability  caused 
by  the  distortion,  as  a  function  of  the  operating  and  geometrical  design 
variables  of  the  machine  which  the  designer  can  control.  Theoretical 
predictions  are  valuable  in  understanding  distortion  phenomena  and  also 
in  making  some  design  decisions;  however,  theory  has  several  shortcomings 
when  used  in  the  design  process.  For  example,  when  trying  to  predict 
distortion  attenuation,  performance  loss,  or  stability  limits  of  a  com¬ 
pressor  or  fan,  theoretical  methods  such  as  the  parallel  compressor  theory 
and  actuator  disc  model  require  the  steady  performance  of  the  compressor 
or  blade  row  to  be  known.  Since  the  performance  characteristics  are 
determined  by  a  number  of  geometric  design  parameters,  these  theoretical 
methods  give  very  little  insight  into  the  selection  of  the  proper  design 
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parameters  to  be  used  in  the  design  process.  Unsteady  airfoil  theories 
and  two-dimensional  cascade  analyses  can  be  used  for  selecting  design 
parameters  to  minimize  the  fluctuating  lift  on  the  blade,  but  these  methods 
generally  apply  to  lightly  loaded  blades  and  neglect  viscous  effects. 

To  aid  in  the  continual  development  of  theoretical  design  techniques 
and  to  help  bridge  the  gap  between  present  day  theories  and  the  needs  of 
designers,  experimental  data  of  a  fundamental  nature  demonstrating  the 
effect  of  geometric  and  operating  variables  are  needed.  Most  of  the  data 
available  to  date  are  not  of  this  nature.  Extensive  studies  have  been 
made  on  modern  aircraft  engines  to  determine  their  response  to  distortions 
and  to  identify  problem  areas.  These  studies,  however,  are  usually  part 
of  a  continuing  development  program  aimed  at  determining  safe  operating 
limits  and  suggesting  modifications  rather  than  providing  the  fundamental 
data  needed  for  designing  the  next  generation  of  engines. 

The  purpose  of  this  investigation  was  to  study  the  unsteady  response 
of  an  axial  flow  turbomachine  blade  row  to  time-mean  spatial  variations 
in  the  rotor  inlet  flow.  This  program  is  a  combination  of  experimental 
and  analytical  efforts  which  are  directed  toward  obtaining  information 
relating  the  influence  of  rotor  geometrical  and  operating  parameters, 
such  as  solidity,  stagger  angle,  and  steady  load,  on  the  unsteady  rotor 
response.  Specifically,  the  effect  of  the  rotor  employed  in  various 
geometric  configurations  on  the  distortion  flow  field  was  studied.  Two 
test  rotor  configurations  were  employed  in  the  experimental  portion  of 
this  investigation.  The  effects  of  rotor  solidity,  stagger  angle,  and 
steady  load  (i.e.,  mean  incidence  angle)  were  examined  using  a  six-bladed 
uncambered  rotor.  The  effects  of  rotor-stator  spacing  and  radial  varia¬ 
tions  of  the  flow  field  were  investigated  with  a  cambered  test  rotor 
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consisting  of  nine  blades.  The  effects  of  radial  variations  were  deter¬ 
mined  by  measuring  the  flow  field  at  various  spanwise  locations  along  the 
blade.  The  approach  taken  to  measure  the  response  of  the  rotor  was  to 
make  time-mean  flow  measurements  upstream  and  downstream  of  the  rotor. 
Detailed  circumferential  surveys  of  the  time-mean  flow  field  upstream  and 
downstream  of  the  rotor  were  made. 

In  addition  to  the  experimental  program,  a  theoretical  analysis  is 
presented  which  predicts  the  time-mean  total  pressure  at  the  exit  of  the 
rotor.  This  theory  is  an  unsteady  cascade  analysis,  and  it  includes  the 
unsteady  blade-to-blade  interaction  effects.  Thus,  a  basis  of  comparison 
is  provided  between  experimental  measurements  and  theoretical  predictions 
using  a  two-dimensional,  unsteady,  inviscid,  cascade  analysis. 

The  objective  of  this  study,  then,  is  to  provide  some  of  the 
necessary  fundamental  data  by  examining  the  response  of  a  turbomachinery 
rotor  as  a  function  of  blade  geometry  and  type  of  distortion.  Comparison 
of  these  results  with  the  predicted  values  of  time-mean  total  pressure 
downstream  of  the  rotor  is  also  used  to  determine  the  validity  of  the 
unsteady  analysis  and  the  practicability  of  the  experimental  technique. 
These  data  will  form  a  basis  from  which  design  data  can  be  derived  which 
should  result  in  axial  flow  fan  configurations  with  improved  performance 
when  operated  in  an  inlet  distortion. 


CHAPTER  II 


THEORETICAL  ANALYSIS 

2. 1  Prediction  of  Time-Mean  Total  Pressure  at  the  Exit  of  a 

Rotating  Blade  Row 

The  operation  of  a  turbomachine  in  a  spatially  varying  inflow  will 
result  in  the  generation  of  time-varying  or  unsteady  pressures  and  forces 
on  the  blades.  This  unsteady  response  will  also  lead  to  a  variation  in 
the  time-mean  flow  as  it  passes  through  the  blade  row.  From  a  knowledge 
of  the  absolute  flow  in  front  of  and  behind  the  blade  row,  it  is  possible 
to  determine  the  unsteady  response  of  the  blades.  Specifically,  it  is 
possible  to  determine  the  unsteady  circulation  on  the  blades  which  can 
then  be  used  to  predict  the  unsteady  effect  of  the  turbomachine  rotor  on 
the  inflow  distortion.  The  characteristics  of  the  distortion  after  its 
passage  through  the  rotor  are  important  for  the  determination  of  the 
unsteady  interaction  of  subsequent  blade  rows. 

The  analyses  of  Horlock  and  Daneshyar  in  Reference  [19]  and  Hawthorne 
in  Reference  [20]  show  that  if  a  moving  blade  row  has  an  unsteady  lift  on 
its  blades,  then  a  spatial  variation  in  time-mean  total  pressure  which  is 
associated  with  the  unsteady  lift  will  exist  downstream  of  the  blade  row. 
This  variation  in  time-mean  total  pressure  can  be  related  to  the  unsteady 
circulation  on  the  blades,  thus  describing  its  unsteady  response.  The 
work  of  Horlock  and  Daneshyar  also  suggests  an  experimental  procedure 
which  permits  the  unsteady  response  of  a  rotating  blade  row  to  spatial 
variations  in  its  inlet  flow  to  be  determined  from  measurements  of  the 
time-mean  total  pressure.  Henderson  [32]  employs  this  experimental 
procedure  to  determine  the  unsteady  circulation  on  a  rotating  blade  row 
exposed  to  a  spatially  varying  inlet  flow. 
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By  reversing  the  procedure  described  by  Henderson  in  Reference  [32], 
it  is  possible  to  predict  the  circumferential  variation  of  time-mean 
total  pressure  at  the  exit  of  the  blade  row.  This  is  precisely  the 
problem  which  was  undertaken  in  this  investigation.  In  the  following 
sections,  the  theoretical  development  is  presented  which  results  in  the 
equation  used  to  predict  the  time-mean  total  pressure  at  the  exit  of  the 
rotor.  Horlock  and  Daneshyar  present  the  derivation  of  the  equation 
describing  the  variation  of  time-mean  total  pressure  in  unsteady  flow. 

They  derive  this  result  using  two  different  -"alyses.  For  a  complete 
derivation  of  the  time-mean  total  pressure  downstream  of  a  row  of  blades 
with  fluctuating  circulation  using  a  vortex  analysis,  refer  to  Reference 
[19].  The  changes  in  total  pressure  produced  by  movement  of  rows,  arrays, 
and  sheets  of  vortices  are  derived.  This  result  can  also  be  derived  from 
energy  considerations,  and  a  detailed  formulation  of  this  procedure  is 
presented  here. 

2.1.1  Steady  Transport  of  a  Distortion  Through  a  Blade  Row. 

Consider  a  two-dimensional  cascade  of  airfoils  operated  in  a  distorted 
flow  which  varies  sinusoidally  along  the  face  of  the  cascade  as  shown  in 
Figure  1.  The  variation  of  total  pressure  shown  in  Figure  1  exists  in  the 
presence  or  absence  of  the  rotor.  The  presence  of  the  rotor  will  induce 
a  circumferential  velocity  upstream  of  the  rotor  which  will  influence  the 
axial  velocity  and  static  pressure.  As  the  distortion  is  carried  through 
the  moving  cascade,  it  is  transported  in  the  circumferential  direction  by 
an  amount  dependent  upon  the  amount  of  mean  swirl  (i.e.,  circumferential 
velocity)  introduced  by  the  blades. 

The  velocity  diagram  of  the  flow  relative  to  the  rotor  is  shown  in 
Figure  2.  Also  shown  are  the  relative  mean  streamlines  of  the  flow  as  it 
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Figure  2.  Flow  Relative  to  Rotor 


passes  through  the  rotor  for  the  cases  of  zero  steady  lift  and  non-steady 
lift  on  the  blades.  Assume  that  the  time-mean  total  pressure  distri¬ 
bution  at  the  cascade  inlet  can  be  expressed  as: 
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From  Figure  2,  it  is  seen  that  the  time  for  a  particle  of  fluid  to  be 
transported  through  the  cascade  is: 


At 


(2) 


where  W  is  the  mean  relative  velocity  through  the  cascade.  In  this 
m  - - - 

length  of  time,  a  fluid  particle  is  transported  in  the  circumferential 
direction  a  distance  equal  to  AWy*At  where  AW^  is  the  change  in  time- 
independent  tangential  or  circumferential-mean  velocity  relative  to  the 
cascade.  Thus,  at  the  exit  of  the  cascade  the  inlet  distortion  given  in 
Equation  (1)  has  the  form: 
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where  y  is  the  circumferential  measurement  location.  Since 

AW  *  W  -  W  »  V  (taniS.  -  tanB.)  if  V  -  constant.  Equation  (3) 
y  y^  y2  x  l  l  a 

becomes: 
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The  steady  lift  coefficient  on  the  blades  can  be  written  as  (see 
Reference  [33]): 

C  *  2 (~)  [tang,  -  tang,]  cosg  -  C  tang  (7) 
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Therefore,  substitution  of  Equation  (7)  into  Equation  (6)  gives: 
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This  is  the  expression  for  the  inlet  distortion  downstream  of  the  rotor 
blade  row  after  it  has  been  transported  through  the  blades  with  the 
circumferential-mean  or  steady  flow  assuming  that  its  amplitude  A^  is 
constant,  i.e.,  there  is  no  mixing  or  diffusion  of  the  distortion. 

2.1.2  Influence  of  Unsteady  Blade  Forces.  The  forces  on  the 
blades  of  an  isolated  turbomachine  rotor  relative  to  the  blade  row  are 
constant  if  the  flow  entering  the  blade  row  is  uniform.  This  implies 
that  the  total  pressure  must  increase  across  the  blade  row  if  work  is 
done  on  the  fluid.  A  total  pressure  probe  located  downstream  of  the 
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rotor  In  the  absolute  reference  frame  (l.e.,  fixed  with  respect  to  the 
casing  wall)  will  sense  an  instantaneous  total  pressure  variation  which 
will  be  constant  when  it  is  averaged  with  time.  Therefore,  the  time- 
mean  total  pressure  downstream  of  a  rotating  blade  row  with  time  invariant 
forces  is  constant. 

If  the  blades  experience  an  unsteady  force,  Horlock  and  Daneshyar 
[19]  and  Hawthorne  [20]  have  shown  that  the  time-mean  total  pressure  will 
vary  in  the  absolute  frame  of  reference.  This  variation  is  a  function  of 
the  time-mean  absolute  flow  angle  at  the  rotor  exit  and  the  unsteady 
circulation  on  the  blades.  Thus,  it  is  possible  to  determine  the  unsteady 
characteristics  of  the  blades,  in  particular  the  unsteady  circulation, 
from  measurements  of  the  time-mean  total  pressure.  This  procedure  was 
demonstrated  in  Reference  [32]. 

While  the  relationship  of  the  time-mean  total  pressure  downstream  of 

the  blades,  P_  ,  can  be  derived  by  a  vortex  analysis  of  the  blades,  as 
2 

was  done  in  References  [19]  and  [20],  it  can  also  be  derived  by  consider¬ 
ing  the  energy  exchange  across  the  blade  row,  as  shown  in  Reference  [19]. 
In  a  turbomachine,  the  flow  can  be  considered  as  adiabatic.  Therefore, 
the  energy  equation,  neglecting  potential  energy  effects  and  considering 
a  compressor  application,  can  be  written  as: 


w  . 


(9) 


This  expression  relates  the  energy  per  unit  mass  through  a  control  volume. 
If  a  control  volume  fixed  in  the  absolute  reference  frame  is  considered. 


as  shown  in  Figure  3,  then  the  mass  flow  rate  is  m  -  pVdy.  Therefore 
Equation  (9)  for  a  length  of  time  dt  becomes: 
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pVh  dydt  -  pVh  dydt  “  dW„, 
°2  °1 
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(10) 


where  dWcv  is  the  work  transferred  to  the  control  volume  in  time  dt.  As 
seen  in  Figure  3,  the  axial  velocity  V  is  assumed  to  be  constant  throughout 
the  control  volume. 

The  blades  shown  in  Figure  3  are  moving  in  the  circumferential  or  y 
direction  with  a  velocity  U  and  have  a  blade-to-blade  spacing  of  S. 

Hence,  for  an  incompressible,  isentropic  flow,  integration  of  Equation 
(10)  over  the  period  T  *  S/U  gives: 
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where  AWCV  is  the  total  work  transfer  to  the  control  volume  in  the  period 
T.  This  work  transfer  occurs  when  a  blade  passes  through  the  control 
volume. 

The  circulation  on  the  blade  is  unsteady  and  assumed  to  have  the 
harmonic  form 


f  -  rss  +  I? |e'i(vt  +  ""W  ,  (12) 

where  is  the  steady-state  circulation,  ( T (  is  the  magnitude  of  the 

unsteady  circulation,  V  is  the  frequency  of  the  variation  of  the  unsteady 

circulation  relative  to  the  moving  blades,  and  4>  is  the  phase  angle 

between  the  steady  and  unsteady  circulation.  <f>  is  referenced  to  the 

J.E 

blade  trailing  edge  since  that  is  where  the  vorticity  (or  circulation) 
is  shed  into  the  absolute  flow  or  wake  and,  hence,  enters  the  control 
volume.  The  addition  of  in  the  exponential  term  of  Equation  (12)  is 
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due  to  the  fact  that  the  quasi-steady  circulation  is  180  degrees,  or  tt 
radians,  out  of  phase  with  the  distortion.  All  terms  in  the  unsteady 
analysis  must  be  referenced  to  the  distortion,  thus  the  factor  of  it  must 
be  included.  From  Figure  3,  the  frequency  V  is  equal  to  -2itU/JI  and  the 
length  of  one  period  I  is  t  =  -(y  +  xtana^)/!).  And  since  e  *  costt  - 
isimr  =  -1,  Equation  (12)  can  be  expressed  as 

2ir 

f  =  rss  -|r|e"ilX  +  xtana2)  "  ^TE1  .  (13) 


It  is  assumed  that  the  lift  on  the  blades  and  their  circulation  are 

in  phase  with  each  other.  This  is  a  valid  assumption  for  low  values  of 

reduced  frequency.  Then  the  work  done,  AW  ,  is  the  product  of  the  y 

cv 

component  of  the  lift,  pVT",  and  the  distance  moved,  dy;  i.e.: 

AW^  -  pVfdy  .  (14) 


Recalling  that  the  period  T  -  S/U  and  combining  Equations  (11),  (13),  and 
(14)  gives: 
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or 
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But, 
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Substitution  of  Equation  (16)  into  Equation  (17)  gives: 
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And  finally. 
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Equation  (19)  is  the  final  expression  for  the  time-mean  total 
pressure  at  the  exit  of  the  rotor  whose  blades  experience  unsteady  forces 
and  circulation.  It  is  the  same  result  as  that  obtained  by  Horlock  and 
Daneshyar  in  Reference  [19]  from  a  vortex  analysis  if,  assuming  as  they 
did,  4>„  =  0.  Thus,  it  is  seen  from  Equation  (19)  that  as  the  flow  is 

li!# 

transported  through  the  moving  blades,  three  time-mean  total  pressure 

fields  are  superimposed  to  give  the  time-mean  total  pressure,  PT  ,  which 

2 

is  recorded  by  a  total  pressure  probe  located  downstream  of  the  blades. 
The  first  term  of  Equation  (19)  represents  the  upstream  disturbance  field 
which  is  transported  through  the  blades  with  the  circumferential-mean  or 
steady  flow,  and  it  is  given  by  Equation  (8);  the  second  term  represents 
the  steady  total  pressure  rise  which  arises  due  to  the  steady  lift  on  the 
blades  and  is  uniform  in  the  direction  of  blade  motion;  the  third  term 
represents  the  time-mean  total  pressure  due  to  the  unsteady  circulation 


on  the  blades. 
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2.2  Solution  of  the  Problem 

In  the  previous  section.  Equation  (19)  was  derived  which  predicts 
theoretically  the  time-mean  total  pressure  downstream  of  a  rotating  blade 
row.  The  equation  was  found  to  be  a  function  of  the  type  of  disturbance 
entering  the  rotor  and  the  geometric  and  operating  variables  of  the  rotor. 
The  purpose  of  this  section  is  to  give  an  explanation  of  how  the 
theoretical  problem  was  solved  and  how  the  various  operating  parameters 
used  in  the  solution  were  determined.  Also,  the  contribution  of  each 
separate  term  in  Equation  (19)  to  the  total  solution  will  be  shown. 

First  consider  the  time-mean  total  pressure  distribution  at  the 
cascade  inlet  given  by  Equation  (1).  This  total  pressure  distribution 
can  be  represented  by  a  standard  Fourier  series  as  a  function  of  0,  the 
circumferential  measurement  position,  as  follows: 

_____  a  00 

P-r  (9)  *  ~7T  +  2  a  cosn0  +  b  sinnG  .  (20) 
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Dividing  and  multiplying  by  the  square  root  of  the  sum  of  the  squares  of 
the  Fourier  coefficients  gives: 
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The  following  definitions  are  made  for  the  Fourier  coefficients: 

a  =  -A  sin4> 
n  n  n 

b  =  A  cosdi  , 
n  n  n 


and 
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2  2 

where  A  » »  a  +  b  and  6  is  the  phase  angle  defined  as: 
n  v  n  n  n 


a 

-1(— ) 

<J)  i  -tan 
n  n 


Using  these  definitions  for  a^  and  b^,  Equation  (21)  becomes: 
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And,  finally,  by  applying  the  trigonometric  relation  for  the  difference 
of  two  angles,  we  arrive  at  the  following  expression  which  describes  the 
sinusoidal  distribution  of  total  pressure  at  the  rotor  inlet: 


_  a  00 

PT  (6)  -  -f  +  l  A  sin(n0-<J>  )  .  (23) 
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Physically,  this  term  is  obtained  from  measurements  of  the  total  pressure 
inlet  distortion  without  the  rotor  operating. 

As  the  total  pressure  distortion,  given  by  Equation  (23),  passes 
through  the  rotor  blade  row,  it  undergoes  a  further  shift  in  phase  angle 
due  to  the  steady  lift  on  the  blades.  The  amount  of  phase  shift  due  to 
the  steady  transport  of  a  distortion  through  a  blade  row  was  discussed  in 
detail  in  Section  2.1.1.  Therefore,  assuming  the  amplitude  of  the 
distortion  is  not  altered,  i.e.,  A^  -  constant,  the  expression  for  the 
inlet  distortion  downstream  of  the  rotor  blade  row  after  it  has  been 
transported  through  the  blades  with  the  steady  flow  is: 
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The  calculation  of  lift  coefficient,  CT ,  was  made  from  the  following 

Ij 

equation: 


C  *  2(S/C) [tan0  -  tang.]  cosB  -  C_tanB 
l  X  L  m  D  m 
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The  angles  8.  and  8„  are  found  from  the  geometry  in  Figure  2.  8  is  then 
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given  by 
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The  drag  coefficient,  C^,  is  given  by: 
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The  steady  lift  on  the  blades  does  not  change  the  magnitude  of  the 

total  pressure  variation  downstream  of  the  rotor;  it  only  transports  the 

distortion  in  the  circumferential  direction.  In  a  real  physical  situation, 

viscous  effects  would  cause  a  change  in  magnitude.  However,  the  theory 

does  not  account  for  viscous  effects,  so  the  magnitude  of  the  amplitude 

of  the  total  pressure  is  unchanged  by  the  steady  lift. 

The  second  term  in  Equation  (19)  represents  the  steady  state  total 

pressure  rise,  P_  ,  -  P_  „  ,  and  it  is  also  a  result  of 

T  downstream  T  upstream 
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the  steady  lift  on  the  blades.  This  term  is  merely  the  differences  in 
total  pressure  between  the  upstream  and  downstream  measuring  locations  in 
a  uniform  flow.  It  does  not  shift  the  phase  angle  of  the  total  pressure 
disturbance;  it  only  adds  to  the  overall  magnitude  of  the  circumferential 
mean  value  inlet  distortion  as  it  passes  through  the  blade  row. 

The  last  term  which  must  be  considered  in  the  solution  of  the  time- 
mean  total  pressure  at  the  rotor  exit  is  the  unsteady  term.  The  unsteadi¬ 
ness  of  the  total  pressure  distortion  is  represented  by  the  magnitude  of 
unsteady  circulation  and  a  phase  angle.  This  phase  angle,  <j>  ,  is 

referenced  to  the  trailing  edge  of  the  blade.  Henderson,  in  Reference 
[21],  uses  measured  values  of  time-mean  total  pressure  to  predict  the 
unsteady  circulation  and  phase  angle  as  a  function  of  blade  stagger  angle, 
space-to-chord  ratio,  and  reduced  frequency.  These  predicted  values  of 
unsteady  circulation  and  phase  angle  are  used  in  the  solution  of  the 
problem  which  is  presently  being  discussed.  It  should  be  noted  that  the 
values  of  magnitude  and  phase  angle  of  the  circulation  which  are  calcu¬ 
lated  from  the  analysis  in  Reference  [21]  are  referenced  to  the  mid-chord 
of  the  blade.  The  following  transformation  is  necessary  to  reference  the 
circulation  and  phase  angle  to  the  blade  trailing  edge: 


'TE 


=  r 


iu> 


MC 


(28) 


[T  +  ir  ]  =  [f  -4-  ir  ]  (cosw  +  isinoj)  (29) 

real  imaginary  TE  real  imaginary  MC 


and 


4>te  =  _tan 


(r  ) 

imaginary  TE 
^real^TE 
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where  uj  is  the  reduced  frequency.  This  phase  angle,  4>  ,  is  measured 

TE 
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relative  to  the  steady,  or  quasi-steady  (u>=u),  circulation.  To  express 

this  phase  angle  relative  to  the  inlet  distortion,  P  ,  it  must  be 

11 

increased  by  ir  radians,  since  the  steady  circulation  is  180  degrees  out 
of  phase  with  the  distortion  field  [21].  Hence: 

[41  =  7T—  [  cb  1 

TE  distortion  TE  steady  circulation  . 

An  additional  phase  shift  arises  in  the  unsteady  term  which  represents  the 
amount  that  the  distortion  is  convected  after  its  passage  through  the 
rotor.  This  distance,  xtanc^,  is  shown  in  Figure  3.  Therefore,  the 
contribution  to  the  time-mean  total  pressure  due  to  the  unsteady  circula¬ 
tion  on  the  blades  is  given  by: 

PT  (6)  =  I  ££iLL  sin(Tr+ne-<)>TE  +  ■—  xtana2)  .  (31) 

2  n=l 

The  final  solution  of  the  time-mean  total  pressure  downstream  of  the 
rotor  blade  row  is  then  obtained  by  considering  all  the  terms  given  in 
Equations  (24)  and  (31)  plus  the  steady-state  total  pressure  rise. 

Equation  (19)  was  derived  in  the  previous  section  for  a  general  case  in 
exponential  form;  however,  only  the  sine  portion  of  the  general  form  needs 
to  be  considered  for  the  solution. 


CHAPTER  III 


EXPERIMENTAL  EQUIPMENT,  INSTRUMENTATION, 

AND  EXPERIMENTAL  STUDIES  CONDUCTED 

3.1  The  Axial  Flow  Research  Fan 

The  entire  experimental  program  was  conducted  in  a  facility  referred 
to  as  the  Axial  Flow  Research  Fan  (AFRF),  which  was  designed  and  con¬ 
structed  at  the  Applied  Research  Laboratory,  Hie  Pennsylvania  State 
University.  This  facility  was  designed  specifically  for  studying  the 
response  of  axial  flow  turbomachinery  blade  elements  and/or  a  single 
stage  to  subsonic  incompressible  unsteady  flows  of  the  type  encountered 
in  practice.  Thus,  the  overall  dimensions  of  the  test  rotor  and  the 
remainder  of  the  facility  were  made  sufficiently  large  so  that  the  blades 
could  be  instrumented  or  flow  surveys  could  be  made  with  few  spatial 
limitations.  Figure  4  is  a  sketch  of  the  AFRF  which  shows  the  important 
components  and  some  overall  dimensions  of  the  facility.  A  highly  detailed 
description  of  the  system  is  given  in  Reference  [34]. 

The  AFRF  can  most  simply  be  described  as  a  subsonic,  incompressible, 
open  circuit  wind  tunnel  with  an  annular  cross  section.  The  facility  has 
a  hub  radius  of  4.75  inches  (12.1  cm),  a  hub-to-tip  radius  ratio  of  0.442, 
and  was  designed  to  operate  in  the  total  pressure  ratio  range  from  1.0  to 
1.1.  The  overall  length  is  19.67  feet  (6.0  m)  and  consists  of  an  annular 
flow  passage  bounded  at  one  end  by  a  bellmouth  inlet  and  at  the  other  end 
by  an  exhaust  throttle.  The  facility  can  be  operated  with  axial  veloc¬ 
ities  up  to  112  ft/sec  (34.14  m/sec)  with  test  rotor  speeds  up  to 
320  ft/sec  (97.54  m/sec).  The  forward  region,  from  the  inlet  to  a  point 
just  upstream  of  the  test  rotor  blade  drive  motor,  is  bounded  by  a  9.5 
inch  (24.13  cm)  diameter  cylindrical  hub  surface  and  a  21.5  inch  (54.6  cm) 
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Figure  4.  The  Axial  Flow  Research  Fan 
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inside  diameter  cylindrical  outer  casing.  The  drive  unit  for  the  rotor 
is  a  70  horsepower  (52.2  kw)  motor  enclosed  in  the  hub  of  the  AFRF.  The 
auxiliary  fan  drive  motor  and  the  rotor  blade  drive  motor  can  be  independ¬ 
ently  regulated  at  any  speed  up  to  3400  rpm  through  the  use  of  two  Borg- 
Warner  Model  No.  BW1200  Solid  State  Adjustable  Frequency  Drive  inverter 
units.  This  gives  a  maximum  relative  blade  tip  Mach  number  of  0.338. 

The  airflow  through  the  facility  is  controlled  by  the  auxiliary  axial  flow 
fan  and  the  throttle  arrangement  provided  at  the  exit  of  the  facility. 

With  a  zero  steady  lift  rotor  of  the  type  being  used  in  the  initial 
research  experiments,  the  final  adjustment  of  the  throughflov  velocity 
can  be  produced  either  by  adjusting  the  auxiliary  fan  rpm  with  a  fixed 
throttle  setting,  or  by  adjusting  the  throttle  position  with  the  auxiliary 
fan  operating  at  a  fixed  speed. 

The  disturbance  generating  section  is  housed  in  the  two  24-inch 
(61.0  cm)  long  outer  casing  segments  positioned  immediately  downstream 
of  the  bellmouth  inlet.  This  48-inch  (122.0  cm)  long  flowpath  was 
provided  to  permit  the  development  of  distorted  flowfields.  In  the 
absence  of  the  disturbance  producing  screens,  the  inlet  of  the  facility 
was  designed  to  provide  a  steady  uniform  flow  to  the  annulus.  For  studies 
involving  distorted  flow,  screens  are  placed  in  the  inlet  which  produce 
sinusoidal  or  square  circumferential  variations  of  total  pressure.  These 
disturbance  generating  screens  of  the  type  shown  in  Figure  5  provide  a 
capability  to  carry  out  experiments  at  reduced  frequencies  from  0.27  to 
2.46,  based  on  the  mid-span  radius,  a  chord  length  of  six  inches,  and  an 
inlet  relative  flow  angle  of  45  degrees.  The  design  and  performance 
evaluation  of  these  screens  is  presented  in  Reference  [35].  Immediately 
downstream  of  the  screens,  as  shown  by  Figure  4,  is  a  section  of  the 


Figure  5.  Four-Cycle  Disturbance  Producing  Screen 
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annulus  which  serves  as  a  settling  section  where  small  velocity  variations 
and  the  static  pressure  gradient  are  removed.  The  section  between  the 
screen  and  test  rotor  also  isolates  the  rotor  from  the  screen.  This  is 
very  important,  since  a  rotor  will  strongly  alter  the  upstream  flow  in  the 
immediate  vicinity  of  the  rotor. 

The  section  just  downstream  of  the  settling  section  contains  the  test 
rotor-stator  stage.  An  important  feature  of  the  facility  with  regard  to 
this  investigation  is  the  ability  of  the  outer  casing  surrounding  the 
rotor-stator  stage  to  be  rotated  freely  about  the  AFRF  centerline. 

Rotation  of  the  casing  is  accomplished  by  electrically  turning  a  wheel 
driven  worm  gear,  and  it  is  automatically  controlled  by  the  data  acqui¬ 
sition  system.  The  worm  gear  is  driven  by  a  1800  inch-ounce  (130  cm-kg) 
torque  Slo-Syn  Type  SS1800-1007  motor  which  is  controlled  by  a  Slo-Syn 
Preset  Indexer.  This  feature  is  provided  to  enable  circumferential 
surveys  of  the  flow  to  be  performed  by  rigidly  attaching  instrumentation 
to  the  casing  and  indexing  the  casing  in  preselected  increments. 

The  AFRF  design  utilizes  split  outer  casing  and  hub  sections  both 
upstream  and  downstream  of  the  stator.  By  interchanging  sections  of 
different  lengths,  it  is  possible  to  move  the  stator  blade  row  relative 
to  the  rotor  blade  row.  Tests  with  rotor-stator  axial  spacings,  rotor 
trailing  edge  to  stator  leading  edge,  of  3,  6,  9,  or  12  inches  (i.e., 

7.62,  15.24,  22.86,  or  30.48'cm)  are  possible.  For  six-inch  (15.24  cm) 
chord  blades,  these  dimensions  correspond  to  non-dimensional  rotor-stator 
spacings  from  0.5  to  2.0  blade  chord  lengths. 

To  determine  the  influence  of  rotor  geometrical  parameters  on  the 
unsteady  rotor  response,  two  rotor  configurations  were  employed  in  this 
investigation.  The  first  configuration  consisted  of  a  zero  lift, 
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uncambered  rotor  which  was  designed  to  produce  a  zero  circumferential- 
mean  total  pressure  rise  at  its  design  operating  condition  in  a  uniform 
flow.  The  rotor  was  designed  with  this  criterion  primarily  to  enable 
accurate  comparisons  to  be  made  with  unsteady  cascade  theories  which  are 
usually  restricted  to  zero  or  lightly  loaded  blades.  The  characteristics 
of  the  zero  steady  lift  test  rotor  are  shown  in  Figure  6.  For  the  tests 
conducted  in  this  study,  the  rotor  was  used  in  its  six-bladed  configura¬ 
tion.  The  second  rotor  configuration  employed  in  the  tests  consisted  of 
a  nine-bladed  cambered  rotor  with  a  fixed  number  of  blades  and  fixed 
stagger  angle.  Both  the  zero  lift,  uncambered  rotor  and  the  lifting,  or 
cambered,  rotor  were  designed  with  identical  blade  thickness  distribution, 
chord  length,  and  blade  span.  The  characteristics  of  the  cambered  rotor 
are  shown  in  Figure  7. 

For  the  tests  conducted  in  this  study,  the  zero  steady  lift  rotor 
was  assembled  with  six  blades,  which  results  in  a  space-to-chord  ratio 
of  1.35  at  the  mean  radius.  Variations  in  blade  stagger  angle  were 
achieved  by  rotating  the  blades  in  the  hub  and  by  using  different  hub 
spacer  blocks.  In  addition  to  the  nominal  45  degree  stagger  angle  for 
which  the  blades  are  designed,  stagger  angles  of  35  degrees  and  55  degrees 
at  the  mean  radius  were  also  evaluated.  Operation  at  these  stagger  angle 
settings  produces  an  off-design  condition  which  will  result  in  three- 
dimensional  flow  effects.  However,  th8  deviation  of  the  blade  angle  from 
a  design  based  on  that  particular  stagger  angle  is  not  large.  Thus,  the 
three-dimensional  effects  should  be  small  and  the  flow  is  still  justi¬ 
fiably  considered  as  two-dimensional  at  the  mean  radius.  The  data  that 
was  measured  in  the  present  investigation  using  the  six-bladed  rotor  in 
its  35  and  55  degree  stagger  angle  configurations  show  that  the 
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fluctuations  (i.e.,  the  local  peak  values)  of  radial  velocity  were  less 
than  ten  percent  of  the  free-stream  velocity.  Thus,  the  assumption  of 
two-dimensional  flow  at  the  mean  radius  is  valid. 

Downstream  of  the  test  rotor-stator  stage  is  another  section  of  the 
annulus  which  isolates  the  rotor  from  the  effects  of  the  auxiliary  fan  and 
throttle.  As  mentioned  previously,  the  mean  axial  velocity  through  the 
facility  is  controlled  by  the  speed  of  the  auxiliary  fan.  The  test  rotor 
is  driven  by  the  70  horsepower  motor,  so  that  control  of  the  test  rotor 
rpm  setting  is  independent  of  the  auxiliary  fan  rpm  setting.  By  increasing 
the  test  rotor  RPM  while  making  minor  adjustments  to  the  auxiliary  fan  rpm 
seating  to  maintain  the  mean  flow  velocity  at  a  constant  value,  it  is 
possible  to  vary  the  flow  incidence  on  the  test  rotor  blades.  The  test 
rotor  can  thus  be  operated  with  an  off-design  incidence  resulting  in  a 
steady  load  on  the  blades.  This  enables  the  facility  to  be  used  to 
investigate  the  effects  of  incidence,  or  steady  loading,  ranging  all  the 
way  to  stall,  or  the  unsteady  response  of  the  blade  row. 

3.2  Instrumentation 

In  many  complex  flow  fields  such  as  those  encountered  in  turbomachines, 
the  experimental  determination  of  the  three-dimensional  characteristics  of 
the  flow  field  is  frequently  required.  If  space  limitations  or  other 
considerations  make  nulling  techniques  impractical,  five-hole  probes  in 
a  non-nulling  mode  can  be  employed  for  measurements  in  low  speed,  incom¬ 
pressible  flows.  In  this  program,  the  approach  taken  to  measure  the 
response  of  the  rotor  was  to  take  time-mean  flow  measurements  upstream 
and  downstream  of  the  rotor.  To  make  these  measurements,  two  five-hole 
pressure  probes  were  used  in  the  non-nulling  mode.  Five-hole  probes  were 
selected  for  the  measurements  because  they  enable  the  total  pressure, 
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static  pressure,  and  axial,  circumferential,  and  radial  components  of 
velocity  to  be  determined.  The  probes  employed  were  commercially  avail¬ 
able  United  Sensor  prism  type  probes.  The  geometry  of  this  type  of  probe 
offers  easy  insertion  and  mounting  of  the  probe  in  the  AFRF.  An  enlarged 
sketch  of  the  sensing  element  and  a  sketch  of  the  entire  probe  showing 
its  dimensions  is  presented  in  Figure  8.  The  probes  were  positioned  at 
stations  referred  to  as  1-D  and  15,  axially  located  0.137  rotor  blade 
chord  length  upstream  of  the  rotor  and  0.820  chord  length  downstream  of 
the  rotor  at  the  mean  radius,  respectively.  The  probes  used  in  this 
investigation  were  calibrated  in  a  known  flow  field  by  Yocum.  A  brief 
description  of  the  procedure  used  in  the  calibration  of  these  probes  is 
presented  here  for  the  sake  of  completeness. 

The  five-hole  probes  employed  in  the  experimental  program  were  used 
in  a  non-nulling  mode,  which  means  that  they  had  a  fixed  orientation  with 
respect  to  the  casing  during  the  tests.  This  means  that  relationships 
must  be  determined  between  the  measured  pressures  at  the  five  holes  and 
the  true,  local  total  and  static  pressure  or  velocity.  These  desired 
relationships  are  usually  expressed  as  dimensionless  pressure  coefficients, 
which  are  functions  of  the  flow  angularity. 

For  operation  in  the  non-nulling  mode,  the  calibration  character¬ 
istics  must  include  data  that  represent  pressure  differences  in  both  the 
pitch  and  yaw  planes,  see  Figure  9,  as  well  as  differences  between  the 
measured  and  the  true  local  total  and  static  pressures.  The  pressure 
coefficients  representing  these  data  must  be  defined  so  that  they  are 
independent  of  velocity  and  are  a  function  only  of  the  flow  angularity. 

Referring  to  Figure  9,  the  dimensionless  yaw,  pitch,  total  pressure, 
and  static  pressure  coefficients  are  defined  as  follows: 


GEOMETRY  OF  THE  PRISM  PROBES 


VECTOR  RESOLUTION  FOR  A  PROBE 
CALIBRATED  IN  Tl£  PITCH-YAW  MODE 


Figure  9.  Geometry  of  the  Prism  Probes 
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P^  and  P^  »  pressures  measured  in  the  pitch  plane. 

These  coefficients,  which  are  a  function  of  the  yaw  and  pitch  angles, 
completely  describe  the  response  of  the  probe  to  the  flow. 

The  two  flow  angles  and  the  true  total  and  static  pressures  are 
calculated  from  the  five  pressures  measured  from  the  probe  through  the 
use  of  a  detailed  calibration.  The  probes  were  calibrated  in  air  with 
the  use  of  an  open  jet  facility.  To  conduct  the  calibration  of  a  given 
probe,  the  test  velocity  was  maintained  at  a  constant  value.  The  cali¬ 
bration  range  was  +  30  degrees  in  both  the  yaw  and  pitch  planes.  The  probe 
was  positioned  at  one  of  the  predetermined  pitch  angles  and  then  moved  in 
prescribed  increments  through  the  yaw  angle  range.  The  pitch  angle  was 
then  changed,  and  the  probe  was  again  moved  through  the  yaw  angle  range. 

At  each  calibration  point,  the  five  pressures  from  the  probe  as  well  as 
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the  known  total  and  static  pressure  of  the  flow  were  recorded  and  then 
processed  on  an  IBM  1130  computer.  The  data  reduction  program  calculated 
the  four  previously  defined  pressure  coefficients,  the  test  section 
velocity,  and  the  Reynolds  number  based  on  the  probe  tip  diameter.  Each 
probe  was  calibrated  three  times  to  verify  the  repeatability  of  the 
results.  To  reduce  deviations  in  the  calibration  data,  the  results  were 
averaged.  All  calibration  data  were  repeatable  to  within  two  percent  of 
the  reference  dynamic  pressure  when  subjected  to  recalibration  [36]. 

These  data  consistently  described  the  response  of  a  particular  probe, 
ai.d  thus  permit  the  successful  application  of  five-hole  probes  in  a 
non-nulling  mode. 

With  all  of  the  above  information  known,  when  the  probe  is  placed  in 
an  unknown  flow,  the  yaw  and  pitch  angles  of  the  unknown  flow  can  be 
determined  through  the  use  of  the  calibration  data.  After  the  two  angles 
are  known,  the  total  and  static  pressure  coefficients  are  interpolated 
from  the  calibration  data  and  the  two  pressures  can  then  be  calculated. 

The  velocity  is  calculated  with  the  Bernoulli  equation  and  then  resolved 
into  its  axial,  circumferential,  and  radial  components  using  the  known 
flow  angles.  The  derivation  of  the  three  velocity  components  as  well  as 
the  resolution  of  the  flow  angles  for  the  cases  where  the  probe  was 
rotated  in  the  yaw  plane  are  presented  in  Appendix  A. 

Meaningful  calibration  data  should  be  independent  of  the  measured 
quantities.  In  most  five-hole  probe  applications,  velocity  is  the  primary 
parameter  to  be  measured.  It  turns  out  that  the  calibration  curves  of  a 
five-hole  probe  are  usually  weakly  dependent  on  velocity  or  Reynolds 
number.  Therefore,  the  effects  of  Reynolds  number  variations  were 
assessed  by  calibrating  the  five-hole  probes  in  air  over  a  Reynolds 
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number  range  of  2000  to  7000.  The  results  of  these  tests  showed  that  the 
total  pressure,  pitch  coefficient,  and  yaw  coefficient  were  essentially 
unaffected  by  the  Reynolds  number  variation  [36].  On  the  other  hand, 
a  measurable  change  in  the  static  pressure  coefficient  was  observed. 
However,  average  correction  factors  were  incorporated  into  the  data 
reduction  procedures  to  account  for  these  variations  in  the  static 
pressure  coefficient  for  the  range  of  Reynolds  numbers  that  the  probe 
would  experience.  The  correction  could  be  neglected  in  most  cases  with 
the  error  in  the  resulting  velocity  being  approximately  two  percent.  The 
correction  was  included,  however,  to  increase  the  accuracy  of  the  results. 

The  accuracy  of  the  measurements  made  with  the  five-hole  probes  was 
checked  by  comparing  the  results  with  data  obtained  using  other  types  of 
instrumentation.  In  the  absence  of  the  rotor,  static  pressure  gradients 
do  not  exist  far  downstream  of  the  distortion  producing  screen.  Circum¬ 
ferential  surveys  of  the  flow  in  the  AFRF  can  be  conducted  using  a  Kiel 
probe  to  measure  total  pressure  and  a  wall  tap  in  the  casing  to  measure 
static  pressure.  The  velocity  measured  in  this  manner  was  compared  with 
data  obtained  with  the  five-hole  probes.  The  agreement  was  found  to  be 
very  good,  which  demonstrated  the  validity  of  the  five-hole  probe  measure¬ 
ments  and  data  reduction  [37]. 

The  accuracy  of  the  probes  in  the  non-nulling  mode,  as  compared  with 
the  nulling  mode,  was  also  checked  to  assure  that  in  the  non-nulling  mode 
the  probe  did  not  miss  any  of  the  details  of  the  flow.  For  this  check, 
measurements  were  made  with  the  probe  nulled  in  the  plane  perpendicular 
to  its  axis  and  the  results  were  compared  with  the  results  obtained  in 
the  non-nulling  mode.  Once  again  the  agreement  was  found  to  be  quite 


good. 
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A  schematic  sketch  showing  the  entire  experimental  system  is  presented 
in  Figure  10.  The  ten  pressure  outputs  from  the  five-hole  pressure  probes 
were  directly  connected  by  flexible  Tigon  tubing  to  a  scanivalve.  Channels 
1  through  5  on  the  scanivalve  read  the  pressure  differences  from  the  up¬ 
stream  or  location  1-D  probe,  and  channels  6  through  10  read  the  output 
from  the  downstream  or  location  15  probe.  Channel  11  was  connected  to  a 
wall  static  pressure  tap  in  the  outer  casing  of  the  AFRF  located  in  the 
settling  section  of  the  facility,  and  it  read  the  static  pressure  in  the 
axial  plane  of  measurement.  Each  channel  was  recorded  individually  and 
then  automatically  advanced  through  the  use  of  an  automatic  control  on  the 
scanivalve.  The  scanivalve  output  was  directly  fed  into  a  Validyne  Model 
DP15  strain  gaged  differential  pressure  transducer.  This  transducer, 
using  the  ambient  pressure  as  its  reference,  produced  an  electrical  signal 
proportional  to  the  difference  between  the  probe  time-mean  total  pressure 
and  the  ambient  pressure.  The  electrical  signal  was  carried,  as  shown 
schematically  in  Figure  10,  through  a  Validyne  Model  CD15  Carrier 
Demodulator,  a  laboratory-constructed  integrating  digital  voltmeter  and 
data  logger,  and  finally  presented  as  a  voltage  output  on  punched  paper 
tape.  The  paper  tape  output  was  then  interpreted  and  analyzed  on  an  IBM 
1130  computer. 

Before  conducting  measurements  in  the  AFRF,  the  data  sensing  and 
recording  system  was  calibrated  to  determine  its  voltage  output  as  a 
function  of  the  applied  pressure  differential.  The  calibration  was  con¬ 
ducted  using  the  Validyne  pressure  transducer  and  carrier  demodulator,  an 
integrating  digital  voltmeter,  and  a  Meriam  Co.  Micromanometer.  With  the 
micromanometer,  a  range  of  applied  pressures,  readable  to  0.001  inch  of 
water,  was  applied  to  the  system,  and  the  corresponding  voltage  output 


Figure  10.  Schematic  of  Test  Apparatus 
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was  recorded.  Figure  11  presents  the  results  of  this  calibration  and 
determines  the  voltage  output  per  inch  of  water  for  the  particular  trans¬ 
ducer  employed  in  the  experimental  tests.  The  calibration  slope  defined 
the  transducer  constant  which  was  then  used  in  the  computer  analysis  of 
the  data.  Shown  in  Figure  11  are  two  independent  calibrations  of  the 
transducer  conducted  over  a  six-week  period.  These  two  calibrations 
yielded  values  of  transducer  constants  which  differed  by  0.2  percent,  thus 
demonstrating  the  repeatability  of  the  pressure  sensing  system. 

In  addition  to  the  pressure  measurements,  it  was  necessary  to  measure 
the  rotational  speed  of  the  rotor.  A  conventional  system  employing  a 
slotted  disc  rotating  with  the  rotor,  a  light  source,  a  photocell,  and  an 
electronic  counter  in  the  data  acquisition  system  was  used  for  this 
purpose.  With  sixty  slots  in  the  disc  and  a  one-second  sample  period, 
the  counter  reading  gave  the  rpm  of  the  rotor  directly  and  displayed  it 
on  the  data  logger. 

3.3  Experimental  Studies  Conducted 

In  order  to  determine  the  influence  of  rotor  geometrical  parameters 
on  the  unsteady  rotor  response,  two  rotor  configurations  were  employed  in 
this  investigation.  The  two  test  rotors  used  in  the  program  were  a  zero 
lift,  uncambered  rotor  which  has  variable  solidity  and  stagger  angle,  as 
shown  in  Figure  6,  and  a  lifting,  or  cambered,  rotor  with  a  fixed  number 
of  blades  and  stagger  angle,  as  shown  in  Figure  7. 

For  each  geometrical  configuration  of  the  test  rotor,  a  test  was 
conducted  which  involved  the  measurement  of  the  steady  performance  of  the 
rotor  with  uniform  inlet  flow.  This  information  was  necessary  if  the 
performance  degradation  caused  by  the  flow  distortions  was  to  be  deter¬ 
mined.  To  determine  the  rotor  performance,  flow  measurements  were  made 
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upstream  and  downstream  of  the  rotor  over  a  range  of  flow  coefficients. 
Incidence  angles  ranging  from  negative  ten  degrees  up  to  the  point  where 
stall  began  were  evaluated.  The  performance  measurements  were  made  at  the 
mean  radius,  since  the  circumferential  surveys  of  the  distorted  flow  were 
made  at  the  mean  radius  and  the  flow  is  assumed  to  be  two-dimensional 
throughout  the  study. 

The  effects  of  rotor  solidity  (i.e.,  chord-to-spacing  ratio)  were 
studied  using  the  zero  lift  rotor  in  its  six-bladed  configuration.  These 
data  can  be  compared  to  the  results  from  the  rotor  in  its  twelve-bladed 
configuration  obtained  by  A.  M.  Yocum  in  Reference  [37].  For  the  six- 
bladed  test  rotor,  stagger  angle  values  of  35,  45,  and  55  degrees  were 
tested.  At  each  of  these  stagger  angles,  experiments  were  made  using  the 
four-  and  six-cycle  sinusoidal  distortion  screens  and  the  90  degree  and 
180  degree  square  wave  distortion  screens.  The  original  intention  was  to 
test  three  values  of  blade  loading  for  each  stagger  angle  and  distortion 
screen.  The  three  blade  loadings  were  to  correspond  to  the  design  point, 
i.e.,  zero  incidence  angle,  the  condition  of  zero  total  pressure  rise 
(=  APT/l/2pU2),  and  a  condition  nearing  stall.  For  the  twelve-bladed  rotor 
configuration,  it  was  found  that  a  zero  total  pressure  rise  requires  a  small 
positive  incidence  so  that  the  pressure  loss  due  to  the  viscous  effects  is 
balanced  by  the  energy  input  by  the  rotor.  However,  the  steady-state  per¬ 
formance  evaluations  on  the  six-bladed  rotor  revealed  that  the  conditions 
of  zero  incidence  and  zero  total  pressure  rise  were  very  nearly  identical. 
This  is  due  to  the  fact  that  the  rotor  with  only  six  blades  does  not 
experience  as  large  a  total  pressure  loss  as  the  twelve-bladed  rotor. 

Figure  12  shows  a  sample  of  the  procedure  used  to  determine  the 
various  operating  conditions  for  the  six-bladed  rotor.  The  values  used 


“  axial  velocity 


U  =  rotor  blade  speed  at 
mean  radius 

8-^  =  inlet  air  angle 

£  =  stagger  angle 

i  ■  incidence  angle 

r  *  mean  radius  =  7.75  in. 


Known: 

Four-cycle  distortion  screen  ■+■  V  ■  65.911  ft/sec 

A 

Stagger  angle  =  45  degrees 

Desired  incidence  angle  =  8  degrees 

Calculate  Rotor  RPM: 

tan  (C+i)  -  U  /Vv 
in  a 

U  *  V  tan  (£+i)  ■  (65 . 91 1  ft/sec)  tan  53  degrees 
m  a 

U  =  87.47  ft/sec 
in 

rpm  »  (87.47  ft/sec)  (12  in. /ft)  (60  sec/min) 

(2tt  rad/rev)  (7.75  in. /rad) 

RPM  -  1293 


Figure  12.  Calculation  of  Test  Conditions  for  the 
Six-Bladed  Uncarabered  Rotor 


for  the  axial  velocity  are  circumferential-mean  values  which  were  deter¬ 
mined  by  Yocum  [37].  They  were  determined  by  making  circumferential 
surveys  of  the  distorted  flow  in  the  absence  of  the  rotor.  Thus,  a  value 
for  the  axial  velocity  was  predetermined  for  each  distortion  screen  used. 
Next,  the  desired  value  of  circumferential  mean  incidence  angle  was 
selected.  Then,  by  knowing  the  axial  velocity,  incidence  angle,  and 
stagger  angle  a  value  of  blade  speed  at  the  mean  radius  is  calculated 
and  converted  to  revolutions  per  minute.  This  rotor  rpm  was  the  variable 
parameter  which  defined  the  different  operating  conditions  for  the  six- 
bladed  rotor. 

A  tabulated  summary  of  the  test  conditions  which  were  conducted  with 

the  six-bladed,  non- lifting  rotor  appears  in  Table  B.l  of  Appendix  B. 

Presented  are  the  type  of  distortion  tested,  the  stagger  angle  of  the 

blades,  the  rotor  rpm,  the  average  flow  coefficient,  V„/U  ,  based  on  the 

Xm 

blade  speed  at  the  mean  radius,  the  average  pressure  rise  coefficient, 
APT/l/2pV^2,  and  a  circumferentially  averaged  value  of  incidence  angle  at 
the  rotor  mean  radius. 

The  next  series  of  experiments  performed  in  this  investigation 
utilized  a  cambered  test  rotor  with  nine  blades  and  a  fixed  stagger  angle 
of  50  degrees  at  the  mean  radius.  Essentially  two  effects  were  studied 
with  the  cambered  rotor.  First  of  all,  the  effect  of  rotor-stator  spacing 
on  the  response  of  the  rotor  was  examined.  Eight  cambered  stator  blades 
of  the  same  design  as  the  rotor  blades  were  installed  in  the  Axial  Flow 
Research  Fan  for  this  purpose.  In  the  first  series  of  tests,  the  stator 
blade  row  was  located  twelve  inches  (30.48  cm),  or  two  rotor  blade  chord 
lengths,  downstream  of  the  rotor  trailing  edge.  Then,  to  look  at  the 
effects  of  rotor-stator  spacing,  another  series  of  tests  was  performed 
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using  the  same  operating  conditions  as  in  the  first  series.  However,  for 
these  tests,  the  stator  row  was  moved  to  a  position  three  inches  (7.62  cm) 
or  one-half  blade  chord  length,  downstream  of  the  rotor  trailing  edge. 

For  all  of  the  tests  just  mentioned,  the  one-,  two-,  four-,  and  six-cycle 
sinusoidal  distortions  and  the  90  degree  and  180  degree  square  wave  dis¬ 
tortions  were  used.  Three  values  of  blade  loading,  one  at  the  design 
point,  one  below  design,  and  one  above  design,  were  evaluated. 

The  final  series  of  tests  with  the  cambered  rotor  studied  the  effect 
of  spanwise  location  on  the  inlet  distortion,  that  is,  to  see  what  the 
distortion  flow  field  looks  like  at  radial  locations  other  than  the  mean 
rrdius  of  the  rotor  blades.  Four  radial  locations  were  tested;  two 
locations  were  between  the  hub  and  the  mean  radius,  and  two  were  between 
the  mean  radius  and  the  tip  of  the  blade.  At  each  radial  location,  the 
two-  and  four-cycle  sinusoidal  distortions  and  the  90  degree  and  180 
degree  square  wave  distortions  were  tested,  and  only  one  incidence  angle 
was  used.  A  tabulated  summary  of  the  operating  conditions  used  for  all 
the  tests  conducted  with  the  nine-bladed  cambered  rotor  appears  in  Tables 
B.2  through  B.7  of  Appendix  B. 

It  should  be  noted  that  a  slightly  different  approach  was  used  to 

determine  the  operating  conditions  for  the  cambered  rotor  as  opposed  to 

the  zero  lift  rotor.  Figure  13  shows  the  method  used  to  determine  the 

blade  loadings  used  for  the  tests  on  the  cambered  rotor.  In  the  tests 

with  the  non-lifting  rotor,  the  values  of  incidence  angle  to  be  used  were 

selected  from  the  steady  performance  data,  and  the  rotor  rpm  was  then 

calculated.  But  for  the  cambered  rotor,  values  of  flow  coefficient, 

V  /U  ,  were  selected  from  the  steady  performance  data.  Then,  by  knowing 
a  m 

the  value  of  axial  velocity,  Vx>  to  be  used  for  each  distortion  screen, 
the  rotor  rpm  was  calculated  to  give  the  proper  flow  coefficient. 
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V  =  axial  velocity 

A 

U  *  rotor  blade  speed  at  mean  radius 
m 

3^  =  inlet  air  angle 

8  '  =  inlet  blade  angle 

E,  =  stagger  angle 

i  =  incidence  angle 

r  =  mean  radius  =  7.75  in. 
m 


Known: 

Two-cycle  distortion  screen  -+  VY  »  67.533  ft/sec 

A 

Desired  flow  coefficient  =  0.586 
Calculate  Rotor  RPM: 

VU.  -  ‘~531 i,f‘/«.c  ■  0.586 
m 

U  »*  115.24  ft/sec 
m 

RPM  m  (115.24  ft/sec)  (12  in. /ft)  (60  sec/min) 
(2tt  rad/rev)  (7.75  in. /rad) 

RPM  -  1704 

Then, 

tan  (8^+i)  -  Um/Vx  -  (1.706) 
i  *  tan  (1.706)  —  56  degrees 
i  *  3.623  degrees 


Figure  13.  Calculation  of  Test  Conditions  for  the 
Nine-Bladed  Cambered  Rotor 


CHAPTER  IV 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 

4. 1  Rotor  Steady  Performance 

When  describing  the  performance  of  an  axial  flow  fan  or  compressor 

in  a  uniform  inflow,  the  performance  characteristic  of  primary  interest 

is  the  relationship  between  the  pressure  rise  and  the  inlet  flow  conditions 

and  rotor  speed.  For  an  incompessible  flow,  this  desired  relationship  is 

adequately  described  by  relating  the  total  pressure  rise  coefficient,  ¥, 

defined  as  (P„  .  -  P„  .  ,  )/l/2pU2  to  the  flow  coefficient,  $,  defined 

1  exit  T  inlet 

as  Vx/U. 

For  a  flow  with  a  zero  circumferential  velocity  component  upstream 
of  the  rotor,  and  assuming  there  are  no  losses,  a  relationship  is 
described  between  ¥  and  $  from  Euler's  pump  equation.  This  equation, 
which  shows  a  direct  relationship  between  ¥  and  $,  is  as  follows: 

¥  -  2  -  2$tan&2  .  (32) 

For  a  particular  rotor  configuration,  B2  will  depend  on  the  flow  incidence 
and  thus  will  also  be  related  to  the  flow  coefficient,  $.  Furthermore, 
the  losses  through  the  blade  row,  which  are  neglected  in  Equation  (32), 
will  depend  primarily  on  the  incidence  or  the  flow  coefficient.  For  these 
reasons,  ¥  is  expected  to  be  a  function  of  $  alone,  which  justifies 
expressing  the  rotor  performance  in  terms  of  the  total  pressure  rise 
coefficient  and  the  flow  coefficient  [37]. 

The  experimental  data  of  ¥  versus  $  measured  at  the  mean  radius  is 
shown  in  Figure  14  for  the  six-bladed,  zero  steady  lift  rotor  configu¬ 
rations  having  blade  stagger  angles  of  33,  45,  and  55  degrees,  and 
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Figure  15  shows  the  data  for  the  cambered  rotor  configuration  having  a 
fixed  stagger  angle  of  50  degrees.  For  these  data,  which  determine  the 
steady  performance  of  the  rotors,  the  flow  rate  was  set  to  give  an  average 
axial  velocity  of  67  ft/sec  (20.4  m/sec).  This  was  done  to  enable  the 
probes  to  operate  in  approximately  the  same  Reynolds  number  range  at  all 
times.  Measurements  were  made  only  at  the  rotor  mean  radius.  The  rotor 
blade  speed  was  then  adjusted  to  give  a  range  of  incidence  angles  in  two 
degree  increments  from  negative  ten  degrees  up  to  the  point  where  the 

rotor  began  to  stall.  This  range  of  incidence  angles,  and  thus  flow 

« 

coefficients,  was  sufficient  to  completely  describe  the  steady  performance 

characteristics  of  the  two  rotors. 

An  important  characteristic  of  the  data  presented  in  Figures  14  and 

15  is  the  slope  of  the  curves.  In  a  distorted  flow  field,  the  distorted 

and  undistorted  sectors  of  the  flow  field  have  different  axial  velocities, 

and,  thus,  operate  at  different  values  of  flow  coefficient.  From  Figure 

14,  it  can  be  seen  that  the  distorted  and  undistorted  sectors  experience 

different  amounts  of  total  pressure  rise  (the  distorted  sector  would  see 

a  higher  total  pressure  rise  than  the  undistorted  sector).  Therefore, 

from  a  quasi-steady  viewpoint  the  slope  of  these  curves  determines  the 

attenuation  of  a  total  pressure  distortion.  It  is  seen  from  Figure  14 

that  the  slope  of  the  performance  curves  becomes  steeper,  i.e.,  more 

negative,  with  increasing  stagger  angle  for  the  uncambered  rotor.  This 

result  can  be  expected  by  examining  Equation  (32).  By  assuming  that  f?2 

dH* 

is  constant  and  differentiating  Equation  (32),  it  is  found  that  ^  * 

-2tan82-  For  the  rotor  with  no  camber,  except  for  a  small  deviation  angle, 

d'P 

@2  is  approximately  equal  to  the  stagger  angle  and,  therefore,  ^  will 
become  more  negative  with  higher  stagger  angles.  This  does  not  hold  true, 
however,  for  the  case  of  the  cambered  rotor.  It  is  observed  from  Figure  15 
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that  the  slope  of  the  curve  for  the  cambered  rotor  with  a  stagger  angle  of 
50  degrees  is  more  negative  (steeper)  than  the  slope  of  the  curve  for  the 
uncambered  rotor  with  a  stagger  angle  of  55  degrees.  For  the  cambered 
rotor,  the  exit  blade  angle,  and  the  stagger  angle  are  no  longer 
approximately  equal.  In  fact,  for  the  nine-bladed  cambered  rotor  used  in 
this  investigation,  the  angle  0^  is  equal  to  56  degrees.  This  accounts 
for  the  steeper  slope  for  the  cambered  rotor  despite  its  smaller  stagger 
angle,  and  the  expected  result  seen  from  the  differentiation  of  Equation 
(32)  is  still  valid. 

The  next  important  concideration  to  be  made  from  the  steady  perfor¬ 
mance  data  is  the  amount  of  swirl  put  into  the  flow  by  the  rotor.  This 
circumferential  component  of  velocity,  or  swirl,  is  important  because  it 
allows  the  energy  added  to  the  flow  by  the  rotor,  or  the  ideal  total 
pressure  rise,  to  be  calculated.  The  ideal  total  pressure  rise  is  required 
to  determine  the  total  pressure  losses  across  the  rotor. 

In  Figures  16  and  17,  experimental  values  of  /U  measured  at 

the  mean  radius  are  presented  as  a  function  of  flow  coefficient,  Vx/U> 
for  the  uncambered  and  cambered  rotors,  respectively.  It  is  observed 
from  Figures  16  and  17  that  the  data  are  almost  perfectly  linear.  This 
linear  nature  is  not  surprising,  however,  since  the  flow  through  a  blade 
row  usually  leaves  nearly  tangent  to  the  blade  trailing  edge  (except  for 
a  small  deviation  angle),  unless  large  blade  loadings  cause  flow  separa¬ 
tion.  With  the  flow  leaving  the  blades  at  a  constant  relative  angle,  it 
can  easily  be  shown  that  the  swirl  should  vary  linearly  with  flow  coeffi¬ 
cient.  At  low  values  of  flow  coefficient  which  correspond  to  large  blade 
loadings,  a  slight  deviation  from  the  linear  behavior  of  the  curves  is 
seen  in  Figures  16  and  17.  This  is  the  point  at  which  the  rotor  begins 
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Figure  17.  A  Comparison  of  V@  exit^  versus  Vjj/U  Obtained 
Experimentally  in  Uniform  Flow  at  the  Mean 
Radius  for  the  Cambered  Rotor 
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to  operate  in  the  stall  region  and  local  flow  separation  on  the  blades 
begins  to  occur. 

To  complete  the  study  of  the  steady  performance  of  the  rotors,  the 
data  of  the  ideal  and  actual  total  pressure  rise  must  be  interpreted. 

From  these  data  the  total  pressure  losses  are  also  obtained.  The  steady 
losses  are  needed  not  only  to  evaluate  the  rotor  performance  in  a  uniform 
flow,  but  are  also  needed  as  a  reference  to  evaluate  the  losses  associated 
with  the  distorted  flow.  The  ideal  total  pressure  rise  is  obtained  from 
Euler's  pump  equation  and  is  expressed  non-dimensionally  as  follows: 

T  ideal  '  0  exit  0  inlet; 

^ ideal  =  l/wF"  =  - U -  *  (33) 

As  stated  previously  and  as  seen  from  Equation  (33),  the  ideal  total 
pressure  rise  is  a  measure  of  the  energy  added  to  the  flow  by  the  rotor, 
and  it  depends  upon  the  amount  of  swirl  put  into  the  flow  by  the  rotor. 
Upstream  of  the  rotor,  the  inlet  flow  has  no  angular  momentum,  and,  there¬ 
fore,  has  no  component  of  circumferential  velocity  or  swirl.  Therefore, 
the  inlet  component  of  circumferential  velocity  was  assumed  to  be  zero  in 
the  calculations  of  ideal  total  pressure  rise  from  Equation  (33).  In 
reality,  since  the  upstream  probe  is  located  very  close  to  the  leading 
edge  of  the  rotor,  there  is  a  small  component  of  circumferential  velocity 
observed  at  the  upstream  location  which  is  induced  by  the  rotor.  However, 
for  the  calculations  of  ideal  pressure  rise,  this  small  value  of  inlet 
swirl  is  assumed  to  be  zero. 

The  actual  pressure  rise  coefficient,  ¥,  is  simply  the  difference 
between  the  total  pressures  measured  upstream  and  downstream  of  the  rotor. 
The  difference  between  the  ideal  and  actual  pressure  rise  represents  the 


total  pressure  losses  across  the  rotor.  A  loss  coefficient,  Cp 
is  defined  in  this  manner  as  it  is  expressed  in  the  following 
equation: 
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The  value  of  loss  coefficient  must  necessarily  be  greater  than  zero 
since  the  ideal  pressure  rise  represents  the  highest  value  of  pressure 
rise  which  can  be  achieved  for  the  given  operating  conditions,  since  any 
losses  across  the  rotor  are  neglected.  However,  it  was  initially  found 
that  in  some  cases  of  the  cambered  rotor  losses,  the  actual  total 
pressure  rise  coefficient  exceeded  the  ideal  total  pressure  rise 
coefficient.  After  a  thorough  investigation  of  the  data,  the  problem  was 
determined  to  be  in  the  components  of  the  measured  circumferential  velocity 
which  were  used  to  determine  the  ideal  pressure  rise  coefficients.  A 
probable  cause  for  the  problem  was  that  the  five-hole  probes  were  not 
properly  aligned  perpendicular  to  the  axial  flow  direction  in  the  yaw 
plane. 

Independent  calibrations  of  the  disturbance  generating  screens  were 
performed  at  certain  intervals  during  this  experimental  investigation. 

These  screen  calibrations  involve  making  circumferential  surveys  of  the 
distorted  flow  produced  by  the  screens  in  the  absence  of  the  rotor.  Since 
there  is  no  rotor  present  to  put  any  turning  into  the  flow,  the  circum¬ 
ferential  component  of  velocity  is  expected  to  be  zero  at  both  upstream 
and  downstream  measuring  locations.  The  periods  at  which  the  screen 
calibrations  were  carried  out  were  carefully  logged,  and  the  tests  were 
made  with  the  probes  aligned  as  they  were  for  the  present  investigation. 
Thus,  it  was  proposed  that  an  examination  of  the  screen  calibration  data 


would  yield  a  correction  wich  could  be  applied  to  the  yaw  angle  alignment 
of  the  probes.  The  data  revealed  that  significant  components  of  circum¬ 
ferential  velocity  occurred  at  both  measurement  locations  used  in  the  tests 
made  with  the  cambered  rotor  and  at  the  upstream  measurement  location  used 
in  the  studies  with  the  non-lifting  rotor.  These  values  of  circumferential 
velocity  determined  from  tests  without  the  rotor  installed  were  found  to 
be  from  two  to  five  percent  of  the  axial  free-stream  velocity,  and  they 
were  consistent  for  a  number  of  tests  done  at  values  of  axial  velocity 
corresponding  to  those  used  in  this  study.  Based  on  the  relationship 

tana.  =  VQ/V„,  the  amount  of  deviation  from  the  normal  direction  in  the 
Z  o  X 

yaw  angle  of  the  probes  could  be  found.  For  the  tests  done  with  the  nine- 
bladed  cambered  rotor,  the  upstream  and  downstream  probes  were  out  of 
alignment  by  1.23  and  2.75  degrees,  respectively.  For  the  tests  done 
with  the  six-bladed  uncarabered  rotor,  the  downstream  probe  was  aligned 
properly  and  the  upstream  probe  deviated  from  the  normal  direction  by  3.49 
degrees.  It  should  be  noted  that  a  correction  of  this  type  need  not  be 
applied  to  the  upstream  probe  in  order  to  calculate  the  ideal  total 
pressure  rise  since  the  inlet  component  of  circumferential  velocity  was 
assumed  to  be  zero. 

This  method  of  arriving  at  a  correction  factor  to  be  applied  to  the 
data  proved  to  be  clean  and  acceptable  for  two  reasons.  First,  the  screen 
calibrations  were  made  at  times  coinciding  with  the  tests  made  in  the 
present  investigation,  and  the  experimental  apparatus,  particularly  the 
alignment  of  the  probes,  remained  unchanged.  Secondly,  the  rotor  was  not 
present  to  disrupt  the  flow,  making  the  determination  of  an  angular 
deflection  in  the  yaw  plane  quite  simple  and  accurate. 

The  effect  of  the  yaw  angle  correction  was  minimal  on  the  axial 
velocity,  radial  velocity,  and  static  pressure;  only  the  circumferential 


component  of  velocity  was  significantly  altered.  The  total  pressure  was 
completed  unchanged  because  the  total  pressure  is  determined  from  the 
center  hole  on  the  five-hole  probe,  and  the  readings  on  the  center  hole 
will  not  change  over  a  range  of  ±10  degrees  in  the  yaw  plane  (as  deter¬ 
mined  from  the  calibration  of  the  five-hole  probe  [37]). 

Figures  18,  19,  and  20  present  the  ideal  (with  yaw  angle  corrections) 
and  actual  total  pressure  rise  coefficient  and  the  total  pressure  loss 
coefficient  as  a  function  of  incidence  angle  for  the  six-bladed  uncambered 
rotor  configurations  having  blade  stagger  angles  of  35,  45,  and  55  degrees, 
respectively.  Figure  21  shows  the  three  pressure  coefficients  as  a 
function  of  incidence  for  the  nine-bladed  cambered  rotor  having  a  blade 
stagger  angle  of  50  degrees.  Typically,  the  losses  are  lowest  at  the  low 
absolute  values  of  incidence.  As  the  incidence  increases  in  both  the 
positive  and  negative  directions  the  losses  also  increase.  The  loss  data 
for  the  non-lifting  rotor  with  a  stagger  angle  of  45  degrees  and  for  the 
cambered  rotor  with  a  50  degree  stagger  angle  show  similar  trends,  but  a 
slight  dip  in  the  loss  curve  is  evident  at  higher  incidence  angles.  For 
the  45  degree  stagger  configuration,  the  dip  occurs  at  approximately  twelve 
degree  incidence;  for  the  50  degree  stagger  configuration,  the  dip  occurs 
at  a  much  lower  value  of  approximately  two-degree  incidence  angle.  The 
loss  curve  for  the  55  degree  stagger  angle  appears  to  begin  to  dip  slightly 
at  an  incidence  angle  of  roughly  eight  degrees,  and  for  the  35  degree 
stagger  configuration,  the  dip  begins  to  occur  slightly  at  an  incidence  of 
15  degrees. 

The  dip  in  the  loss  data  which  is  most  evident  in  Figures  19  and  21 
for  the  45  degree  and  50  stagger  configurations  was  also  observed  in  the 
steady  performance  data  measured  with  the  twelve-bladed  uncambered  rotor 
presented  in  Reference  [37],  As  pointed  out  in  Reference  [37],  the  dip 
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Figure  21.  The  Ideal  and  Actual  Total  Pressure  Rise  Coefficient  and 
the  Total  Pressure  Loss  Coefficient  as  a  Function  of 
Incidence  Angle  for  the  Cambered  Rotor  with  a  Blade 
Stagger  Angle  of  50  Degrees 
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in  the  loss  data  was  not  important  enough  to  warrant  the  additional  experi¬ 
mentation  which  would  be  required  to  find  the  cause.  However,  several 
possible  causes  were  proposed  to  explain  the  existence  of  the  dip  in  the 
loss  curves.  For  the  nine-bladed  rotor  with  camber,  the  dip  in  the  loss 
curve  can  be  explained  by  the  fact  that  the  dip  occurs  at  a  point  which 
corresponds  to  the  minimum- loss  incidence  angle  for  a  circular  arc  blade 
blade  in  a  two-dimensional  cascade  [38].  It  can  also  be  seen  in 
both  Figures  19  and  21  that  the  decrease  in  the  losses  at  the  higher 
incidence  angles  is  due  primarily  to  a  change  in  the  ideal  pressure  rise, 
which  is  calculated  from  the  measured  circumferential  velocity  components. 
Referring  to  Figures  16  and  17,  where  the  measured  circumferential  velocity 
components  at  the  exit  of  the  rotor  are  plotted,  a  small  dip  is  evident 
in  the  data  of  Vq  exitVu  versus  0  at  the  values  of  $  corresponding  to  the 
incidences  where  the  dip  is  found  in  the  loss  curves.  A  change  in  the 
Vg  component  caused  by  experimental  error  could  be  small  and  still  cause 
the  observed  dip  in  the  loss  data.  But  since  the  dip  is  observed  for  both 
the  45  and  50  degree  stagger  configurations  and  also  appears  to  begin  to 
occur  for  the  35  and  55  degree  stagger  configurations,  it  is  unlikely  that 
this  phenomenon  can  be  attributed  to  experimental  error. 

As  the  Incidence  angle  increases,  the  pressure  gradient  across  the 
blade  row  becomes  more  adverse  and  the  losses  in  the  boundary  layer  on  the 
blades  become  larger.  This  trend  will  eventually  lead  to  separation  of 
the  flow  and  a  very  rapid  increase  in  the  pressure  losses.  Therefore,  a 
proposed  explanation  for  the  observed  behavior  in  the  loss  data  is  that 
the  flow  initially  undergoes  a  laminar  separation  causing  an  increase  in 
losses,  but  at  higher  incidence,  the  flow  either  reattaches  or  the  laminar 
separation  is  eliminated  completely  due  to  the  increased  turbulence  with 
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the  larger  relative  velocities  at  higher  rotor  speeds.  Reference  [37] 
suggests  that  laminar  separation  is  not  likely  to  have  occurred  since  the 
Reynolds  numbers  corresponding  to  the  operating  conditions  at  which  the 
dips  occur  in  the  loss  data  are  not  within  the  Reynolds  number  range  where 
laminar  separation  may  take  place.  The  lower  Reynolds  number  limit  where 
laminar  separation  may  occur  is  between  3  x  105  and  2  x  105,  where  the 
Reynolds  number  is  based  on  the  chord  length  of  the  blades  and  the  relative 
flow  velocity.  In  Figure  18,  for  example,  the  losses  first  become  large 
at  an  incidence  angle  of  approximately  9.5  degrees.  The  Reynolds  number 
corresponding  to  this  operating  condition  is  2.9  x  105,  which  suggests 
that  laminar  separation  is  not  likely  to  have  occurred.  Then,  since 
laminar  separation  does  not  occur,  another  possible  explanation  for  the 
deflection  in  the  loss  curve  is  turbulent  separation,  but  verification 
of  turbulent  separation  would  require  additional  experimentation. 

4.2  Distortion  Measurements  With  the  Rotor  Installed 

The  major  phase  of  the  experiments  consisted  of  measuring  the  dis¬ 
torted  flow  field  upstream  and  downstream  of  the  operating  rotor.  As 
discussed  in  Section  3.3,  a  parametric  study  was  conducted  in  which  the 
response  of  the  rotor  to  an  inlet  distortion  was  obtained  as  a  function  of 
various  design  and  operating  conditions.  These  distortion  measurements 
could  then  be  compared  to  the  measurements  made  by  Yocum  in  Reference 
[37]  of  the  distorted  flow  field  in  the  absence  of  the  rotor.  This  is 
an  important  comparison  since  it  shows  how  the  rotor  alters  the  inlet  and 
exit  flows.  Aside  from  comparing  the  distortion  measurements  made  with 
and  without  the  rotor  installed,  this  section  also  demonstrates  the  type 
of  data  which  was  obtained  and  explains  some  of  the  more  noticeable 
features  of  the  distortion  shapes  both  upstream  and  downstream  of  the 
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rotor.  Due  to  the  excessively  large  amount  of  data  obtained  in  this 
study,  it  was  not  practical  to  document  all  of  the  data  here.  A  complete 
collection  of  all  the  data  can  be  found  in  Reference  [39],  which  includes 
all  the  circumferential  survey  measurements  for  both  the  six-bladed 
uncambered  and  nine-b laded  cambered  rotors.  In  Appendix  B,  tables  of 
the  run  conditions  for  all  the  tests  conducted  in  this  investigation  are 
presented. 

The  results  of  the  distortion  measurements  made  without  the  rotor 
installed  for  a  four-cycle  sinusoidal  distortion  are  presented  in 
Reference  [37].  Circumferential  flow  surveys  were  conducted  at  the 
same  axial  locations  as  the  experiments  when  the  rotor  was  installed 
(i.e.,  stations  1-D  and  15),  thus  providing  a  direct  reference  with  which 
to  determine  the  changes  caused  by  the  rotor.  As  mentioned  previously, 
the  determination  of  the  flow  field  characteristics  without  the  rotor  is 
very  important  since  the  rotor  will  affect  the  flow  field.  A  true  measure 
of  the  velocity  distortion  attenuation  requires  a  knowledge  of  the  flow 
in  the  absence  of  the  rotor. 

The  results  of  the  measurements  made  in  the  absence  of  the  rotor  are 
presented  in  Reference  [37]  as  graphs  of  the  dimensionless  axial,  circum¬ 
ferential,  and  radial  velocities  as  well  as  the  dimensionless  total  and 
static  pressures  as  a  function  of  the  circumferential  location.  It  is 
seen  from  these  results  that  no  changes  occur  in  the  flow  between  the  two 
measurement  locations  and  there  are  also  no  static  pressure  variations. 

For  these  reasons,  it  can  be  concluded  that  any  flow  redistribution  caused 
by  the  screen  occurs  upstream  of  the  measurement  locations,  and  the 
viscous  effects  due  to  the  shearing  nature  of  the  flow  are  not  important 
over  the  short  distance  between  the  measurement  locations  [37]. 
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Therefore,  any  change  observed  when  the  rotor  is  installed  will  be  due 
to  the  rotor  alone. 

4.2.1  Results  for  the  Six-Bladed  Zero  Steady  Lift  Rotor.  Typical 
results  from  the  measurements  made  of  the  distorted  flow  field  with  the 
six-bladed  uncambered  rotor  operating  in  a  four-cycle  distortion  appear 
in  Figures  22  through  27,  wh^re  the  axial,  circumferential,  and  radial 
velocities,  total  and  static  pressures,  and  incidence  angle  are  plotted, 
respectively.  The  data  for  the  two  axial  locations  are  presented  on  the 
same  figure,  location  1-D  being  0.137  chord  length  upstream  of  the  rotor 
and  location  15  being  0.82  chord  length  downstream  of  the  rotor.  These 
figures  reveal  the  usual  characteristics  observed  in  distortion  measure¬ 
ments.  A  comparison  of  the  axial  velocity  with  and  without  the  rotor 
present  shows  that  the  axial  velocity  distortion  was  attenuated  upstream 
of  the  rotor.  Downstream  of  the  rotor,  the  axial  velocity  distortion 
is  seen  to  have  undergone  a  further  reduction  in  amplitude.  It  is  seen 
that  the  rotor  creates  a  variation  in  the  circumferential  component  of 
velocity  upstream  of  the  rotor.  Downstream  of  the  rotor,  the  circum¬ 
ferential  velocity  is  nearly  uniform.  The  radial  velocity  distribution 
presented  in  Figure  24  shows  that  the  radial  component  of  velocity  is 
less  than  ten  percent  of  the  axial  free-stream  velocity  for  this  condition 
of  large  incidence  angle,  8.5  degrees.  This  demonstrates  the  basis  of  the 
assumption  made  throughout  the  study  that  the  flow  field  is  two-dimensional. 
At  lower  values  of  incidence  angle,  the  radial  velocities  were  less.  It 
is  observed  in  Figure  26  that  a  static  pressure  gradient  arises  at  the 
upstream  location.  This  is  due  to  the  presence  of  the  rotor  since  a 
gradient  in  static  pressure  at  the  upstream  location  was  not  observed 
without  the  rotor  operating.  Downstream  of  the  rotor,  as  seen  in  Figure 
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Figure  22.  Dimensionless  Axial  Velocity  Variation  for  the 
Four-Cycle  Distortion  with  the  Six-Bladed 
Uncambered  Rotor 
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Figure  23.  Dimensionless  Circumferential  Velocity  Variation 
for  the  Four-Cycle  Distortion  with  the  Six-Bladed 
Uncambered  Rotor 
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Figure  24.  Dimensionless  Radial  Velocity  Variation  for  the 
Four-Cycle  Distortion  with  the  Six-Bladed 
Uncambered  Rotor 
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on  the  Blades  for  the  Four-Cycle  Distortion 
with  the  Six-Bladed  Uncambered  Rotor 
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26,  the  static  pressure  is  constant.  A  comparison  of  the  total  pressure 
distortion  data  with  and  without  the  rotor  operating  reveals  that  the 
total  pressure  distortion  upstream  of  the  rotor  is  basically  unchanged 
since  no  work  has  been  done  on  the  flow.  However,  as  the  flow  passes 
through  the  rotor,  the  amplitude  of  the  total  pressure  distortion  is 
significantly  reduced.  Finally,  Figure  27  shows  the  incidence  angle  on 
the  blades  as  a  function  of  circumferential  location.  From  a  quasi-steady 
viewpoint,  the  local  incidence  determines  the  fluctuating  forces  on  the 
blades . 

Presented  in  Figures  28  through  33  are  the  distortion  data  obtained 
in  the  AFRF  for  a  90  degree  square  distortion  with  the  six-bladed  rotor 
operating.  Once  again,  the  dimensionless  velocity  and  pressure  components 
as  well  as  the  local  blade  incidence  angle  are  plotted  as  functions  of 
circumferential  location.  In  Reference  [37],  the  results  are  included  of 
the  distortion  measurements  made  with  the  90  degree  square  distortion  with¬ 
out  the  rotor  installed.  It  is  found  from  a  comparison  of  these  data  with 
Figures  28,  29,  and  32  that  the  changes  in  the  axial  velocity,  circum¬ 
ferential  velocity,  and  static  pressure  caused  by  the  rotor  are  the  same 
as  those  just  discussed  for  a  four-cycle  sinusoidal  distortion.  However, 
the  total  pressure  distortion  behaves  differently  since  it  is  found  from 
a  comparison  of  Figure  31  and  the  data  in  Reference  [37]  that  the  magnitude 
of  the  distortion  upstream  of  the  rotor  has  increased.  Although  the  change 
was  small,  no  change  was  expected  because,  from  a  two-dimensional  view¬ 
point,  the  total  pressure  at  the  mean  radius  can  only  be  changed  by  doing 
work  on  the  fluid  or  if  a  change  in  the  mixing  losses  occurs.  The  rotor 
obviously  cannot  do  work  on  the  fluid  (or  remove  energy)  upstream  from 
itself.  It  is  possible  that  the  mixing  losses  could  increase  due  to  the 
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Figure  28.  Dimensionless  Axial  Velocity  Variation  for  the 
90  Degree  Square  Distortion  with  the  Six-Bladed 
Uncambered  Rotor 
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Figure  29.  Dimensionless  Circumferential  Velocity  Variation 
for  the  90  Degree  Square  Distortion  with  the  Six- 
Bladed  Uncambered  Rotor 
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flow  redistribution  which  takes  place.  However,  the  symmetry  of  the  dis¬ 
tortion  and  the  magnitude  of  the  change  makes  this  unlikely.  Therefore, 
as  pointed  out  in  Reference  [37],  the  most  likely  cause  of  the  observed 
change  in  the  total  pressure  distortion  amplitude  is  a  radial  total 
pressure  gradient  produced  by  the  screens,  followed  by  a  radial  shift  in 
the  streamlines  due  to  the  flow  redistribution  which  takes  place  near  the 
rotor.  This  radial  shift  would  apparently  invalidate  the  two-dimensional 
flow  assumption  made  throughout  this  study,  but  the  change  in  the  total 
pressure  distortion  is  sufficiently  small  so  that  the  radial  shift  may  be 
neglected. 

4.2.2  Results  for  the  Nine-Bladed  Cambered  Rotor.  A  series  of 
distortion  measurements  was  also  made  using  a  nine-bladed  cambered 
rotor  with  a  fixed  stagger  angle  of  50  degrees  at  the  mean  radius. 

In  order  to  determine  the  effect  of  the  rotor  on  the  inlet  distortion  at 
different  radial  or  spanwise  locations  along  the  blade,  tests  were  con¬ 
ducted  with  the  probes  located  at  the  following  radial  positions:  9.333 
inches  (23.708  cm),  8.417  inches  (21.381  cm),  6.583  inches  (11.829  cm), 
and  5.667  inches  (14.396  cm).  The  mean  radius  of  the  blade  is  7.75  inches 
(19.685  cm),  so  the  tests  were  made  at  two  radii  greater  than  the  mean 
radius  and  two  radii  less  than  the  mean  radius.  As  before,  the  probes 
were  axially  located  at  stations  1-D  and  15. 

Figures  34,  35,  36,  and  37  show  the  axial  velocity,  circumferential 
velocity,  total  pressure,  and  static  pressure  variations,  respectively, 
for  the  measurements  made  for  a  four-cycle  distortion  with  the  probes 
located  at  the  radial  position  closest  to  the  outer  casing  wall.  In 
order  to  examine  the  effects  of  different  radial  locations,  Figures  38 
through  41  show  the  velocity  and  pressure  distributions  for  approximately 
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Figure  34.  Axial  Velocity  Distribution  for  the  Four-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor, 
Measured  at  a  Radial  Distance  of  11.642  cm  from 
the  Hub  Surface 
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Figure  35.  Circumferential  Velocity  Distribution  for  the 
Four-Cycle  Distortion  with  the  Nine-Bladed 
Cambered  Rotor,  Measured  at  a  Radial  Distance 
of  11.642  cm  from  the  Hub  Surface 
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Figure  38.  Axial  Velocity  Distribution  for  the  Four-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor, 
Measured  at  a  Radial  Distance  of  2.328  cm  from 
the  Hub  Surface 
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the  same  run  conditions  as  tested  in  Figures  34  through  37,  but  these 
measurements  were  taken  at  the  radial  probe  location  nearest  to  the  hub 
surface.  Once  again,  the  same  distortion  shapes  and  phenomena  are 
observed  that  were  discussed  in  the  case  of  the  six-bladed  uncambered 
rotor  in  Section  4.2.1.  By  comparing  the  two  sets  of  graphs,  it  is 
noticed  that  the  magnitudes  of  the  circumferential-mean  axial  and  circum¬ 
ferential  velocities  and  total  and  static  pressures  are  similar  but  not 
identical  upstream  of  the  rotor  for  the  various  radial  locations  at  which 
the  measurements  were  taken.  The  slight  changes  observed  in  the  magnitudes 
are  due  to  the  radial  amplitude  variations  of  the  distortion  screens  as 
seen  in  Reference  [35].  The  probes  were  positioned  far  enough  away  from 
the  solid  boundaries  so  that  the  measurements  were  not  affected  by  the 
boundary  layers  from  the  casing  wall  and  the  hub  surface.  Therefore,  the 
changes  in  magnitude  of  the  velocity  and  pressure  distributions  upstream 
of  the  rotor  are  caused  by  the  design  of  the  screens,  and  they  exist  in 
the  absence  of  the  rotor  as  well. 

By  examining  Figures  36  and  40,  it  is  seen  that  at  the  different 
radial  locations,  there  is  a  minor  difference,  ten  percent,  in  the  magni¬ 
tude  of  the  circumferentially  averaged  total  pressure  downstream  of  the 
rotor.  It  is  also  noticed  that  the  fluctuations  in  total  pressure  are 
greater  at  the  tip  than  at  the  hub,  and  the  frequency  or  harmonic  content 
of  the  distortions  is  quite  different.  The  fact  that  the  magnitudes  of 
the  circumferentially  averaged  total  pressure  distributions  downstream  of 
the  rotor  are  different  at  the  two  extreme  spanwise  locations  is  not  too 
surprising  since  the  solidity  of  the  blades  is  different  at  these  two 
locations,  which  also  is  a  possible  explanation  for  the  observed  changes 
in  the  total  pressure  variations  at  the  two  radial  measurement  locations. 
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Figure  39  shows  an  increase  in  magnitude  of  the  local  values  of 
circumferential  velocity  downstream  of  the  rotor  near  the  hub  over  the 
circumferential  velocity  near  the  outer  casing  as  shown  in  Figure  35. 

This  merely  indicates  that  the  rotor  turns  the  flow  more  near  the  hub 

than  at  the  tip.  The  higher  turning  near  the  hub  is  a  result  of  the 

rotor  design  which  utilizes  a  free  vortex  spanwise  loading. 

Comparing  Figures  34  and  38,  it  is  seen  that  the  magnitude  of  the 
circumferential-mean  axial  velocity  distribution  is  greater  at  the  measure¬ 
ment  location  near  the  tip  than  near  the  hub  downstream  of  the  rotor, 
location  15.  This  is  due  to  the  hub  wall  boundary  layer  since  ^  is 

tne  circumferentially  averaged  axial  velocity  at  station  1-D.  Figures  34 

and  38  also  show  a  noticeable  difference  in  the  harmonic  content  of  the 
axial  velocity  distortion  downstream  of  the  rotor.  The  more  random  nature 
of  the  V  variation  in  Figure  38  could  be  due  to  greater  turbulent  mixing 

A 

near  the  hub  of  the  rotor. 

It  is  seen  from  Figures  37  and  41  that  the  magnitude  of  the  static 
pressure  downstream  of  the  rotor  is  larger  near  the  tip  than  near  the  hub 
of  the  rotor.  The  gradient  in  static  pressure  between  the  two  radial 
locations  is  due  simply  to  the  spanwise  loading  on  the  blades. 

A  second  series  of  distortion  measurements  was  made  with  the  nine- 
bladed  cambered  rotor  to  identify  the  effects  of  rotor-stator  spacing  on 
the  inlet  disturbance.  For  this  purpose,  eight  cambered  stator  blades 
were  installed  in  the  AFRF,  and  a  study  was  performed  employing  various 
distortion  types  and  blade  loadings.  The  measurements  in  this  series  of 
tests  were  all  taken  at  the  mean  radius  of  the  blades.  One  set  of  data 
was  obtained  with  the  leading  edge  of  the  stator  blade  row  located  two 
blade  chord  lengths  (30.48  cm)  behind  the  trailing  edge  of  the  rotor  blade 
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row.  As  in  the  tests  conducted  previously,  the  probes  were  located  at 
axial  stations  1-D  and  15.  Another  set  of  data  was  obtained  with  the 
stator  row  positioned  one-half  of  a  blade  chord  length  (7.62  cm)  down¬ 
stream  of  the  rotor  trailing  edge.  For  this  series  of  tests,  the  upstream 
probe  was  located  at  station  1-D,  but  the  downstream  probe  was  moved 
forward  to  a  new  axial  position  designated  as  location  15-A.  Location 
15-A  was  positioned  at  an  axial  distance  of  0.40  chord  length  downstream 
of  the  rotor  trailing  edge.  This  change  in  probe  location  was  necessary 
so  that  the  downstream  probe  remained  upstream  of  the  stator  blade  row. 

Figures  42  through  45  present  the  axial  velocity,  circumferential 
velocity,  total  pressure,  and  static  pressure  distributions  as  a  function 
of  circumferential  location  for  the  two-cycle  sinusoidal  distortion  with 
the  stator  row  located  two  blade  chord  lengths  behind  the  trailing  edge 
of  the  rotor.  In  contrast  to  these.  Figures  46  through  49  show  the  axial 
and  circumferential  velocities  and  the  total  and  static  pressures  measured 
using  the  same  operating  conditions  with  the  stator  blades  located  one- 
half  of  a  blade  chord  length  behind  the  rotor. 

It  is  quickly  observed  that  the  different  rotor-stator  spacings  had 
no  effect  whatsoever  on  the  data  measured  upstream  of  the  rotor,  location 
1-D.  The  energy  added  to  the  flow  was  the  same  in  both  cases  since  the 
measurements  were  all  made  at  the  mean  radius  and  the  same  operating 
conditions  were  tested.  Thus,  the  total  pressure  variation  is  unchanged 
downstream  of  the  rotor  as  seen  in  Figures  44  and  48.  The  most  important 
characteristic  of  the  flow  field  which  is  observed  by  changing  the  rotor- 
stator  spacing  is  the  shapes  of  the  axial  velocity,  circumferential 
velocity,  and  static  pressure  distributions  in  Figures  46,  47,  and  49, 
respectively.  A  periodic  fluctuation  in  these  distributions  occurs 
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Figure  42.  Axial  Velocity  Distribution  for  the  Two-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor 
for  a  Rotor-Stator  Spacing  of  Two  Blade  Chord 
Lengths  (30.48  cm) 
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Figure  43.  Circumferential  Velocity  Distribution  for  the 
Two-Cycle  Distortion  with  the  Nine-Bladed 
Cambered  Rotor  for  a  Rotor-Stator  Spacing  of 
Two  Blade  Chord  Lengths  (30.48  cm) 
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Figure  44.  Total  Pressure  Distribution  for  the  Two-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor 
for  a  Rotor-Stator  Spacing  of  Two  Blade  Chord 
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Figure  45.  Static  Pressure  Distribution  for  the  Two-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor 
for  a  Rotor-Stator  Spacing  of  Two  Blade  Chord 
Lengths  (30.48  cm) 
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Figure  46.  Axial  Velocity  Distribution  for  the  Two-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor 
for  a  Rotor-Stator  Spacing  of  0.5  Blade  Chord 
Lengths  (7.62  cm) 


Vs/ Vx  «v».  V#/Vx  m. 


9  BLADES  AVG.  FLON  COEF.  =  0-632 

50  DEG-  STAGGER  ANGLE  AVG.  P-RISE  COEF-  =  1-210 
2  CYCLE  DISTORTION  AVG.  INCIDENCE  =  1-56  DEG- 

RPM  =  1529 


Figure  47.  Circumferential  Velocity  Distribution  for  the 
Two-Cycle  Distortion  with  the  Nine-Bladed 
Cambered  Rotor  for  a  Rotor-Stator  Spacing  of 
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Figure  48.  Total  Pressure  Distribution  for  the  Two-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor 
for  a  Rotor-Stator  Spacing  of  0.5  Blade  Chord 
Lengths  (7 . 62  cm) 
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Figure  A9.  Static  Pressure  Distribution  for  the  Two-Cycle 
Distortion  with  the  Nine-Bladed  Cambered  Rotor 
for  a  Rotor-Stator  Spacing  of  0.5  Blade  Chord 
Lengths  (7.62  cm) 
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downstream  of  the  rotor.  It  is  also  important  to  note  that  the  fluctua¬ 
tion  goes  through  eight  cycles  in  a  complete  circumferential  survey  of 
360  degrees.  By  looking  at  Figures  46  and  49,  the  axial  velocity  and 
static  pressure  variations  downstream  of  the  rotor  are  seen  to  be  directly 
out  of  phase  which  is  expected  from  Bernoulli's  equation.  The  eight-cycle 
periodic  type  of  fluctuation  is  attributed  to  a  blockage  effect  due  to  the 
close  proximity  of  the  downstream  probe  to  the  stator  blades  when  the 
stator  blades  are  located  in  the  forward  position,  only  one-half  of  a 
chord  length  behind  the  rotor.  With  the  stator  row  in  the  forward  position, 
the  probe  was  placed  an  axial  distance  of  only  0.6  inch  (1.524  cm)  in 
front  of  the  stator  blades.  The  eight  troughs  in  the  axial  velocity 
fluctuation  (i.e.,  the  peaks  in  the  static  pressure  distribution)  represent 
the  points  at  which  the  probe  was  positioned  directly  in  front  of  the 
stator  blades.  Between  the  blades,  the  velocity  is  accelerated  and  this 
is  shown  by  the  peaks  in  the  axial  velocity  variation  at  location  15-A 
in  Figure  46. 

4.3  Distortion  Attenuation 

The  purpose  of  this  section  is  to  examine  some  of  the  overall  trends 
observed  in  the  distortion  data  as  a  function  of  the  type  of  distortion 
screen  and  rotor  operating  variables.  Because  of  the  large  amount  of  data 
associated  with  the  objective  of  providing  a  data  base  for  subsequent 
studies,  it  is  not  possible  to  examine  and  evaluate  all  of  the  distortion 
characteristics  in  this  paper.  For  this  reason,  only  the  attenuation  (or 
amplification)  of  the  axial  velocity  and  total  pressure  will  be  discussed. 

An  attempt  has  been  made  to  point  out  the  most  significant  parameters 
involved,  and  this  should  serve  as  an  aid  in  conducting  future  studies 


in  this  area. 
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Information  concerning  the  attenuation  of  an  inlet  distortion  through 
a  blade  row  is  important  when  designing  a  multi-stage  turbomachine,  since 
the  latter  stages  will  be  subjected  to  the  exit  flow  of  the  upstream 
stages.  It  may  be  desirable,  if  possible,  to  completely  eliminate  a  dis¬ 
tortion  in  the  first  few  stages  of  a  compressor  so  that  the  latter  stages 
can  operate  in  a  uniform  flow  with  improved  performance  and  stability. 

The  rate  and  amount  of  attenuation,  however,  may  be  limited  by  unsteady 
forces  or  losses  in  the  early  stages. 

The  fact  that  a  total  pressure  distortion  will  be  attenuated  by  a 
rotor  is  easily  understood  from  a  quasi-steady  viewpoint.  Higher  incidence 
angles  in  the  low  velocity  region  will  result  in  a  larger  total  pressure 
rise  (unless  separation  occurs)  than  in  the  remaining  regions  and  thus 
make  up  some  of  the  total  pressure  defect.  However,  the  effects  of  the 
distortion  frequency,  rotor  design,  and  operating  conditions  on  the  dis¬ 
tortion  attenuation  are  not  easily  assessed  theoretically  and  therefore, 
the  experiments  of  this  study  can  provide  some  of  this  information. 

In  order  to  quantify  the  attenuation  of  a  distortion,  a  means  of 
representing  the  amplitude  of  the  distortion  is  necessary.  One  method  of 
quantifying  the  size  of  a  distortion  is  to  represent  the  data  as  a  Fourier 
series  so  that  the  amplitudes  of  the  various  harmonics  can  be  examined  and 
compared.  For  this  reason,  the  dimensionless  axial  velocity  and  total 
pressure  data  were  Fourier  analyzed  by  computing  the  coefficients  of  a 
conventional  Fourier  series  by  using  a  numerical  integration  scheme. 

The  Fourier  analysis  routine  was  written  and  programmed  by  Yocum  for 
use  on  the  IBM  1130  Computing  System.  The  coefficients  of  the  series  were 
then  converted  to  an  amplitude  and  a  phase  angle  so  that  the  data  could 
be  represented  in  the  following  form: 
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The  use  of  Equation  (35)  for  representing  the  data  is  especially 
convenient  for  describing  the  sinusoidal  distortions,  where  one  harmonic 
will  be  dominant  and  can  be  used  to  represent  the  distortion  amplitude. 
Figure  50  contains  plots  of  the  amplitudes,  A^,  of  the  dimensionless  total 
pressure  and  axial  velocity  distortions  versus  the  harmonic  number,  n,  for 
a  four-cycle  distortion.  Three  sets  of  data  are  presented  on  each  graph; 
one  set  describes  the  distortion  without  the  rotor  installed,  while  the 
c:her  two  describe  the  distortions  at  locations  1-D  and  15  when  the 
uncambered  rotor  is  operating.  This  type  of  plot  is  useful  for  demon¬ 
strating  the  harmonic  content  of  the  distortion  and  for  identifying  the 
dominant  harmonics.  The  figure  also  shows  how  each  harmonic  is  affected 
by  the  presence  of  the  rotor.  As  expected  for  the  four-cycle  distortion. 
Figure  50  shows  that  the  fourth  harmonic  is  very  dominant,  and,  thus,  the 
amplitudes  of  the  fourth  harmonics  can  be  used  to  describe  the  distortion 
size.  Similarly,  the  fundamental  harmonic  of  the  one-,  two-,  and  six- 
cycle  sinusoidal  distortions  can  be  used  to  represent  their  distortion 
amplitudes. 

The  square-wave  distortions  were  also  represented  in  the  form  of 
Equation  (35)  so  that  the  harmonic  content  could  be  examined  and  the 
attenuation  of  the  individual  harmonics  compared  with  the  attenuation 
observed  for  the  sinusoidal  distortions.  Plots  similar  to  those  in 
Figure  50  are  presented  for  a  90  degree  square  distortion  in  Figure  51. 

In  this  case,  it  is  found  that  the  first  two  harmonics  are  significant. 
Although  both  of  these  harmonics  could  be  used  to  analyze  the  attenuation 
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Figure  50.  Dimensionless  Total  Pressure  and  Axial  Velocity 

Harmonic  Amplitudes  as  a  Function  of  the  Harmonic 
Number,  n,  for  a  Four-Cycle  Distortion 
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of  the  square  distortions,  no  single  parameter  from  the  Fourier  series 
can  be  used  to  properly  represent  the  distortion  size.  Therefore,  in 
addition  to  Fourier  analyzing  the  data,  a  single  parameter  representing 
the  amplitude  of  the  square  distortions  was  calculated  by  taking  the 
difference  between  the  maximum  and  minimum  values  of  the  axial  velocity 
and  total  pressure.  Thus,  the  amplitudes  of  the  square  inlet  distortions 
were  defined  as: 


and 


P  -  p 

T  max  T  min 

l/2pV* 

X  avg 
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X  max 


-  V. 


X  min 


X  avg 


(36) 
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These  definitions  may  appear  to  be  obvious  ones  since,  for  a  perfect 
square  distortion,  A  is  simply  the  difference  in  magnitude  between  the 
distorted  and  undistorted  sectors.  However,  it  is  necessary  to  define  A 
with  the  maximum  and  minimum  values  because  downstream  of  the  rotor  the 
remaining  distortion  is  no  longer  square.  This  fact  is  shown  in  Figure 
31  where  the  total  pressure  variation  for  a  square  inlet  distortion  is 
plotted  and  the  amplitude  as  defined  by  Equation  (36)  is  indicated. 

For  both  the  sinusoidal  and  the  square  distortions,  the  attenuation 
is  presented  as  a  ratio  of  the  distortion  amplitude  with  the  rotor 
operating  to  the  amplitude  of  the  distortion  at  the  same  axial  location 
without  the  rotor  operating.  With  the  data  presented  in  this  manner,  the 
amplitude  ratio  can  range  from  a  value  of  zero,  indicating  the  distortion 
was  completely  eliminated,  to  a  value  of  one,  indicating  no  attenuation. 
Values  greater  than  cne  mean  that  the  distortion  was  amplified. 
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4.3.1  Results  for  the  Slx-Bladed  Uncambered  Rotor.  Figure  52 
presents  the  attenuation  for  the  six-cycle  sinusoidal  distortion  in  the 
form  of  the  amplitude  ratio  of  the  sixth  harmonic  versus  the  circum¬ 
ferential  average  incidence  angle  for  the  various  stagger  angles  tested. 
Both  the  axial  velocity  and  total  pressure  amplitude  ratios  are  presented 
for  both  upstream  and  downstream  locations,  with  the  three  symbols  for 
each  location  indicating  the  different  blade  stagger  angles.  It  is  seen 
from  Figure  52  that  both  the  total  pressure  and  axial  velocity  attenuation 
upstream  of  the  rotor  are  not  strongly  influenced  by  a  change  in  incidence, 
or  mean  blade  loading.  Downstream  of  the  rotor,  the  V  and  P  distortions 

A  1 

show  increasing  attenuation  (i.e.,  a  decreasing  amplitude  ratio)  as  the 
incidence  angle  is  increased  for  the  rotor  configurations  having  all 
three  blade  stagger  angles.  In  Reference  [37],  incidence  angle  was  found 
to  have  very  little  effect  on  the  distortion  attenuation  with  the  twelve- 
bladed  uncambered  rotor.  Thus,  it  is  believed  that  an  influence  due  to 
the  space- to-chord  ratio  is  causing  the  observed  dependence  of  incidence 
on  the  attenuation  for  the  six-bladed  rotor.  In  the  present  study,  there 
is  not  a  sufficient  amount  of  data  to  prove  this  hypothesis. 

The  most  important  characteristic  of  the  data  revealed  in  Figure  52 
is  the  influence  of  the  blade  stagger  angle  on  the  amount  of  attenuation 
occurring  downstream  of  the  rotor.  The  rotor  configurations  with  the 
higher  blade  stagger  angles  produce  a  larger  attenuation  of  both  the 
total  pressure  and  axial  velocity.  The  increase  in  attenuation  of  the 
total  pressure  distortion  can  be  partially  attributed  to  the  steepening 
of  the  pressure  rise  versus  incidence  performance  curve  which  occurs  for 
the  higher  stagger  angles. 

With  the  exception  of  the  effect  of  incidence  on  the  distortion 
attenuation,  the  general  trends  observed  and  discussed  here  for  the 
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six-bladed  uncambered  rotor  were  also  observed  for  the  uncambered  rotor 
with  twelve  blades,  as  discussed  in  Reference  [37].  By  comparing  the  two 
sets  of  data,  it  is  found  that  the  magnitudes  of  both  the  total  pressure 
and  axial  velocity  amplitude  ratios  are  greater  for  the  six-bladed  rotor 
than  for  the  twelve-bladed  rotor.  Thus,  the  effect  of  decreasing  the 
rotor  space- to-chord  ratio  from  1.356  to  0.676  is  to  increase  the  amount 
of  attenuation  of  the  total  pressure  and  axial  velocity  inlet  distortions 
downstream  of  the  rotor. 

The  next  set  of  attenuation  data  to  be  presented  for  the  uncambered 
rotor  is  the  amplitude  ratios  of  the  total  pressure  and  axial  velocity 
distortions  as  a  function  of  incidence  angle  for  the  two  square-shaped 
inlet  distortions  tested.  These  data  are  presented  in  Figures  53  and  54 
for  the  90  degree  and  180  degree  distortions,  respectively.  The  total 
pressure  distortion  amplitude  is  defined  by  Equation  (36),  and  the  axial 
velocity  distortion  amplitude  is  defined  by  Equation  (37).  As  done 
previously  for  the  sinusoidal  distortion,  the  data  on  the  figures  are 
presented  with  the  symbols  representing  the  various  blade  stagger  angles. 

Figures  53  and  54  show  that  the  effect  of  incidence  or  blade  loading 
on  the  attenuation  of  the  axial  velocity  and  total  pressure  distortions 
is  minor.  The  distortions  undergoe  approximately  a  five  to  ten  percent 
change  in  attenuation  over  the  range  of  incidence  angles  tested.  It  is 
also  found  that  the  stagger  angle  does  not  significantly  affect  the 
attenuation  upstream  or  downstream  of  the  rotor  for  the  square  distortions. 

One  important  observation  to  be  made  in  Figure  53  for  the  90  degree 
square  distortion  is  the  change  in  the  total  pressure  distortion  amplitude 
upstream  of  the  rotor.  It  was  pointed  out  in  Section  4.2.1  that  the 
magnitude  of  the  total  pressure  distortion  upstream  of  the  rotor  had 
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Figure  53.  Measured  Axial  Velocity  and  Total  Pressure 

Distortion  Attenuation  versus  Incidence  Angle 
for  a  90  Degree  Square  Distortion  with  the 
Uncambered  Rotor 
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Distortion  Attenuation  versus  Incidence  Angle 
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increased  with  the  rotor  installed  in  the  case  of  a  90  degree  square 
distortion.  This  phenomenon  is  shown  clearly  in  Figure  53  where  the 
total  pressure  amplitude  ratios  upstream  of  the  rotor  are  greater  than 
one.  It  was  suggested  in  Section  4.2.1  that  the  increase  in  amplitude 
is  due  to  a  radial  total  pressure  gradient  produced  by  the  screen  followed 
by  a  radial  shift  in  the  streamlines  upstream  of  the  rotor.  However,  for 
the  sinusoidal  distortions,  the  variation  is  not  significantly  large  and 
the  mean  amplitude  ratio  is  approximately  1.0  (see  Figure  52). 

The  main  observation  which  can  be  made  from  Figures  53  and  54  is  the 
influence  of  the  circumferential  extent  or  the  width  of  the  distortion. 

It  is  seen  that  the  magnitudes  of  the  axial  velocity  amplitude  ratios 
upstream  and  downstream  of  the  rotor  are  quite  similar  for  both  the  90 
degree  and  180  degree  square  distortions.  The  total  pressure  amplitude 
ratios  downstream  of  the  rotor  also  have  approximately  the  same  magnitude. 
However,  it  appears  from  the  figures  that  the  90  degree  total  pressure 
distortion  has  undergone  a  greater  amplification  upstream  of  the  rotor 
than  the  180  degree  distortion.  The  attenuation  of  the  two  square  dis¬ 
tortions  across  the  rotor  is  approximately  the  same  since  the  axial 
velocity  and  total  pressure  amplitude  ratios  downstream  of  the  rotor  are 
the  same.  Due  to  the  hypothesis  that  three-dimensional  flows  are  involved 
in  amplifying  the  inlet  distortion,  it  cannot  be  concluded  from  this  data 
that  a  larger  amplification  of  the  total  pressure  distortion  upstream  of 
the  rotor  will  always  be  observed  for  smaller  width  distortions. 

From  the  data  taken  with  the  sinusoidal  distortions,  the  variation 
in  the  attenuation  was  found  to  be  much  greater  as  a  function  of  distortion 
frequency  than  the  variation  observed  with  changes  in  incidence  angle. 
Therefore,  the  remaining  data  with  the  uncambered  rotor  for  the  sinusoidal 
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distortions  will  be  presented  as  a  function  of  the  reduced  frequency, 
which  is  a  dimensionless  frequency  parameter  defined  as: 


2? r  £ 

4  '  2  * 


(38) 


This  definition  stems  from  unsteady  aerodynamic  theory  and  is  represented 
as  a  ratio  of  the  blade  chord  length  to  the  wavelength  of  the  disturbance 
along  the  chord.  The  typical  definition  of  reduced  frequency  is  based  on 
the  half-chord,  not  the  chord  length,  as  shown  in  Equation  (38).  The 
wavelength  of  the  disturbance  along  the  chord,  is  related  to  the  wave¬ 
length  along  the  pitch,  £,  in  the  following  manner: 


o  =  — - — 
c  sin  £  ’ 


(39) 


where  £  is  the  stagger  angle.  The  relation  in  Equation  (39)  assumes  that 
the  gust  or  disturbance  is  convected,  i.e. ,  it  moves  past  the  blade  at  the 
same  velocity  as  the  free-stream  velocity.  Substituting  Equation  (39) 
into  Equation  (38),  the  reduced  frequency  becomes: 


w  =  -j-  sin  £  .  (40) 

Equation  (40)  is  the  definition  of  to  which  is  used  throughout  this  study. 
As  a  physical  interpretation  of  w,  it  can  be  shown  that  the  reduced 
frequency  represents  the  number  of  distortion  cycles  across  the  blade 
chord. 

Figures  55  and  56  present  the  amplitude  ratio  of  the  total  pressure 


and  axial  velocity  distortions  versus  to  for  constant  values  of  incidence 
angle  of  zero  and  eight  degrees,  respectively.  It  is  seen  from  both 
Figures  55  and  56  that  an  increase  in  reduced  frequency  decreases  the 
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attenuation  of  the  total  pressure  and  axial  velocity  distortions  down¬ 
stream  of  the  rotor.  In  the  case  of  an  eight-degree  incidence  angle, 
there  is  a  slight  increase  in  the  attenuation  of  the  total  pressure  dis¬ 
tortion  upstream  of  the  rotor  configurations  having  blade  stagger  angles 
of  35  and  55  degrees.  The  decrease  in  the  distortion  attenuation  down¬ 
stream  of  the  rotor  with  increasing  values  of  w  is  contrary  to  the  results 
obtained  with  the  twelve-bladed  rotor  in  Reference  [37], 

It  is  suggested  that  this  decrease  in  attenuation  occurs  because  the 
reduced  frequency  is  changing  to  give  the  situation  where  the  number  of 
blades  on  the  rotor  and  the  number  of  distortion  cycles  are  equal,  i.e., 
the  six-bladed  rotor  and  six-cycle  distortion.  In  this  case,  the  wave¬ 
length  of  the  disturbance  along  the  pitch  direction  is  equal  to  the  blade 
spacing,  and  each  blade  will  experience  the  same  disturbance  velocity  at 
a  given  time.  It  is  hypothesized  that  if  higher  values  of  reduced 
frequency  (for  instance,  nine-cycle  and  twelve-cycle  sinusoidal  distor¬ 
tions)  had  been  tested  in  this  investigation,  an  increase  in  attenuation 
would  be  observed  at  the  higher  values  of  w.  It  is  further  hypothesized 
that  if  a  twelve-cycle  sinusoidal  distortion  had  been  tested  in  the  study 
in  Reference  [37],  a  decrease  in  attenuation  would  have  been  observed  at 
the  corresponding  value  of  reduced  frequency.  These  hypotheses  seem  to 
be  supported  by  the  attenuation  data  for  the  nine-bladed  cambered  rotor 
which  will  be  presented  in  the  following  section,  but  further  experimen¬ 
tation  with  the  six-bladed  rotor  would  be  required  to  completely  verify 


the  assumption. 
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4-3.2  Results  for  the  Nine-Bladed  Cambered  Rotor.  The  next 
series  of  attenuation  data  to  be  presented  shows  the  results  obtained 
with  the  nine-bladed  cambered  rotor.  These  data  are  presented  in  a 
manner  similar  to  the  data  for  the  uncambered  rotor.  The  effect  of 
stagger  angle  can  no  longer  be  considered  since  the  cambered  rotor  has  a 
fixed  stagger  angle  of  50  degrees  at  its  mean  radius.  The  effects  of 
rotor-stator  spacing  and  radial  probe  location  were  studied  with  the 
cambered  rotor;  therefore,  the  influence  of  these  variables  on  the  atten¬ 
uation  of  the  axial  velocity  and  total  pressure  distortions  will  now  be 
discussed. 

Figure  57  shows  the  attenuation  for  the  four-cycle  sinusoidal  dis¬ 
tortion  in  the  form  of  the  amplitude  ratio  of  the  fourth  harmonic  versus 
the  circumferential-mean  incidence  angle  for  the  two  rotor-stator  spacings 
investigated.  Both  the  axial  velocity  and  total  pressure  amplitude  ratios 
are  presented  for  both  locations,  with  the  circles  indicating  a  rotor- 
stator  spacing  of  two  blade  chord  lengths  (30.48  cm)  and  the  squares 
denoting  a  rotor-stator  spacing  of  0.5  blade  chord  length  (7.62  cm).  All 
measurements  were  taken  at  the  mean  radius  of  the  blade.  The  first  thing 
to  be  noticed  in  Figure  57  is  that  rotor-stator  spacing  has  no  effect  on 
the  total  pressure  and  axial  velocity  distortions  downstream  of  the  rotor. 
This  indicates  that  the  inlet  distortion  is  attenuated  within  at  least 
0.4  chord  length  of  the  rotor.  Once  again,  it  is  seen  that  the  total 
pressure  distortion  is  slightly  amplified  upstream  of  the  rotor,  and  this 
amplification  is  very  small.  Another  important  observation  to  be  made 
from  Figure  57  is  the  large  amount  of  attenuation  experienced  between  the 
upstream  and  downstream  locations  for  both  the  axial  velocity  and  total 
pressure  distortions.  The  cambered  rotor  is  seen  to  almost  completely 
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eliminate  the  four-cycle  sinusoidal  inlet  distortion.  It  can  be  pointed 
out  that  this  large  amount  of  attenuation  downstream  of  the  rotor  was 
also  observed  for  the  six-cycle  distortion;  however,  it  was  not  found  to 
occur  for  the  one-  and  two-cycle  sinusoidal  distortions.  In  all  cases, 
however,  it  was  commonly  observed  that  both  the  incidence  angle  and  the 
rotor-stator  spacing  had  a  negligible  effect  on  the  distortion  attenua¬ 
tion.  The  negligible  effect  of  incidence  or  blade  loading  on  the  attenua¬ 
tion  is  attributed  to  the  linear  nature  of  the  performance  curve  corre¬ 
sponding  to  the  range  of  operating  conditions  used  in  this  series  of  tests. 

The  next  set  of  attenuation  data  for  the  cambered  rotor  shows  the 
total  pressure  and  axial  velocity  amplitude  ratios  as  a  function  of 
incidence  angle  for  the  two  square-shaped  inlet  distortions  tested. 

Figures  58  and  59  present  the  results  for  the  90  degree  and  180  degree 
distortions,  respectively.  It  can  be  seen  that  the  incidence  angle  again 
does  not  significantly  influence  the  attenuation  of  the  axial  velocity 
and  total  pressure  distortions.  In  the  case  of  the  square  distortions, 
it  is  observed  that  the  effect  of  rotor-stator  spacing  on  the  distortion 
attenuation  downstream  of  the  rotor  is  no  longer  insignificant.  Figure 
58  shows  approximately  a  20  percent  difference  between  the  attenuation 
of  the  axial  velocity  distortion  produced  with  a  rotor-stator  spacing  of 
2.0  chord  lengths  and  that  produced  with  a  rotor-stator  spacing  of  0.5 
chord  length.  A  similar  difference  in  the  total  pressure  distortion 
downstream  of  the  rotor  is  shown  in  Figure  59  for  the  180  degree  distor¬ 
tion.  The  reason  for  this  effect  of  rotor-stator  spacing  is  not  known. 

The  observation  was  made  for  the  uncambered  rotor  that  the  total  pressure 
distortion  upstream  of  the  rotor  had  been  amplified  more  for  the  90  degree 
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distortion  than  for  the  180  degree  distortion.  This  is  also  true  for  the 
case  of  the  cambered  rotor  as  can  be  seen  in  Figures  58  and  59 

One  final  important  characteristic  of  the  data  in  Figures  58  and  59 
concerns  the  amplitude  ratios  upstream  and  downstream  of  the  rotor.  The 
total  pressure  amplitude  ratio  for  both  the  90  degree  and  180  degree  dis¬ 
tortions  show  that  the  distortion  is  amplified  upstream  of  the  rotor  and 
then  undergoes  a  considerable  amount  of  attenuation  downstream  of  the 
rotor.  However,  it  is  seen  from  the  figures  that  the  attenuation  of  the 
axial  velocity  distortion  for  a  rotor-stator  spacing  of  0.5  chord  length 
actually  decreases  as  the  distortion  passes  through  the  rotor.  As  dis¬ 
cussed  in  Section  4.2.2,  periodic  fluctuations  in  the  axial  velocity  dis¬ 
tribution  appeared  downstream  of  the  rotor  due  to  a  blockage  effect  when 
the  stator  row  was  in  its  forward  position.  The  amplitude  of  the  square 
distortions  is  defined  as  the  difference  between  the  maximum  and  minimum 
values  of  the  distortion.  Since  the  axial  velocity  distribution  upstream 
of  the  rotor  did  not  experience  high  peaks  and  low  valleys  due  to  a 
blockage  effect,  it  is  concluded  that  the  amplitude  of  the  axial  velocity 
distortion  is  greater  downstream  than  upstream  of  the  rotor.  This 
phenomenon  was  not  observed  with  the  sinusoidal  distortions  because  the 
amplitudes  were  defined  in  terms  of  the  fundamental  harmonic  of  the  Fourier 
series. 

In  a  manner  similar  to  the  case  of  the  uncambered  rotor,  the  attenu¬ 
ation  data  for  the  cambered  rotor  was  also  plotted  as  a  function  of  reduced 
frequency.  Figures  60  and  61  present  the  amplitude  ratio  of  the  total 
pressure  and  axial  velocity  distortions  versus  w  for  two  constant  values 
of  incidence,  i  =  0  degrees  and  i  =  5  degrees,  respectively.  The  values 
of  reduced  frequency  correspond  to  the  one-,  two-,  four-,  and  six-cycle 
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Figure  60.  Amplitude  Ratio  versus  Reduced  Frequency  for  a 
Constant  Mean  Incidence  Angle  of  Zero  Degrees 
with  the  Cambered  Rotor 
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sinusoidal  distortions.  By  comparing  the  magnitudes  of  the  amplitude 
ratios  in  the  two  figures,  it  is  easily  seen  again  that  incidence  angle 
has  little  effect  on  the  distortion  attenuation.  Also,  the  effect  of 
rotor-stator  spacing  is  again  found  to  be  insignificant.  Downstream  of 
the  rotor,  the  attenuation  of  the  total  pressure  and  axial  velocity  dis¬ 
tortions  is  increased  with  increasing  values  of  reduced  frequency. 

The  final  set  of  attenuation  data  for  the  cambered  rotor  presents  the 
response  of  the  rotor  on  the  distortion  where  the  measurements  were  made 
at  radial  locations  other  than  the  mean  radius  of  the  blade.  Two  radial 
or  spanwise  locations  were  less  than  the  mean  radius  and  two  were  greater 
than  the  mean  radius.  Figures  62  and  63  present  the  amplitude  ratio  of 
the  total  pressure  and  axial  velocity  distortions  as  a  function  of  the 
radial  distance  of  the  five-hole  probe  from  the  surface  of  the  hub  for  the 
sinusoidal  and  square  distortions,  respectively.  Figure  64  shows  the 
attenuation  data  versus  reduced  frequency  for  a  constant  mean  incidence 
angle  of  five  degrees.  The  four  curves  on  the  graph  represent  the  four 
radial  measurement  locations  that  were  tested.  Figure  64  shows  the 
expected  result  of  increased  attenuation  downstream  of  the  rotor  with 
increasing  reduced  frequency.  The  main  observation  to  be  made  from  all 
three  figures  is  that  the  effect  of  different  radial  or  spanwise  measure¬ 
ment  locations  on  the  attenuation  of  the  total  pressure  and  axial  velocity 
distortions  is  small.  In  general,  the  attenuation  of  the  total  pressure 
and  axial  velocity  distortions  downstream  of  the  rotor  for  both  the 
sinusoidal  and  square-shaped  distortions  increases  closer  to  the  hub.  It 
appears,  however,  that  the  increase  is  only  approximately  ten  to  fifteen 
percent  between  the  two  extreme  radial  locations.  It  was  pointed  out  in 
Section  4.2.2  that  the  total  pressure  flow  field  was  only  slightly  changed 
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Figure  62.  Amplitude  Ratio  versus  Radial  Location  for  a 
Constant  Mean  Incidence  Angle  of  Five  Degrees 
with  the  Cambered  Rotor 
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Figure  63.  Amplitude  Ratio  versus  Radial  Location  for  a 
Constant  Mean  Incidence  Angle  of  Five  Degrees 
with  the  Cambered  Rotor 
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and  the  axial  velocity  flow  field  experienced  a  larger  change  in  magnitude 
at  different  radial  locations  along  the  blade.  These  observations  appear 
to  be  supported  by  the  attenuation  data  shown  in  Figures  62,  63,  and  64. 


4.4  Unsteady  Total  Pressure  Losses 

The  operation  of  a  rotating  blade  row  in  either  a  uniform  or  a  dis¬ 
torted  flow  results  in  losses  in  total  pressure.  As  a  consequence  of  these 
losses,  the  total  pressure  rise  predicted  by  ideal  fluid  theories  is  not 
realized.  In  an  actual  turbomachine,  these  losses  are  the  result  of  skin 
friction  and  pressure  drag  on  the  blades,  viscous  drag  on  the  casing  walls, 
and  secondary  flows  due  to  three-dimensional  effects.  While  the  contri¬ 
bution  of  each  of  these  losses  is  well  documented  for  a  rotor  in  a  steady, 
uniform  flow,  there  is  little,  if  any,  data  available  regarding  the  losses 
in  a  distorted  inflow. 

In  Section  4.1,  the  total  pressure  loss  coefficient  for  steady  flow 

was  defined  by  Equation  (34).  The  same  definition  is  used  to  represent 

the  unsteady  losses  in  this  section,  except  that  the  circumferential 

averaged  values  of  P„  and  VQ  are  used.  Thus,  the  unsteady  actual  and 

i  y 

ideal  and  total  pressure  rise  and  the  total  pressure  loss  coefficients 
are  defined  by  the  following  equations: 


V  = 


^T  exit  ^T  inlet 
l/2pU2 


(41) 


f  = 


2  (V  -  V  ) 

9  exit  0  inlet2 


(42) 


and 


CpT  loss  '  W  '  ' 


(43) 


where  the  bar  over  P_  and  VQ  indicates  circumferential-mean  values.  A 

1  U 

ratio  of  the  unsteady  loss  coefficient  to  the  steady  loss  coefficient  at 
the  same  mean  incidence  angle  can  then  be  formed. 

The  studies  which  have  been  conducted  in  this  area  have  been  con¬ 
cerned  with  the  effects  of  upstream  unsteadiness  in  turbine  cascades 
[40,  41].  These  studies  were  conducted  with  a  stationary  turbine  cascade 
and  a  set  of  rotating  cylinders  located  upstream  of  the  cascade.  As  a 
result,  the  cascade  experienced  a  time  varying  inflow.  The  conclusions 
from  these  studies  are  that  the  presence  of  an  unsteady  flow  increases 
the  total  pressure  loss  as  compared  with  that  experienced  in  a  uniform 
flow.  In  both  investigations,  the  loss  associated  with  the  unsteadiness 
increased  with  the  speed  of  rotation  of  the  cylinders,  i.e.,  reduced 
frequency  and,  in  some  cases,  were  twice  the  magnitude  of  the  steady 
loss . 

A  portion  of  the  data  obtained  in  this  study  can  be  analyzed  to 
determine  the  total  pressure  losses  associated  with  a  distorted  inflow. 
Since  most  of  the  flow  surveys  were  conducted  at  the  mean  blade  radius, 
the  losses  experienced  are  primarily  those  due  to  skin  friction  and 
pressure  drag  on  the  blades.  The  losses  due  to  the  viscous  drag  on  the 
casing  walls  are  not  present  and  the  effects  of  secondary  losses  are 
minimal  for  the  studies  conducted  at  the  mean  radius  [42].  Considering 
only  the  influence  of  skin  friction  and  pressure  drag  to  contribute 
to  these  losses,  it  is  concluded  on  a  quasi-steady  basis  that  the  loss 
in  the  distorted  inflow  will  be  greater  than  that  experienced  in  a 
uniform  flow  since  the  blades  are  locally  experiencing  higher  angles  of 
incidence.  This  follows  from  data  of  cascades  in  a  uniform  or  steady  flow 
showning  an  increase  in  cascade  loss  with  increasing  angle  of  incidence. 
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Furthermore,  one  would  suspect  the  presence  of  an  unsteady  loss  which  is 
a  function  of  the  type  of  distortion  in  which  the  blade  row  is  operated, 
i.e.,  reduced  frequency. 

In  a  manner  similar  to  that  described  in  Section  4.1  for  uniform 
flow,  the  ideal  and  actual  circumferential-mean  total  pressure  rise 
coefficients  can  be  obtained  for  the  blade  row  operated  in  a  distorted 
inflow.  The  ideal  circumferential-mean  total  pressure  rise  for  a  sinus¬ 
oidal  distorted  inflow  is  identical  to  that  in  a  uniform  flow  if  the 
pressure  rise  characteristic  is  linear.  Higher  blade  stagger  angles 
result  in  a  larger  component  of  circumferential  velocity  downstream  of  the 
bxade  row  since  the  amount  of  swirl  or  turning  of  the  flow  is  increased 
as  t,  is  increased.  However,  to  maintain  equal  values  of  incidence  or 
blade  loading  at  different  stagger  angles,  the  blade  speed,  U,  must  also 
change  (for  a  constant  axial  velocity).  Therefore,  within  experimental 
accuracy,  the  ideal  total  pressure  rise  in  uniform  flow  does  not  change 
for  the  various  blade  stagger  angles.  Figure  65  presents  some  typical 
data  for  the  six-bladed  uncambered  rotor  operated  in  a  four-  and  six- 
cycle  sinusoidal  distortion  flow  field.  From  Equation  (34),  the  difference 
between  the  actual  and  ideal  circumferential-mean  total  pressure  rise 
coefficient  represents  the  total  pressure  loss  associated  with  operation 
in  the  distorted  inlet  flow. 

The  ratio  of  the  unsteady  total  pressure  loss  coefficient  to  the 
steady  total  pressure  loss  coefficient  calculated  from  the  data  for  the 
sinusoidal  distortions  is  presented  as  a  function  of  reduced  frequency 
for  the  uncambered  and  cambered  test  rotors  in  Figures  66  and  67, 
respectively.  The  reduced  frequency  was  previously  defined  in  equation 
(40);  in  a  uniform  flow,  the  distortion  wavelength  is  infinite  and  hence 
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Figure  66.  Ratio  of  Unsteady  Loss  to  Steady  Loss  versus 
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Figure  67.  Ratio  of  Unsteady  Loss  to  Steady  Loss  versus 

Reduced  Frequency  for  the  Nine  Bladed  Cambered 
Rotor 
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The  studies  of  Yocum  [37]  and  Kirillov  et  al.  [40]  reveal  that  the 
unsteady  losses  are  highly  dependent  on  reduced  frequency.  Their  data 
show  that  the  unsteady  losses  reach  a  maximum  at  low  values  of  reduced 
frequency  and,  in  some  cases,  the  peak  values  of  unsteady  losses  are 
greater  than  three  times  the  losses  in  steady  flow.  No  such  peak  is 
observed  in  the  data  of  this  study  as  shown  in  Figures  66  and  67.  For 
the  case  of  the  uncambered  rotor  in  Figure  66,  no  experimental  data  was 
obtained  to  determine  whether  the  unsteady  losses  reach  a  maximum  at  low 
values  of  reduced  frequency.  It  is  necessary  to  perform  experiments  at 
low  reduced  frequencies,  i.e.,  in  a  distortion  flow  field  with  one-  and 
tw,  — cycle  sinusoidal  variations  of  total  pressure,  in  order  to  verify  the 
conclusions  reached  in  References  [37]  and  [40]. 

Based  on  the  data  available,  an  attempt  was  made  to  correlate  the 
total  pressure  loss  data  obtained  in  this  study  for  the  six-bladed 
uncambered  rotor  and  the  data  obtained  in  Reference  [37]  for  the  twelve- 
bladed  uncambered  rotor.  This  was  done  to  try  to  prove  the  existence  of 
peak  unsteady  losses  at  low  distortion  frequencies.  To  do  this,  the 
ratio  of  unsteady  loss  to  steady  loss  was  plotted  as  a  function  of  the 
ratio  of  rotor  blade  spacing  to  distortion  wavelength  (the  reason  for 
using  this  ratio  will  be  discussed  in  Chapter  V),  and  the  results  are 
shown  in  Figure  68. 

Part  of  the  change  of  unsteady  losses  with  respect  to  frequency 
observed  by  Kirillov  [40]  was  due  to  a  change  in  the  flow  nonuniformity 
as  the  frequency  changed.  Additional  experimental  and  analytical  investi¬ 
gations  are  required  to  identify  the  causes  of  the  unsteady  losses  and 
to  explain  the  dependency  on  frequency. 

Figure  69  presents  the  ratio  of  unsteady  to  steady  total  pressure 
loss  as  a  function  of  reduced  frequency  for  the  four  different  radial 
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Figure  69,  Ratio  of  Unsteady  Loss  to  Steady  Loss  versus 

Reduced  Frequency  at  the  Four  Radial  Locations 
for  the  Cambered  Rotor 
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locations  at  which  the  circumferential  flow  surveys  were  taken  using  the 
nine-b laded  cambered  rotor.  It  can  be  seen  that  the  unsteady  losses  are 
much  greater  near  the  hub  of  the  blade  than  near  the  tip.  In  fact,  the 
data  show  that  the  loss  ratio  is  close  to  zero  near  the  tip  for  the  values 
of  reduced  frequency  that  were  tested.  It  is  possible  that  local  flow 
separation  near  the  tip  could  cause  large  values  of  steady  total  pressure 
loss. 

From  an  unsteady  viewpoint,  additional  losses  in  total  pressure  may 
be  caused  by  the  vortex  shedding  associated  with  the  unsteady  circulation 
on  the  blades.  The  unsteady  pressure  distribution  on  the  blades  will  also 
affect  the  growth  and  stability  of  the  boundary  layer,  and  the  effect  of 
the  unsteadiness  on  the  boundary  layer  is  expected  to  be  dependent  on  the 
reduced  frequency.  The  data  in  Figure  69  indicate  that  the  losses  due  to 
the  boundary  layer  at  the  hub  surface  are  much  higher  than  the  losses  due 
to  the  annulus  wall  boundary  layer.  Also,  the  losses  due  to  secondary 
flows  near  the  hub  seem  to  be  important  and  may  account  for  the  increase 
in  the  loss  curves.  Furthermore,  it  is  suspected  that  a  radial  shift  in 
streamlines  could  seriously  affect  the  unsteady  losses.  More  detailed 
and  extensive  radial  and  circumferential  flow  surveys  are  required  in  a 
flow  where  a  radial  shift  in  the  streamline  occurs,  and  the  data  should 
then  be  energy  or  mass  averaged.  Although  some  data  at  different  radial 
locations  have  been  presented  here,  future  studies  with  regard  to  a 
proper  averaging  technique  should  be  conducted  to  identify  the  causes  of 
the  unsteady  losses.  The  contributions  to  the  unsteady  losses  from  various 
sources  such  as  hub  and  annulus  wall  boundary  layer,  secondary  flows, 
vortex  shedding,  and  viscous  drag  should  also  be  examined  in  more  detail. 
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Radial  surveys  in  a  uniform  flow  could  also  be  conducted  to  determine  the 
effects  of  the  hub  and  annulus  wall  boundary  layers  on  the  steady  total 
pressure  losses. 


CHAPTER  V 


COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 


In  this  section,  the  results  of  the  theoretical  analysis  discussed 
in  Chapter  II  are  presented  and  compared  with  experimental  results 
obtained  in  the  AFRF.  The  results  for  the  six-bladed  and  twelve-bladed 
uncambered  rotors  and  the  nine-bladed  cambered  rotor  with  space-to-chord 
ratios  of  1-353,  0.676,  and  0.90,  respectively,  are  presented.  The 
data  for  the  twelve-bladed  rotor  was  provided  by  Yocum  [37].  These 
data  are  then  compared  with  the  predicted  results  of  the  time-mean 
total  pressure  downstream  of  the  rotor  obtained  by  the  use  of  Equation 
(19). 

In  addition  to  the  reduced  frequency,  the  unsteady  response  of  a 
cascade  or  blade  row  is  also  dependent  upon  a  pitchwise  frequency 
parameter,  the  intrablade  frequency  (or  intrablade  phase  angle)  which  is 
defined  as 


T  = 


2ttS 
£  » 


(44) 


where  S  is  the  rotor  blade  spacing  and  £  is  the  wavelength  of  the  dis¬ 
tortion  in  the  pitchwise  direction.  This  frequency  parameter  relates  the 
unsteady  response  of  the  neighboring  blades  of  the  blade  row  or  cascade 
to  that  of  the  reference  blade  at  a  given  instant  in  time.  For  a  rigid, 
non-vibrating  cascade  of  blades,  the  intrablade  frequency  is  related  to 
the  reduced  frequency,  u),  as  follows: 

T  =  -  2w  |  CSC£  .  (45) 


The  experimental  and  theoretical  results  presented  in  this  chapter 
are  plotted  as  the  ratio  of  the  distortion  amplitude  with  the  rotor 
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operating  to  the  distortion  amplitude  at  the  rotor  inlet  without  the 
rotor  operating  as  a  function  of  the  ratio  of  blade  spacing  to  distortion 
circumferential  wavelength.  The  blade  spacing  and  distortion  circum¬ 
ferential  wavelength  are  defined,  respectively,  as: 


and 


S 


2ttt 

B 


SL 


2rrr 

N 


(46) 


(47) 


where  r  is  the  radius  at  which  the  measurements  were  taken  (in  most  cases, 
the  mean  radius),  B  is  the  number  of  rotor  blades,  and  N  is  the  number  of 
distortion  cycles.  Combining  Equations  (46)  and  (47),  the  ratio  of  blade 
spacing  to  distortion  wavelength  can  be  written  as: 


S  _  2trr/B  _  N 
Z  ~  2nr/N  B  ' 


(48) 


Figures  70  and  71  present  the  experimental  values  of  the  distortion 
amplitude  ratio  versus  the  ratio  s /£.  Both  the  axial  velocity  and  total 
pressure  distortions  upstream  and  downstream  of  the  rotor  are  plotted  in 
the  same  figure.  Figure  70  shows  the  experimental  data  for  the  six-bladed 
and  twelve-bladed  rotors  for  blade  stagger  angles  of  35  and  45  degrees. 
Figure  71  shows  the  experimental  data  for  the  six-bladed  and  twelve-bladed 
rotors  for  stagger  angles  of  55  degrees  and  also  the  data  for  the  nine- 
bladed  rotor  with  a  stagger  angle  of  50  degrees.  It  was  determined  in 
Section  4.3.1  that  stagger  angle  is  an  important  parameter  which  influences 
the  attenuation  of  a  distortion.  Therefore,  the  data  on  Figures  70  and  71 
are  grouped  together  according  to  nearly  equivalent  stagger  angles. 
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Figure  71.  Measured  Values  of  Amplitude  Ratio  versus  Ratio 
of  Blade  Spacing  to  Distortion  Wavelength  for 
Stagger  Angles  of  50  and  55  Degrees 
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The  most  important  observation  to  be  made  in  Figures  70  and  71  is 
the  effect  on  the  distortion  amplitude  which  produces  a  maximum  attenua¬ 
tion  when  the  blade  spacing  is  equal  to  one-half  the  wavelength  of  the 
distortion.  This  maximum  attenuation  downstream  of  the  rotor  corresponds 
to  a  value  of  intrablade  phase  angle  equal  to  -u  and  different  values  of 
u>  dependent  upon  the  value  of  S/C  and  £.  As  seen  from  the  figures,  the 

attenuation  then  decreases  as  S/l  approaches  one,  i.e.,  as  the  intrablade 

phase  angle  approaches  -2tt.  It  should  be  noted  here  that  the  negative 
signs  on  tt  and  2ir  come  from  the  definition  of  X  in  Equation  (44)  which  is 
is  the  result  of  the  coordinate  system  which  is  chosen  for  the  cascade 

[21].  In  the  present  system,  the  cascade  is  moving  to  the  right,  and 

the  adjacent  blades  to  the  right  of  the  reference  blade  are  designated 
as  positive.  Thus,  the  blades  to  the  left  of  the  reference  blade,  which 
are  designated  as  negative,  see  the  same  disturbance  as  the  reference 
blade  at  time  t>o.  Therefore,  a  negative  sign  occurs  simply  because 
of  the  way  in  which  the  blades  in  the  cascade  (or  blade  row)  were 
counted. 

The  conclusion  that  the  ratio  of  blade  spacing  to  distortion  wave¬ 
length  is  a  primary  factor  in  the  attenuation  of  the  inlet  distortion  and 
that  a  value  of  S/l  =  0.5  (i.e.,  x  =  -it)  produces  a  maximum  attenuation 
is  very  significant.  This  is  a  condition  when  the  unsteady  response  of 
adjacent  blades  in  the  blade  row  are  180  degrees  out  of  phase;  at  a  given 
instant  of  time,  if  a  blade  experiences  a  maximum  distortion  velocity  in 
a  sinusoidal  distortion,  the  adjacent  blades  are  experiencing  a  minimum 
distortion  velocity.  Correspondingly,  at  S/l  **  1.0  (t  =  -2tt),  all  of  the 
blades  experience  a  similar  portion  of  the  sinusoidal  distortion  at  a 
given  instant  of  time  and  the  amount  of  attenuation  is  minimized.  Thus, 
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the  unsteady  cascade  effect  on  the  distortion  is  very  significant  and,  if 
the  blade  spacing  and  distortion  wavelength  are  properly  matched  to  give 
S/£  =0.5,  a  maximum  attenuation  of  the  distortion  amplitude  will  occur. 

The  data  for  all  three  rotors  shown  in  Figure  71  appear  to  collapse 
onto  the  same  curve  downstream  of  the  rotor.  It  is  also  seen  that  the 
effect  of  stagger  angle  is  consistent  between  the  six-  and  twelve-bladed 
rotors.  Downstream  of  the  rotor,  however,  the  data  for  the  six-bladed 
rotor  at  zero  degree  incidence  does  not  collapse  on  the  curve;  these  data 
show  a  significant  decrease  in  the  level  of  attenuation.  The  most  probable 
cause  for  this  is  an  effect  due  to  the  space-to-chord  ratio.  For  a  higher 
space-to-chord  ratio,  there  is  less  variation  in  flow  angle  between  the 
blades.  This  could  result  in  less  turbulent  mixing  and,  therefore,  less 
attenuation  downstream  of  the  rotor.  However,  to  verify  this  effect 
requires  data  showing  the  distortion  amplitude  ratio  as  a  function  of 
S/C  for  constant  values  of  S/Z.  At  present,  such  data  do  not  exist, 
although  the  two  points  at  S/Z  =  0.666  in  Figure  71  exhibit  the 
anticipated  trend. 

Figures  72  and  73  show  the  results  of  the  total  pressure  amplitude 
downstream  of  the  rotor  predicted  by  the  theoretical  analysis  described 
in  Chapter  II.  Also  shown  in  these  figures  are  the  experimental  data 
measured  downstream  of  the  rotor.  An  important  feature  of  these  two 
graphs  is  that  the  amount  of  attenuation  measured  experimentally  is  found 
to  be  much  greater  than  the  attenuation  predicted  by  theory.  In  fact, 
the  theory  predicts  an  amplification  of  the  distortion.  This  is  due  to 
the  magnitude  of  the  unsteady  blade  circulation  which  increases  the  dis¬ 
tortion  amplitude  above  the  amplitude  experienced  in  the  absence  of  the 
rotor.  The  magnitude  of  unsteady  blade  circulation  plays  a  very 
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significant  role  in  determining  the  amplitude  of  the  distortion  in  the 
theoretical  analysis.  As  the  magnitude  of  unsteady  blade  circulation 
increases,  the  amplitude  ratio  of  the  total  pressure  distortion  also 
increases  since  the  theory  does  not  account  for  any  viscous  losses.  An 
increase  in  the  amplitude  ratio  is  equivalent  to  a  decrease  in  the 
attenuation  of  the  distortion.  The  large  difference  between  the  amount 
of  attenuation  predicted  by  theory  and  that  measured  experimentally  is 
attributed  to  blade  thickness,  boundary  layer,  and  turbulence  mixing 
effects  which  are  assumed  negligible  in  the  theory.  There  is  also  the 
possibility  of  unsteady  radial  flows  which  result  from  a  radial  gradient 
of  unsteady  circulation  and  are  not  included  in  the  two-dimensional 
theoretical  analysis. 

It  was  stated  previously  that  the  experimental  results  show  a  maximum 
attenuation  at  S /l  =0.5  and  then  a  decrease  in  attenuation  as  the  ratio 
of  spacing  to  wavelength  approaches  on?.  From  Figures  72  and  73,  it  is 
seen  that  the  predicted  results:  also  demonstrate  an  increasing  attenuation 
as  a  value  of  S/£.  =  0.5  is  approached.  However,  as  S/£  is  increased  above 
0.5  for  the  six-  and  twelve-bladed  rotors,  the  predicted  attenuation  also 
shows  a  slight  increase.  The  data  obtained  in  this  study  for  the  nine- 
bladed  cambered  rotor  employed  the  one-,  two-,  four-,  and  six-cycle 
sinusoidal  distortion  screens.  These  distortions  correspond  to  values  of 
S/l  equal  to  0.111,  0.222,  0.444,  and  0.667,  respectively.  For  these 
points,  as  seen  on  Figure  73,  the  trend  of  the  predicted  results  agree 
quite  well  with  the  trend  of  the  measured  values.  In  an  attempt  to  verify 
t  '•!*•  increase  in  the  experimental  attenuation  which  is  observed  for  the 
c  urbered  rotor  as  S/l  approaches  one,  a  set  of  AFRF  data  was  obtained 

■  r. i. e  for  the  nine-bladed  rotor  in  a  nine-cycle  sinusoidal 


distortion  flow  field.  This  yields  a  ratio  of  blade  spacing  to  distortion 
wavelength  of  one.  From  this  set  of  data,  as  shown  in  Figure  73,  it  was 
again  found  that  the  predicted  attenuation  increased  as  S/l  goes  to  one. 

It  should  be  noted  that  the  points  shown  on  the  theoretical  curves  in 
Figures  72  and  73  are  the  points  at  which  the  theory  was  applied  and  are 
based  on  experimental  measurements  of  the  distortion  at  the  rotor  inlet. 
The  shape  of  the  curves  between  these  discrete  points  is  not  defined. 
Additional  data  utilizing  more  distortion  screens  would  be  required  to 
determine  the  complete  shape  of  the  predicted  curve. 

One  final  observation  to  be  made  from  the  figures  just  discussed  is 
t..e  effect  of  incidence  angle  on  the  attenuation.  It  is  seen  that  the 
measured  attenuation  is  increased  at  higher  values  of  incidence  angle, 
and  the  effect  is  most  prominent  with  the  six-bladed  rotor.  This  result 
was  observed  previously  and  discussed  in  Section  4.3. 

When  dealing  with  inlet  distortions,  two  parameters  are  used  to 
describe  the  unsteadiness  of  the  flow;  these  two  parameters  are  the 
amplitude  and  the  phase  angle.  The  amplitude  has  been  discussed  in  the 
preceding  paragraphs.  The  phase  angle  merely  shows  how  the  flow  is 
displaced  circumferentially  due  to  both  the  steady  and  unsteady  effects 
of  the  rotor.  In  Figures  74  and  75,  the  phase  angle  upstream  of  the  rotor 
is  plotted  as  a  function  of  the  ratio  of  blade  spacing  to  distortion  wave¬ 
length.  The  phase  angle  of  both  the  axial  velocity  and  total  pressure 
distortions  are  presented.  The  data  in  Figure  74  represent  the  six-  and 
twelve-bladed  rotors  having  stagger  angles  of  35  and  45  degrees.  Figure 
75  presents  the  data  for  the  six-  and  twelve-bladed  rotors  with  a  stagger 
angle  of  55  degrees  and  for  the  nine-bladed  cambered  rotor  with  a  stagger 
angle  of  50  degrees.  Since  these  figures  show  the  phase  angle  at  the 


Figure  75.  Phase  Angle  versus  Ratio  of  Blade  Spacing  to 
Distortion  Wavelength  Upstream  of  Rotor  for 
Stagger  Angles  of  50  and  55  Degrees 


inlet  before  its  passage  through  the  rotor,  all  the  data,  as  expected, 
show  a  constant  phase  angle  of  approximately  95  degrees.  This  represents 
the  phase  angle  of  the  distortion  only,  and  it  is  not  influenced  by  the 
action  of  the  rotor.  The  one-,  two-,  four-,  and  six-cycle  distortion 
screens  were  all  placed  in  the  inlet  of  the  AFRF  such  that  a  minimum  value 
of  velocity  occurred  at  a  circumferential  location  of  zero  degrees.  Since 
the  screens  were  all  placed  in  the  AFRF  in  the  same  way,  the  phase  angles 
of  the  fundamental  harmonics  at  the  inlet  to  the  rotor  were  all  approxi¬ 
mately  the  same  and  equal  to  approximately  90  degrees. 

As  the  distortion  is  transported  through  the  rotor,  several  factors 
contribute  to  the  change  of  phase  angle.  Among  these  are  the  steady  lift 
on  the  rotor  blades,  the  magnitude  and  phase  angle  of  the  unsteady  circu¬ 
lation  referenced  to  the  trailing  edge  of  the  blade,  and  the  convection 
of  the  distortion  downstream  of  the  rotor.  These  factors  were  all 
included  in  the  theoretical  analysis  to  give  a  predicted  value  of  the 
phase  angle  of  the  total  pressure  distortion  downstream  of  the  rotor. 

These  predicted  phase  angles  are  plotted  and  compared  to  the  phase  angles 
measured  experimentally  downstream  of  the  rotor  in  the  AFRF  in  Figures 
76  and  77.  Once  again,  the  large  discrepancy  between  the  measured  and  the 
predicted  values  may  be  due  to  the  effects  of  blade  thickness  and  viscosity 
(i.e.,  boundary  layer)  which  are  neglected  in  the  theory.  The  presence 
of  a  boundary  layer  on  the  blade  acts  to  delay  the  response  of  the  blade 
tj  the  disturbance  flow.  The  existence  of  an  unsteady  radial  or  spanwise 
flow  may  also  contribute  to  the  large  difference  between  the  theoretical 
and  experimental  results.  An  increase  in  the  phase  angle  is  noticed  in 
the  experimental  data  at  the  higher  incidence  angles,  and  it  is  most 
noticeable  at  larger  values  of  S/l.  Higher  mean  incidence  angles  result 


Figure  77.  Phase  Angle  versus  Ratio  of  Blade  Spacing  to 
Distortion  Wavelength  Downstream  of  Rotor 
Showing  the  Predicted  Values  of  Total  Pressure 
for  Stagger  Angles  of  50  and  55  Degrees 
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in  greater  steady  lift  on  the  rotor  blades  (up  to  a  point),  thus  increasing 
the  amount  of  phase  shift  in  the  circumferential  direction.  This  effect 
would  also  be  reflected  in  the  theoretical  values  of  phase  angle;  however, 
the  theory  was  only  applied  to  cases  with  an  incidence  angle  of  zero 
degrees. 

The  conclusions  reached  in  this  chapter  are  very  important,  and  they 
provide  a  successful  tool  for  the  designer  which  can  be  applied  to  the 
unsteady  design  of  an  axial  flow  fan  or  compressor.  The  data  of  this  study 
show  that,  if  the  type  of  inlet  distortion  is  known,  the  maximum  amount  of 
attenuation  of  the  distortion  can  be  achieved  by  designing  a  blade  row 
w.th  a  blade  spacing  of  one-half  the  distortion  wavelength.  Both  experi¬ 
mental  and  analytical  efforts,  however,  have  only  reached  this  conclusion 
for  sinusoidal  inlet  distortions.  This  conclusion  was  not  observed  for 
the  90  and  180  degree  square  distortions  since  the  values  of  reduced 
frequency  were  too  low  for  the  square  distortions. 


CHAPTER  VI 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 
FOR  FUTURE  RESEARCH 

6. 1  Summary  and  Conclusions 

A  parametric  study  of  the  effects  of  blade  stagger  angle,  mean  blade 
loading,  rotor-stator  spacing,  blade  spanwise  location,  and  distortion 
characteristics  on  the  flow  fields  upstream  and  downstream  of  two  test 
rotors  has  been  conducted.  The  investigation  involved  the  collection  of 
a  large  amount  of  experimental  data  to  describe  the  effects  of  the  various 
parameters  just  listed.  These  data  can,  hopefully,  be  used  in  future 
studies  to  obtain  a  better  understanding  of  distortion  phenomena  and 
to  improve  current  prediction  techniques.  It  is  also  hoped  that  data 
such  as  the  type  provided  in  this  study  will  aid  turbomachinery  designers 
in  the  design  of  axial  flow  fans  and  compressors  with  improved  perfor¬ 
mance.  A  complete  set  of  the  data  obtained  in  this  investigation  can 
be  found  in  Reference  [39].  The  study  has  also  revealed  some  important 
parameters  affecting  distortion  attenuation  and  unsteady  losses,  which 
should  provide  guidance  for  conducting  future  experimental  and  analytical 
research  in  this  area. 

On  the  basis  of  the  results  obtained  from  the  experimental  and 
analytical  research  conducted  in  this  investigation,  the  following  con¬ 
clusions  can  be  made: 

1.  The  incidence  angle  or  mean  blade  loading  was  found  to 
have  very  little  effect  on  the  distortion  attenuation 
for  the  nine-bladed  cambered  rotor.  This  is  probably 
due  partially  to  the  linearity  of  the  steady  performance 
characteristics  of  the  rotor.  For  the  tests  made  with 
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the  six-bladed  uncambered  rotor,  the  effect  of  incidence 
angle  was  more  significant.  The  attenuation  of  the 
distortion  downstream  of  the  rotor  was  found  to  increase 
as  the  incidence  angle  increased.  This  dependence  on 
blade  loading  is  due  to  the  non-linear  behavior  of  the 
performance  curve  at  high  angles  of  incidence  and  also 
to  greater  turbulent  mixing. 

2.  The  effect  of  rotor-stator  spacing  was  found  to  be 
negligible.  Although  a  blockage  effect  was  seen  to 
cause  a  drastic  change  in  the  velocity  and  static 
pressure  profiles  with  the  stator  row  located  in  the 
forward  position,  the  amplitudes  of  the  axial  velocity 
and  total  pressure  distortions  were  essentially 
unchanged  for  the  rotor-stator  spacings  of  2.0  and 
0.5  blade  chord  lengths  that  were  tested.  This  is 
because  only  the  fundamental  harmonic  of  the  distortion 
was  considered  in  the  calculation  of  the  distortion 
amplitude. 

3.  Although  a  decrease  in  the  distortion  attenuation  was 
observed  downstream  of  the  rotor  as  the  radial  measure¬ 
ment  distance  from  the  hub  was  increased,  this  effect 
was  small.  A  difference  of  only  approximately  ten  to 
fifteen  percent  was  found  between  the  attenuation 
measured  closest  to  the  hub  and  that  measured  closest 
to  the  outer  casing.  The  most  noticeable  effect  of  the 
different  spanwise  locations  occurred  in  the  unsteady 
total  pressure  losses;  the  losses  were  significantly 
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increased  as  the  measurement  locations  approached  the 
surface  of  the  inner  hub.  This  increase  in  losses  near 
the  hub  could  be  due  to  the  growth  of  the  hub  wall 
boundary  layer  and  also  secondary  losses  occurring 
near  the  hub. 

4.  Increasing  the  reduced  frequency  was  found  to  result  in 
an  increased  attenuation  of  the  sinusoidal  distortions 
that  were  tested  using  the  nine-bladed  rotor  as  S/SL 
approached  0.5.  The  one  case  of  S/£  >  0.5  showed  the 
attenuation  to  decrease  from  the  maximum  value. 

Increasing  the  blade  stagger  angle  for  a  fixed  value 

of  reduced  frequency  with  the  six-bladed  rotor  also 
resulted  in  an  increase  in  the  distortion  attenuation. 

It  should  be  noted  that  the  stagger  angle  is  contained 
in  the  definition  of  reduced  frequency,  and  an  increase 
in  E,  will  increase  u>.  Stagger  angle  also  changes  the 
disturbance  velocities  that  are  parallel  and  perpen¬ 
dicular  to  the  axial  direction  (i.e.,  the  chordwise 
and  transverse  gusts).  The  stagger  angle,  therefore, 
appears  to  be  a  very  important  parameter  influencing 
the  attenuation  of  a  distortion. 

5.  The  result  obtained  in  References  [37]  and  [40]  that 
the  unsteady  losses  reached  a  maximum  at  low  values  of 
reduced  frequency  is  not  observed  in  this  investigation. 

In  most  cases,  the  operation  of  a  rotor  in  a  distorted 
flow  field  was  found  to  result  in  total  pressure  losses 
greater  than  those  measured  with  the  rotor  operating 
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in  a  uniform  flow.  However,  no  peak  values  of  unsteady 
losses  occurred  at  low  values  of  us  for  the  cambered 
rotor.  For  the  uncambered  rotor,  such  a  peak  may  exist, 
but  experimental  measurements  at  low  values  of  reduced 
frequency  are  needed  to  confirm  this  suspicion. 

An  important  conclusion  drawn  from  the  data  with  six- 
and  nine-bladed  rotors  employed  in  this  study  and  with 
the  twelve-bladed  rotor  used  in  Reference  [37]  is  that 
the  inlet  distortion  is  most  greatly  attenuated  at  a 
value  of  rotor  blade  spacing  equal  to  one-half  the 
distortion  wavelength.  As  the  ratio  of  blade  spacing 
to  distortion  wavelength  is  further  increased  and 
approaches  one  (i.e.,  the  number  of  blades  equals  the 
number  of  distortion  cycles),  the  attenuation  is 
decreased.  This  explains  the  decrease  in  attenuation 
for  increasing  values  of  reduced  frequency  which  was 
observed  with  the  six-bladed  rotor;  a  decrease  in 
attenuation  was  observed  for  the  case  of  a  six-bladed 
rotor  in  a  six-cycle  distortion.  This  conclusion  is 
particularly  important  since  it  shows  that,  knowing 
the  shape  of  the  inlet  distortion,  the  maximum  attenu¬ 
ation  can  be  achieved  by  designing  a  rotor  such  that 
the  blade  spacing  is  one-half  the  distortion  wave¬ 
length.  This  conclusion  also  appears  to  be  supported 
by  the  unsteady  theory  used  in  this  investigation  with 
one  exception;  no  decrease  in  attenuation  is  observed 
for  the  uncambered  rotor  as  the  ratio  of  spacing  to 
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to  wavelength  approaches  one.  In  addition,  the  measured 
attenuation  was  found  to  be  much  greater  than  the  atten¬ 
uation  predicted  by  the  theoretical  analysis.  This  is 
not  surprising,  however,  since  the  theory  does  not 
account  for  the  effect  of  viscosity. 

6. 2  Recommendations  for  Future  Research 

Besides  providing  data  of  a  fundamental  nature  which  demonstrate  the 
effects  of  various  design  and  operating  variables  on  the  response  of  a 
rotor  to  inlet  distortions,  this  investigation  also  uncovered  several  areas 
where  future  research  should  be  conducted.  Some  of  the  future  research 
will  serve  to  verify  or  explain  the  results  of  the  present  study.  Other 
research  should  expand  on  the  experimental  and  theoretical  studies  made 
in  this  investigation,  thus  providing  a  more  complete  and  detailed  under¬ 
standing  of  turboraachinery  unsteady  aerodynamics  and  inlet  distortion 
phenomena . 

A  more  thorough  investigation  should  be  undertaken  to  determine  the 
effects  of  the  radial  flows  in  the  AFRF.  It  has  been  pointed  out  that 
radial  shifts  in  the  streamlines  will  alter  the  distortion  amplitude  at 
the  mean  radius  and  can  yield  misleading  results  with  regard  to  the  dis¬ 
tortion  attenuation  and  the  unsteady  losses. 

Another  recommendation  for  future  research  is  to  conduct  a  more 
extensive  experimental  program  employing  a  larger  number  of  distortion 
screens.  This  would  give  a  larger  number  of  reduced  frequencies  and 
intrablade  frequencies,  thus  allowing  both  experimental  and  theoretical 
attenuation  data  and  unsteady  total  pressure  losses  to  be  more  accurately 
analyzed.  Further  experimental  studies  should  be  conducted  with  the  nine- 
bladed  cambered  rotor  in  distortions  with  higher  values  of  reduced 
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frequency  or  intrablade  phase  angle.  This  would  help  to  verify  the  shape 
of  the  curve  of  the  amplitude  ratio  as  a  function  of  blade  spacing  over 
distortion  wavelength  observed  for  the  uncambered  rotor.  Also,  low  values 
of  u  and  x  should  be  tested  experimentally  with  the  six-bladed  rotor  since 
only  the  four-  and  six-cycle  sinusoidal  distortions  were  evaluated  in  the 
present  study. 

A  targe  discrepancy  was  observed  between  the  amount  of  attenuation 
measured  experimentally  and  the  amount  of  attenuation  predicted  by  a  two- 
dimensional,  unsteady  theoretical  analysis.  It  was  suggested  that  the 
discrepancy  was  caused  by  the  fact  that  the  effects  of  viscosity  and 
turbulent  mixing  were  neglected  in  the  theory.  Further  study  should  be 
aimed  at  discovering  the  causes  of  the  attenuation  observed  in  the  measured 
data  and  also  at  providing  an  analytical  method  which  includes  viscous 
effects  in  the  prediction  of  the  total  pressure  distribution  downstream 
of  the  rotor.  Mokelke  in  Reference  [30]  has  analytically  investigated  the 
contribution  of  turbulent  mixing,  and  a  similar  analysis  could  be  used 
as  a  starting  point  for  a  study  of  the  observed  phenomena  in  the  AFRF. 

From  Reference  [37],  it  is  known  that  no  significant  turbulent  mixing 
occurs  without  the  rotor  installed;  if  turbulent  mixing  is  significant, 
it  roust  be  caused  by  the  rotor. 

The  last  topic  to  be  included  in  recommendations  for  future  research 
is  the  area  of  unsteady  losses.  A  maximum  value  of  unsteady  losses  at 
low  values  of  reduced  frequency  was  observed  in  References  [37]  and  [40]. 

To  verify  such  a  peak  for  the  six-bladed  test  rotor  used  in  this  study, 
measurements  at  lower  values  of  reduced  frequency  are  required.  The 
effect  of  radial  measurement  location  on  the  unsteady  losses  was  found  to 
be  very  important;  however,  no  definitive  explanation  of  the  observed 
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behavior  is  available.  Possible  sources  of  unsteady  losses  are  local 
regions  of  flow  separation,  turbulent  mixing,  vortex  shedding  due  to  the 
unsteady  blade  circulation,  pressure  and  viscous  drag,  and  hub  and  annulus 
wall  boundary  layers.  The  contributions  of  these  sources  to  the  losses 


should  be  investigated. 
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APPENDIX  A 


RESOLUTION  OF  THE  VELOCITY  COMPONENTS  FOR  THE  FIVE-HOLE  PROBES 

In  this  appendix,  an  analytical  development  of  the  vector  resolution 
of  the  three  velocity  components  measured  by  the  five-hole  probes  is 
presented.  It  is  desired  to  know  the  axial,  circumferential,  and  radial 
components  of  a  velocity  vector  V  at  angles  y  and  y  with  respect  to  the 
local  probe  axes.  Actually,  these  components  depend  on  the  definition  of 
the  pitch  angle  y  and  the  yaw  angle  y.  Calibration  of  the  probe  where 
the  first  rotation  was  in  the  yaw  plane  followed  by  rotation  in  the  pitch 
plane  (THE  YAW-PITCH  MODE)  yield  one  set  of  values  for  the  velocity 
components.  However,  the  alternate  rotational  procedure,  that  is,  the 
pitch-yaw  mode,  will  result  in  distinctly  different  values  for  the  velocity 
components.  The  particular  five-hole  probes  used  in  this  study  were 
calibrated  in  the  yaw-pitch  mode.  Thus,  only  the  vector  resolution  for 
probes  calibrated  in  the  yaw-pitch  mode  will  be  presented.  This  derivation 
of  the  velocity  components  exactly  follows  that  presented  in  Reference 
[36],  and  is  reproduced  here  for  convenience. 

Consider  any  vector  V  whose  components  are  V^,  V^,  and  V^  in  the 
original  or  unprimed  coordinate  system  and  V^,  ,  and  V in  the  primed 

or  rotated  coordinate  system  (see  Figure  A.l).  It  can  be  shown  that  the 
primed  components  can  be  computed  from 


V  '  J  V*  • 


where  a^'s  are  direction  cosines  defined  by: 


ij 


(A.l) 


cos  (Xi'»  Xj)  . 


(A. 2) 
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The  symbol  (X^',  X^)  represents  the  angle  between  the  positive  direction 
of  the  axes  X^'  and  X^ . 


Figure  A-l.  Rotated  Coordinate  System 


Let  V  be  the  reference  velocity  in  the  positive  X^  direction  of  the 
original  coordinate  system;  that  is,  V.^  -  0.0,  V 2  =  0.0,  and  *  V.  The 
first  rotation  is  a  yawing  motion  about  the  X^  axis  of  magnitude  y.  The 
geometry  is  illustrated  in  Figure  A.l. 


Thus 
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all 

as 

COS 

(X1% 

V 

« 

cos 

(0°)  =  1.0, 

a21 

= 

cos 

(X2', 

V 

= 

cos 

(90°)  =  0.0, 

a31 

= 

cos 

(X3', 

V 

= 

cos 

(90°)  =  0.0, 

a12 

= 

cos 

(X1% 

X2> 

= 

cos 

(90°)  =  0.0, 

a22 

= 

cos 

(x/. 

x2) 

= 

cos 

(y)  =  cos  y , 

a32 

= 

cos 

(x3'. 

x2) 

= 

cos 

(90®+y)  =  -siny 

al3 

= 

cos 

( x 1', 

X3) 

= 

cos 

(90°)  =  0.0, 

a23 

= 

cos 

(x2\ 

V 

= 

cos 

(90*— y)  =  siny. 

a33 

= 

cos 

(x3'» 

x3) 

= 

cosy 

,  , 

(A.  3) 


and 

Fr  m  Equation  (A.l),  the  components  of  V  in  the  rotated  axis  system  are: 

V  ’  *11  V1  +  “12  V2  +  “13  V3  ’  °-0’  (A‘4> 


and 


V2  =  a21  V1  +  a22  V2  +  a23  V3  =  V  sin  y’ 


(A.  5) 


V3  =  a31  V1  +  a32  V2  +  a33  V3  V  sin 


(A.  6) 


For  the  yaw-pitch  calibration,  the  next  rotation  is  a  pitching  motion 
about  the  X 2  axis  of  magnitude  a.  This  final  axis  orientation  will  be 
denoted  by  the  double  primed  symbols.  The  geometry  is  shown  in  Figure 
A. 2. 
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Figure  A-2.  Final  Coordinate  System 

all  =  cos  (Xi"’  xi  ^  =  cos  Y» 

a21  =  COS  ^X2 xi =  cos  (90®)  =  0.0, 

a31  =  cos  &  ",  X  ')  «  cos  (90*+y)  =  -siny  , 

a12  -  cos  (X1",  X2')  «  cos  (90°)  -  0.0, 

a22  =  C0S  (X2 X2  ^  =  cos  “  1*0»  (A. 7) 

a32  =  cos  (X  X  ')  =  cos  (90®)  =  0.0, 

a13  =  cos  (Xi~'  xj  )  *  cos  (90®-y)  =  siny, 

a23  -  cos  (X2",  X3')  -  cos  (90®)  -  0.0, 

and 

a33  *  cos  (X3  "•  X3 ”)  ■  cosy* 

The  components  of  V,"",  Vn",  and  V  ""  of  the  original  vector  V 

J-  A.  J 

relative  to  the  final  axis  orientation  are  obtained  by  combining  Equations 
A. 4  through  A. 7,  using  Equation  A.l. 


(A. 8) 


V'  =  an  Y  +  a!2  V  +  al3  V  =  Vcos^sinY» 

y 2  '  =  a21  V  +  a22  ^ 2  +  a23  V  =  ?s±n}J’  (A'9) 

and 

V3"  a^  V1"  +  a^2  V2'  +  =  VcosycosY*  (A.  10) 

In  many  applications,  these  three  velocity  components  are  identified  as 
follows: 


Radial:  =  V  ""  ~  VcosusinY*  (A. 11) 

Circumferential:  V~  =  V."  =  Vsiny.  (A.  12) 

0  L 

Axial:  =  ^2"  =  VcosucosY*  (A. 13) 

For  surface  ship  applications,  V  is  positive  when  directed  radially 

inward  along  the  longitudinal  probe  axis;  VQ  is  positive  when  directed 

w 

counterclockwise  when  looking  upstream;  and  V  is  positive  in  the  axial 

A 

direction  as  shown  in  Figure  A. 3. 

In  some  cases  where  the  rotor  was  operated  at  large  blade  loadings, 
there  were  large  components  of  circumferential  velocity,  or  swirl,  present 
downstream  of  the  rotor.  In  these  cases,  the  values  of  the  downstream 
yaw  angle  sensed  by  the  probe  were  large  and  fell  outside  the  calibration 
range  of  ±30  degrees.  Since  the  velocity  components  were  calculated  using 
the  probe  calibration  data,  the  velocities  calculated  for  these  particular 
cases  were  meaningless.  Thus,  it  was  necessary  to  rotate  the  probe  in  the 
counterclockwise  direction  by  an  angle  a  in  order  that  the  probe  would 


sense  yaw  angle  values  within  the  calibration  range.  When  the  probe  is 
rotated  in  this  manner,  the  equations  for  the  yaw  and  pitch  angles  used 
in  the  data  reduction  program  must  be  changed.  The  derivation  of  these 
new  angles  is  presented  here. 
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Referring  to  Figure  A. 3,  the  angles  y  and  y  are  the  pitch  and  yaw 
angles  relative  to  the  probe;  PA  and  YA  are  the  pitch  and  yaw  angles 
relative  to  the  axial  reference,  and  they  are  the  angles  used  in  the  data 
reduction  program,  a  is  the  angle  that  the  probe  is  rotated. 

The  three  velocity  components  are  unchanged: 

Vx  =  VcosYAcosPA,  (A. 11) 

VQ  =  VsinYA,  (A. 12) 

U 

and 

V  =  VcosYAsinPA.  (A. 13) 

r 

Now,  with  the  probe  rotated,  we  have: 


and 


Vx  =  V'  cos  (YA-y)  cosPA, 

V  =  V'cosYcos a, 

A. 

V  =  V'siny, 
r 

=  V'cos  (YA-y)  sinPA, 

cos  (YA-y)  cosPA  =  cosycoso, 


siny  =  cos  (YA-y)  sinPA  . 


(A. 14) 

(A. 15) 

(A. 16) 

(A. 17) 

(A.  18) 

(A. 19) 


Substitution  of  Equation  (A. 19)  into  Equation  (A. 18)  gives: 
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COSPA  =  COSYCO80  . 


Then,  by  rearranging,  we  have: 


and 


cosPA  _  cosy 
sinPA  sinY 


co  so 


tanPA  = 

sxna 


Therefore, 


PA 


tan 


-1 


lcosCT  ' 


(A. 20) 


From  equation  (A. 18), 


cos  (YA-y) 


cosYcosQ 

cosPA 


Therefore, 


YA  =  cos 


'cosPA 


cosa)  +  y  . 


(A. 21) 


Equations  (A. 20)  and  (A. 21)  are  the  new  angles  which  are  used  in  the  cases 
where  the  downstream  probe  is  rotated  counterclockwise  in  the  yaw  plane. 


APPENDIX  B 


SUMMARY  OF  EXPERIMENTAL  TEST  CONDITIONS 


Table  B.l.  Table  of  Run  Conditions  Completed 
for  the  Six-Bladed  Rotor 


Stagger 

Av,  V* 

AP 

Screen 

(Deg.) 

RPM 

Airrr 

Avg.  Incidence 

m 

AVg  • 

l/2pVx2avg. 

(Deg.) 

35 

682 

1.303 

0.023 

2.29 

35 

941 

0.992 

0.454 

10.12 

4-Cycle 

35 

1203 

0.802 

1.004 

16.24 

45 

975 

0.985 

0.014 

0.40 

Sinusoidal 

Distortion 

45 

1293 

0.739 

0.660 

8.50 

55 

1392 

0.675 

0.093 

0.57 

55 

1622 

0.526 

0.848 

6.78 

55 

1913 

0.493 

1.732 

8.35 

35 

714 

1.354 

0.001 

1.20 

35 

1020 

1.010 

0.424 

9.58 

6-Cycle 

35 

1260 

0.814 

0.911 

15.79 

45 

1020 

0.971 

0.029 

0.76 

Sinusoidal 

45 

1413 

0.696 

0.851 

10.13 

55 

1457 

0.679 

0.131 

0.19 

55 

1698 

0.586 

0.708 

4.25 

55 


2002 


0.487 


1.767 


8.52 
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Table  B.l.  Table  of  Run  Conditions  Completed 

for  the  Six-Bladed  Rotor  (Continued) 


Screen 

Stagger 
(Deg . ) 

RPM 

vx 

Avg.  Incidence 
(Deg.) 

Av8-  u 

m 

A 

AVg  • 

l/2pV  2 

X  avg. 

35 

696 

1.363 

0.009 

1.29 

35 

993 

0.963 

0.479 

11.07 

90  Degree 

35 

1143 

0.817 

0.847 

15.75 

45 

994 

0.990 

0.027 

0.34 

Square 

Distortion 

45 

1375 

0.722 

0.800 

9.18 

55 

1419 

0.615 

0.145 

3.02 

55 

1653 

0.552 

0.785 

5.57 

55 

1950 

0.483 

1.787 

8.84 

35 

692 

1.396 

0.017 

0.80 

35 

955 

0.977 

0.435 

10.78 

180  Degree 

35 

1138 

0.866 

0.766 

14.18 

45 

989 

0.951 

0.054 

1.57 

Square 

Distortion 

45 

1328 

0.717 

0.755 

9.41 

55 

1412 

0.634 

0.148 

2.18 

55 

1646 

0.533 

0.801 

6.47 

55 


1941 


0.474 


1.794 


9.35 
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Table  B.2.  Table  of  Run  Conditions  Completed  for  the 
Nine-Bladed  Rotor  With  a  Rotor-Stator 
Spacing  of  12  Inches  (30.48  cm) 


Screen 

Stagger 
(Deg. ) 

RPM 

vx 

Avg.  r 
m 

AP 

T 

Avg.  Incidence 
(Deg.) 

AVg. 

l/2pV  2 

X  avg. 

1-Cycle 

50 

1525 

0.647 

1.104 

0.65 

Sinusoidal 

50 

1699 

0.583 

1.775 

3.36 

Distortion 

50 

1837 

0.532 

2.444 

5.65 

2-  ’ycle 

50 

1529 

0.638 

1.188 

1.18 

Sinusoidal 

50 

1704 

0.564 

1.960 

4.28 

Distortion 

50 

1842 

0.517 

2.623 

6.40 

4-Cycle 

50 

1492 

0.620 

1.281 

1.88 

Sinusoidal 

50 

1663 

0.536 

2.248 

5.47 

Distortion 

50 

1798 

0.492 

2.951 

7.50 

6-Cycle 

50 

1462 

0.626 

1.244 

1.65 

Sinusoidal 

50 

1741 

0.527 

2.316 

5.97 

Distortion 

50 

1882 

0.530 

2.515 

5.81 

90  Degree 

50 

1521 

0.616 

1.260 

2.08 

Square 

50 

1695 

0.574 

1.800 

3.84 

Distortion 

50 

1833 

0.523 

2.494 

6.09 

180  Degree 

50 

1515 

0.605 

1.295 

2.56 

Square 

50 

1688 

0.553 

2.011 

4.79 

Distortion 

50 

1825 

0.507 

2.679 

6.85 

Table  B.3.  Table  of  Run  Conditions  Completed  for  the 
Nine-Bladed  Rotor  With  a  Rotor-Stator 
Spacing  of  3  Inches  (7.62  cm) 


Screen 

Stagger 

(Deg.) 

RPM 

Avg-  o* 

m 

apt 

Avg. 

Avg. 

l/2pV  2 

X  avg. 

1-Cycle 

50 

1525 

0.644 

1.117 

Sinusoidal 

50 

1699 

0.577 

1.825 

Distortion 

50 

1837 

0.532 

2.364 

2-Cycle 

50 

1529 

0.632 

1.210 

Sinusoidal 

50 

1704 

0.565 

1.910 

Distortion 

50 

1842 

0.514 

2.639 

Ur-  Cycle 

50 

1492 

0.610 

1.356 

Sinusoidal 

50 

1663 

0.534 

2.252 

Distortion 

50 

1798 

0.489 

2.968 

6-Cycle 

50 

1562 

0.618 

1.301 

Sinusoidal 

50 

1741 

0.536 

2.133 

Distortion 

50 

1882 

0.541 

2.305 

90  Degree 

50 

1521 

0.604 

1.360 

Square 

50 

1696 

0.573 

1.796 

Distortion 

50 

1833 

0.522 

2.482 

180  Degree 

50 

1515 

0.606 

1.268 

Square 

50 

1688 

0.579 

1.762 

Distortion 

50 

1825 

0.530 

2.356 

Incidence 

(Deg.) 

0.95 

3.73 
5.72 
1.56 
4.37 
6.64 
2.41 

4.70 

7.74 
2.09 
5.61 
5.40 

2.70 
3.99 
6.28 
2.69 

3.75 
4.94 
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Table  B.4.  Table  of  Run  Conditions  Completed  for  the 
Nine-Bladed  Rotor  With  the  Probe  Located 
at  a  Radial  Distance  of  4.583  Inches 
(11.64  cm)  from  the  Hub  Surface* 


Screen 

Stagger 

(Deg.) 

RPM 

— 

A  Vx 

A  .  APT 

Avg.  Incidence 
(Deg.) 

Avg.  u 

m 

AVg . 

l/2pV2 

X  avg. 

2 -Cycle 

50 

1704 

0.537 

2.289 

5.59 

4 “Cycle 

50 

1663 

0.534 

2.313 

5.77 

90  Deg.  Sq. 

50 

1695 

0.531 

2.255 

5.86 

180  Deg.  Sq. 

50 

1688 

0.530 

_ 

2.272 

5.95 

Table  B.5.  Table  of  Run  Conditions  Completed  for  the 
Nine-Bladed  Rotor  With  the  Probe  Located 
at  a  Radial  Distance  of  3.667  Inches 
(9.31  cm)  from  the  Hub  Surface 


. . 

Screen 

Stagger 
(Deg. ) 

RPM 

— 

A  ^ 

Avg.  — 

m 

AP^ 

Aver  ^ 

Avg.  Incidence 
(Deg . ) 

AVg. 

1/2”V»vg. 

2- Cycle 

50 

1704 

0.547 

2.131 

5.22 

4- Cycle 

50 

1663 

0.529 

2.321 

5.97 

90  Deg.  Sq. 

50 

1695 

0.540 

2.118 

5.53 

180  Deg.  Sq. 

50 

1688 

0.583 

1-807 

3.66 

*NOTE:  The  mean  radius  is  2.75  inches  (6.99  cm)  from  the 
Hub  Surface. 
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Table  B.6.  Table  of  Run  Conditions  Completed  for  the 
Nine-Bladed  Rotor  With  the  Probe  Located 
at  a  Radial  Distance  of  1.833  Inches 
(4.66  cm)  from  the  Hub  Surface 


Stagger 
(Deg . ) 

Av«  V* 

Avg  A1?T 

KrW 

Avg'  u 

m 

AVg. 

l/2pV  2 

X  avg. 

AV^t  LllLAUCllCC 

(Deg.) 

2-Cycle 

50 

1704 

0.567 

1.882 

4.22 

4-Cycle 

50 

1663 

0.552 

1.997 

4.87 

90  Deg.  Sq. 

50 

1695 

0.541 

2.062 

5.36 

180  Deg.  Sq. 

50 

1688 

0.528 

2.165 

5.95 

Table  B.7.  Table  of  Run  Conditions  Completed  for  the 
Nine-Bladed  Rotor  With  the  Probe  Located 
at  a  Radial  Distance  of  0.917  Inch 
(2.33  cm)  from  the  Hub  Surface 


2-Cycle 
4-Cycle 
90  Deg.  Sq. 
180  Deg.  Sq. 


Stagger 

(Deg.) 

RPM 

VX 

Avg.  — 

ID 

APp 

Avg. 

1/2»V.vg. 

50 

1704 

0.533 

2.059 

50 

1663 

0.519 

2.198 

50 

1695 

0.555 

1.917 

50 

1688 

0.547 

1.949 

Avg .  Incidence 
(Deg.) 
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